Benoit Jacob | 183bf54 | 2008-01-13 23:38:48 +0000 | [diff] [blame] | 1 | #include <Eigen/Core> |
| 2 | USING_PART_OF_NAMESPACE_EIGEN |
| 3 | using namespace std; |
| 4 | |
| 5 | template<typename Scalar, typename Derived> |
Benoit Jacob | f12e9c5 | 2008-02-29 13:56:40 +0000 | [diff] [blame^] | 6 | Eigen::Block<Derived, 2, 2> |
Benoit Jacob | 183bf54 | 2008-01-13 23:38:48 +0000 | [diff] [blame] | 7 | topLeft2x2Corner(MatrixBase<Scalar, Derived>& m) |
| 8 | { |
Benoit Jacob | f12e9c5 | 2008-02-29 13:56:40 +0000 | [diff] [blame^] | 9 | return Eigen::Block<Derived, 2, 2>(m.ref(), 0, 0); |
Benoit Jacob | 183bf54 | 2008-01-13 23:38:48 +0000 | [diff] [blame] | 10 | } |
| 11 | |
| 12 | template<typename Scalar, typename Derived> |
Benoit Jacob | f12e9c5 | 2008-02-29 13:56:40 +0000 | [diff] [blame^] | 13 | const Eigen::Block<Derived, 2, 2> |
Benoit Jacob | 183bf54 | 2008-01-13 23:38:48 +0000 | [diff] [blame] | 14 | topLeft2x2Corner(const MatrixBase<Scalar, Derived>& m) |
| 15 | { |
Benoit Jacob | f12e9c5 | 2008-02-29 13:56:40 +0000 | [diff] [blame^] | 16 | return Eigen::Block<Derived, 2, 2>(m.ref(), 0, 0); |
Benoit Jacob | 183bf54 | 2008-01-13 23:38:48 +0000 | [diff] [blame] | 17 | } |
| 18 | |
| 19 | int main(int, char**) |
| 20 | { |
| 21 | Matrix3d m = Matrix3d::identity(); |
| 22 | cout << topLeft2x2Corner(4*m) << endl; // calls the const version |
| 23 | topLeft2x2Corner(m) *= 2; // calls the non-const version |
| 24 | cout << "Now the matrix m is:" << endl << m << endl; |
| 25 | return 0; |
| 26 | } |