std::size_t cyclic_eigen_jacobi( const Matrix1& A, Matrix2& V, std::valarray<T>& Lambda, std::size_t const max_rot = 80, const T eps = T( 1.0e-10 ) )
{
Lambda.resize( A.row() );
return cyclic_eigen_jacobi( A, V, Lambda.begin(), max_rot, eps );
}
void eigen_hermitian( const Complex_Matrix1& A, Complex_Matrix2& V, std::valarray<T>& Lambda, const T eps = T( 1.0e-10 ) )
{
Lambda.resize( A.row() );
return eigen_hermitian_impl( A, V, Lambda.begin(), eps);
}
std::size_t eigen_jacobi( const Matrix1& A, Matrix2& V, std::valarray<T>& Lambda, const T eps = T( 1.0e-10 ) )
{
Lambda.resize( A.row() );
return eigen_jacobi( A, V, Lambda.begin(), eps);
}