When solving the full eigenvalue problem, the solve method should store all the eigenvalues and eigenvectors in a convenient format.
We need an iterate_solver method which increases the resolution until the first M modes have converged to a specified tolerance.
Right now we are wasting computational power when we call solver.solve() in a loop over different modes.
When solving the full eigenvalue problem, the solve method should store all the eigenvalues and eigenvectors in a convenient format.
We need an iterate_solver method which increases the resolution until the first M modes have converged to a specified tolerance.
Right now we are wasting computational power when we call solver.solve() in a loop over different modes.