1 Matrix2f m
= Matrix2f::Random();
2 m
= (m
+ m
.adjoint()).eval();
3 JacobiRotation
<float> J
;
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
;