krypy
krypy copied to clipboard
andrenarchy
Krylov subspace methods package for Python
Hi, there are some errors in newer versions of NumPy, e.g. 1.24.0 or later. ```python ============================= test session starts ============================== platform linux -- Python 3.10.9, pytest-7.2.1, pluggy-1.0.0 rootdir: /tmp/makepkg/python-krypy/src/krypy-2.2.0 collected...
x=(N,) or (N,1) is different, (N,) is a vector, and (N,1) is an array, when we use matrix-vector multiply for an array with shape (74121,74121), this should be matrix-free. it...
Using Gmres, the right preconditioning is not working properly. The error might be on line 948 of 'linsys.py'. return self.x0 + self.linear_system.Mr * yk Multiplying yk by the right preconditioner...
`setup.cfg` doesn't know about which tools I should have for development, e.g., flake8, black, isort, nose. I think we should re-introduce `requirements.txt` and add these tools. In the workflow file,...
Analogous to the linter, we should add an isort check to the workflow file.
In the deflation module we use svd although we're only interested in the 2-norm of the matrix. We should verify that `numpy.linalg.norm` correctly computes that 2-norm and remove the svd...
Hi, It seems that find_common_dtype returns wrong dtype when passed an operator (maybe a matrix too). Here is the minimal "not working" example: ``` import numpy as np import scipy.sparse.linalg...
We should add an example to the documentation for catching a `ConvergenceError`. A lot of people are unfamiliar with try/except so we can probably help them here a bit. @zimoun...
For harmonic Ritz pairs, the QR decomposition of the extended Hessenberg matrix should be used in order to solve an eigenvalue problem with better conditioning. See algorithm 5.4 in Freund,...