search

NLESC-JCER / Fortran_Davidson

Davidson eigensolver implemented in Fortran
Apache License 2.0
16 stars 4 forks source link

Create a matrix free version of the algorithm #11

Closed felipeZ closed 5 years ago

felipeZ commented 5 years ago

Currently the matrix to be diagonalized and the optional matrix for the general eigenvalue problem are real densed matrices. The interface to call the davidson eigensolver is:

 generalized_eigensolver(mtx, eigenvalues, ritz_vectors, lowest, method, max_iters, &
        tolerance, iters, max_dim_sub, stx)

We can overload the function call to allow the mtx and the optional argument stx to be functions allowing to have a matrix-free version of the algorithm.

felipeZ commented 5 years ago

Profile the matrix-free implementation