Gael Guennebaud | 16c7b1d | 2009-09-03 11:17:16 +0200 | [diff] [blame] | 1 | Matrix2f m = Matrix2f::Random(); |
| 2 | m = (m + m.adjoint()).eval(); |
Benoit Jacob | e259f71 | 2010-10-19 21:56:26 -0400 | [diff] [blame] | 3 | JacobiRotation<float> J; |
Gael Guennebaud | 16c7b1d | 2009-09-03 11:17:16 +0200 | [diff] [blame] | 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); |
Aaron Franke | 5c22c7a | 2020-03-21 21:30:06 -0400 | [diff] [blame^] | 8 | cout << "Here is the matrix J' * m * J:" << endl << m << endl; |