blob: a6025231edace4bfd8dc55524f18cb2093e5c8aa [file] [log] [blame]
Benoit Jacob183bf542008-01-13 23:38:48 +00001#include <Eigen/Core>
2USING_PART_OF_NAMESPACE_EIGEN
3using namespace std;
4
5template<typename Scalar, typename Derived>
Benoit Jacobf12e9c52008-02-29 13:56:40 +00006Eigen::Block<Derived, 2, 2>
Benoit Jacob183bf542008-01-13 23:38:48 +00007topLeft2x2Corner(MatrixBase<Scalar, Derived>& m)
8{
Benoit Jacobf12e9c52008-02-29 13:56:40 +00009 return Eigen::Block<Derived, 2, 2>(m.ref(), 0, 0);
Benoit Jacob183bf542008-01-13 23:38:48 +000010}
11
12template<typename Scalar, typename Derived>
Benoit Jacobf12e9c52008-02-29 13:56:40 +000013const Eigen::Block<Derived, 2, 2>
Benoit Jacob183bf542008-01-13 23:38:48 +000014topLeft2x2Corner(const MatrixBase<Scalar, Derived>& m)
15{
Benoit Jacobf12e9c52008-02-29 13:56:40 +000016 return Eigen::Block<Derived, 2, 2>(m.ref(), 0, 0);
Benoit Jacob183bf542008-01-13 23:38:48 +000017}
18
19int 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}