A (sparse) matrix solver for python.
Solving Ax = b should be as easy as:
Ainv = Solver(A)
x = Ainv * bIn pymatsolver we provide a number of wrappers to existing numerical packages. Nothing fancy here.
All solvers work with scipy.sparse matricies, and a single or multiple right hand sides using numpy:
- L/U Triangular Solves
- Wrapping of SciPy matrix solvers (direct and indirect)
- Pardiso solvers
- Mumps solvers
Often, there are faster solvers available for your system than the default scipy factorizations available.
pymatsolver provides a consistent interface to both MKL's Pardiso routines and the MUMPS solver package. To
make use of these we use intermediate wrappers for the libraries that must be installed separately.
The Pardiso interface is recommended for Intel processor based systems. The interface is enabled by
the pydiso python package, which can be installed through conda-forge as:
conda install -c conda-forge pydiso
Mumps is available for all platforms. The mumps interface is enabled by installing the python-mumps
wrapper package. This can easily be installed through conda-forge with:
conda install -c conda-forge python-mumps
Code: https://github.com/simpeg/pymatsolver
Tests: https://github.com/simpeg/pymatsolver/actions
Bugs & Issues: https://github.com/simpeg/pymatsolver/issues
License: MIT