[SciPy-user] solve for symetric matrix ?
Thu Apr 26 10:43:02 CDT 2007
On 26/04/07, Nils Wagner <email@example.com> wrote:
> Strange. I have increased the number of unknowns to 5000. What happens
> with the residual ??
> Elapsed time LU 28.8819189072
> Residual 7.64156112039e-14
> Elapsed time Cholesky 15.6276760101
> Residual 141.421356237
> >>> n
> Any idea ?
Well, roundoff error is always a problem in numerical linear algebra.
Perhaps use the svd to check your condition number?
More information about the SciPy-user