## "Fossies" - the Fresh Open Source Software Archive

About: Armadillo is a C++ linear algebra library (matrix maths) aiming towards a good balance between speed and ease of use.

skipping to change at line 53 skipping to change at line 53
Col<eT> fac_f; // residual in the Arnoldi factorisation Col<eT> fac_f; // residual in the Arnoldi factorisation
Mat<eT> ritz_vec; // ritz vectors Mat<eT> ritz_vec; // ritz vectors
Col<eT> ritz_est; // last row of ritz_vec Col<eT> ritz_est; // last row of ritz_vec
std::vector<bool> ritz_conv; // indicator of the convergence of ritz values std::vector<bool> ritz_conv; // indicator of the convergence of ritz values
const eT eps; // the machine precision const eT eps; // the machine precision
// eg. ~= 1e-16 for double type // eg. ~= 1e-16 for double type
const eT eps23; // eps^(2/3), used in convergence test const eT eps23; // eps^(2/3), used in convergence test
// tol*eps23 is the absolute tolerance // tol*eps23 is the absolute tolerance
const eT near0; // a very small value, but 1/near0 does not overf low const eT near0; // a very small value, but 1/near0 does not overf low
std::mt19937_64 local_rng; // local random number generator
inline void fill_rand(eT* dest, const uword N, const uword seed_val);
// Arnoldi factorisation starting from step-k // Arnoldi factorisation starting from step-k
inline void factorise_from(uword from_k, uword to_m, const Col<eT>& fk); inline void factorise_from(uword from_k, uword to_m, const Col<eT>& fk);
// Implicitly restarted Arnoldi factorisation // Implicitly restarted Arnoldi factorisation
inline void restart(uword k); inline void restart(uword k);
// Calculate the number of converged Ritz values // Calculate the number of converged Ritz values
inline uword num_converged(eT tol); inline uword num_converged(eT tol);
// Return the adjusted nev for restarting // Return the adjusted nev for restarting
skipping to change at line 93 skipping to change at line 97
inline uword num_iterations() { return niter; } inline uword num_iterations() { return niter; }
//! Returning the number of matrix operations used in the computation. //! Returning the number of matrix operations used in the computation.
inline uword num_operations() { return nmatop; } inline uword num_operations() { return nmatop; }
//! Returning the converged eigenvalues. //! Returning the converged eigenvalues.
inline Col<eT> eigenvalues(); inline Col<eT> eigenvalues();
//! Returning the eigenvectors associated with the converged eigenvalues. //! Returning the eigenvectors associated with the converged eigenvalues.
inline Mat<eT> eigenvectors(uword nvec); inline Mat<eT> eigenvectors(uword nvec);
//! Returning all converged eigenvectors. //! Returning all converged eigenvectors.
inline Mat<eT> eigenvectors() { return eigenvectors(nev); } inline Mat<eT> eigenvectors() { return eigenvectors(nev); }
}; };
} // namespace newarp } // namespace newarp
End of changes. 2 change blocks.
0 lines changed or deleted 5 lines changed or added