2009-07-28 00:50:39 +08:00
|
|
|
#warning deprecated
|
|
|
|
/*
|
2008-07-21 08:34:46 +08:00
|
|
|
Matrix3d m = Matrix3d::Zero();
|
2008-12-20 21:36:12 +08:00
|
|
|
m.part<Eigen::UpperTriangular>().setOnes();
|
2008-05-29 11:12:30 +08:00
|
|
|
cout << "Here is the matrix m:" << endl << m << endl;
|
2008-07-21 08:34:46 +08:00
|
|
|
Matrix3d n = Matrix3d::Ones();
|
2008-12-20 21:36:12 +08:00
|
|
|
n.part<Eigen::LowerTriangular>() *= 2;
|
2008-05-29 11:12:30 +08:00
|
|
|
cout << "Here is the matrix n:" << endl << n << endl;
|
|
|
|
cout << "And now here is m.inverse()*n, taking advantage of the fact that"
|
|
|
|
" m is upper-triangular:" << endl
|
2008-12-20 21:36:12 +08:00
|
|
|
<< m.marked<Eigen::UpperTriangular>().solveTriangular(n);
|
2009-07-28 00:50:39 +08:00
|
|
|
*/
|