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; 9