00001 Matrix2f m = Matrix2f::Random(); 00002 m = (m + m.adjoint()).eval(); 00003 JacobiRotation<float> J; 00004 J.makeJacobi(m, 0, 1); 00005 cout << "Here is the matrix m:" << endl << m << endl; 00006 m.applyOnTheLeft(0, 1, J.adjoint()); 00007 m.applyOnTheRight(0, 1, J); 00008 cout << "Here is the matrix J' * m * J:" << endl << m << endl;