Jacobi_makeGivens.cpp 236 B

123456
  1. Vector2f v = Vector2f::Random();
  2. JacobiRotation<float> G;
  3. G.makeGivens(v.x(), v.y());
  4. cout << "Here is the vector v:" << endl << v << endl;
  5. v.applyOnTheLeft(0, 1, G.adjoint());
  6. cout << "Here is the vector J' * v:" << endl << v << endl;