Nils Wagner a écrit : > If your matrix is also positive definite you may use > linalg.cho_factor > linalg.cho_solve > Hi Nils, Do you mean that using Cholesky decomposition plus solve methods is faster than simply using solve() ? Cheers, -- http://scipy.org/FredericPetit