jan/meanas
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

[eigensolvers] Increase number of lanczos vectors (ncv) based on number of requested eigenvalues

Browse Source
This commit is contained in:
jan 2025-12-10 23:07:28 -08:00
parent c46bed8298
commit be647658d3
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
meanas/eigensolvers.py
View File

@ -135,7 +135,7 @@ def signed_eigensolve(
shifted_operator = operator + spalg.LinearOperator(shape=operator.shape, shifted_operator = operator + spalg.LinearOperator(shape=operator.shape,
matvec=lambda v: shift * v) matvec=lambda v: shift * v)
shifted_eigenvalues, eigenvectors = spalg.eigs(shifted_operator, which='LM', k=how_many, ncv=50) shifted_eigenvalues, eigenvectors = spalg.eigs(shifted_operator, which='LM', k=how_many, ncv=2 * how_many + 50)
eigenvalues = shifted_eigenvalues - shift eigenvalues = shifted_eigenvalues - shift
k = eigenvalues.argsort() k = eigenvalues.argsort()