Gael Guennebaud | 7579360 | 2009-07-27 18:50:39 +0200 | [diff] [blame] | 1 | #warning deprecated |
| 2 | /* |
Benoit Jacob | 9b06e07 | 2009-01-27 17:42:04 +0000 | [diff] [blame] | 3 | Matrix3d m = Matrix3d::Zero(); |
Benoit Jacob | 9e00d94 | 2008-12-20 13:36:12 +0000 | [diff] [blame] | 4 | m.part<Eigen::StrictlyUpperTriangular>().setOnes(); |
Benoit Jacob | 486fdb2 | 2008-05-29 03:12:30 +0000 | [diff] [blame] | 5 | cout << "Here is the matrix m:" << endl << m << endl; |
| 6 | cout << "And let us now compute m*m.adjoint() in a very optimized way" << endl |
| 7 | << "taking advantage of the symmetry." << endl; |
| 8 | Matrix3d n; |
| 9 | n.part<Eigen::SelfAdjoint>() = (m*m.adjoint()).lazy(); |
| 10 | cout << "The result is:" << endl << n << endl; |
Gael Guennebaud | 7579360 | 2009-07-27 18:50:39 +0200 | [diff] [blame] | 11 | */ |