solveForConvergence()

template<typename Solver >

template<typename MatrixScalar , typename Vectors >

SRDef syten::Krylov::ClassicImpl::EigenSolver< Solver >::solveForConvergence ( Solver const &  global_solver, IterativeConfig const &  conf, std::pair< DenseTensor< 2, MatrixScalar >, Vectors > &&  mat_v, SRDef  , SRDef &  eval, Vector &  vec  )

inline

operator () solve the eigen problem, compute error to previous solution

Parameters

mat_v	pair of the matrix to solve and an object that can apply the result vector to the Krylov vectors
eval	the eigenvalue is stored here
global_solver	the solver instance
conf	configuration of the global solver
vec	vector type, only used for temporary storage

Returns

the residuum

References syten::abs(), syten::conj(), syten::DenseTensor< rank, Scalar >::dim(), syten::DenseEigensolver::DenseSymEigensolverBase< Type >::evalues, syten::DenseEigensolver::DenseSymEigensolverBase< Type >::evectors, syten::isnan(), std::move(), syten::Krylov::norm(), syten::Krylov::ClassicImpl::EigenSolver< Solver >::prev_coeffs, syten::Krylov::ClassicImpl::EigenSolver< Solver >::prev_norm, syten::real(), and syten::vsplit().

Here is the call graph for this function: