russell_sparse 1.15.0

Solvers for large sparse linear systems (wraps MUMPS and UMFPACK)
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14

import numpy as np
import scipy.linalg as la

a = np.array([
    [ 2, 3, 0, 0, 0],
    [ 3, 0, 4, 0, 6],
    [ 0,-1,-3, 2, 0],
    [ 0, 0, 1, 0, 0],
    [ 0, 4, 2, 0, 1]
], dtype=float)

print(a)
print(f"determinant(a) = {la.det(a)}")