Krylov subspace methods package for Python
Python Makefile
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Failed to load latest commit information.
docs docs: update bibliography and intro Jan 18, 2018
krypy Bump version to 2.1.7 Jan 18, 2018
.gitignore use rst as longdesc Oct 17, 2015
.travis.yml travis: install dependencies via pip and cache them Jan 17, 2018
Makefile use rst as longdesc Oct 17, 2015 add instructions for releasing to README Oct 17, 2015
README.txt add basic package structure Oct 26, 2012 [doc] more debian packaging Jul 24, 2013
example.png [readme] update example Mar 14, 2014
requirements.txt remove mocks Oct 18, 2015 Bump version to 2.1.7 Jan 18, 2018


Build Status Documentation Status doi Pypi version Pypi downloads

KryPy is a Python (versions 2 and 3) module for Krylov subspace methods for the solution of linear algebraic systems. This includes enhanced versions of CG, MINRES and GMRES as well as methods for the efficient solution of sequences of linear systems.


KryPy gives you an easy-to-use yet flexible interface to Krylov subspace methods for linear algebraic systems. Compared to the implementations in SciPy (or MATLAB), KryPy allows you to supply additional arguments that may help you to tune the solver for the specific problem you want to solve. The additional arguments may also be of interest if you are doing research on Krylov subspace methods.

Some features of KryPy are:

  • User-defined inner products - useful when solving a linear algebraic system whose operator is self-adjoined in a non-Euclidean inner-product. This way, CG or MINRES can be applied to self-adjoined (but non-symmetric/non-Hermitian) operators easily.
  • Full control of preconditioners - the order of applying preconditioners matters. This is why you can supply two left preconditioners (one of whom implicitly changes the inner product and thus has to be positive definite) and one right preconditioner. Take a look at the arguments M, Ml and Mr.
  • Get the Arnoldi/Lanczos basis and Hessenberg matrix - you want to extract further information from the generated vectors (e.g. recycling)? Just pass the optional argument store_arnoldi=True.
  • Explicitly computed residuals on demand - if you do research on Krylov subspace methods or preconditioners, then you sometimes want to know the explicitly computed residual in each iteration (in contrast to an updated residual which can be obtained implicitly). Then you should pass the optional argument explicit_residual=True.
  • Compute errors - if you have (for research purposes) the exact solution at hand and want to monitor the error in each iteration instead of the residual, you can supply the optional argument exact_solution=x_exact to the LinearSystem.



The documentation is hosted at


The above convergence history is obtained with the following example where the Gmres method is used to solve the linear system A*x=b with the diagonal matrix A=diag(1e-3,2,...,100) and right hand side b=[1,...,1].

import numpy
from krypy.linsys import LinearSystem, Gmres

# create linear system and solve
linear_system = LinearSystem(A=numpy.diag([1e-3]+range(2, 101)),
                             b=numpy.ones((100, 1)))
sol = Gmres(linear_system)

# plot residuals
from matplotlib import pyplot

Of course, this is just a toy example where you would not use GMRES in practice. KryPy can handle arbitrary large matrices - as long as the (hopefully sparse) matrices and the generated basis of the Krylov subspace fit into your memory. ;) Furthermore, in actual applications, you definitely want to adjust Gmres' parameters such as the residual tolerance.


Help can be optained via Python's builtin help system. For example, you can use the ? in ipython:

from krypy.linsys import Gmres


pip / PyPi

Simply run pip install krypy.


There's an Ubuntu PPA with packages for Python 2 and Python 3.

Installing from source

KryPy has the following dependencies:

  • NumPy
  • SciPy


KryPy is currently maintained by André Gaul. Feel free to contact André. Please submit feature requests and bugs as github issues.

KryPy is developed with continuous integration. Current status: Build Status


To create a new release

  1. bump the __version__ number,

  2. create a Git tag,

    $ git tag -a v0.3.1
    $ git push --tags


  3. upload to PyPi:

    $ make upload


KryPy is free software licensed under the MIT License.


KryPy evolved from the PyNosh package (Python framework for nonlinear Schrödinger equations; joint work with Nico Schlömer) which was used for experiments in the following publication: