Pb with Python solver interface
The pb reported by @lpoirel is that the the following code in not working after a few iteration:
import pypastix as pastix import scipy.sparse as sps import numpy as np # Set the problem for n in range(5, 100): print(n) A = sps.spdiags([np.ones(n)*i for i in [4, -1, -1, -1, -1]], [0, 1, 3, -1, -3], n, n) x0 = np.ones(n) b = A.dot(x0) # Hack to make sure that the mkl is loaded tmp = np.eye(2).dot(np.ones(2)) # Factorize solver = pastix.solver(A, verbose=False, thread_nbr=1) # Solve x = solver.solve(b, x0=x0, check=True) solver.finalize()
The problem is the corruption of the spm structure that is forwarded to PaStiX.
To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information