The following code

from scipy.linalg.sparse import solve as spsolve
foo = spsolve(AMatrixT, dist)

prints into stdout

dgstrf info 10200

And if I inspect foo, it’s all NaN. The shape of dist happens to be (10200, 1). I have uploaded the sparse AMatrixT matrix and the dense dist array as .npz files here. What does the output mean, and how can I stop it?

# Name                    Version                   Build  Channel
python                    3.11.4               h955ad1f_0 
scipy                     1.12.0          py311h08b1b3b_0