Brian Silverman | 72890c2 | 2015-09-19 14:37:37 -0400 | [diff] [blame] | 1 | Matrix2f m = Matrix2f::Random(); |
2 | m = (m + m.adjoint()).eval(); | ||||
3 | JacobiRotation<float> J; | ||||
4 | J.makeJacobi(m, 0, 1); | ||||
5 | cout << "Here is the matrix m:" << endl << m << endl; | ||||
6 | m.applyOnTheLeft(0, 1, J.adjoint()); | ||||
7 | m.applyOnTheRight(0, 1, J); | ||||
8 | cout << "Here is the matrix J' * m * J:" << endl << m << endl; |