namespace ietl {The rayleigh function takes five arguments:
template <class MATRIX, class GEN, class IT, class VS>
std::pair<typename vectorspace_traits<VS>::magnitude_type,
typename vectorspace_traits<VS>::vector_type>
rayleigh(const MATRIX& matrix, GEN& gen, const SOLVER& solver, ITER& iter, const VS& vec);
}
solver(matrix,rho,v,y)The matrix stays the same throughout the iteration and is the same as passed to inverse. This can (and should) be used by the solver to optimize the calculation.