[SciPy-user] solve for symetric matrix ?

Nils Wagner nwagner@iam.uni-stuttgart...
Thu Apr 26 06:31:53 CDT 2007


Joachim Dahl wrote:
>
>
> On 4/26/07, *fred* <fredmfp@gmail.com <mailto:fredmfp@gmail.com>> wrote:
>
>     Hi,
>
>     I use scipy.linalg.solve() to solve Au=b, where A is symetric (and
>     diagonal filled with 0).
>
>
> If the diagonal is 0, then the matrix is not definite (so, neither is
> -A)...
>
>
>
I missed the main information (zero diagonal)
Are you aware of a simple proof that a symmetric matrix with zero
diagonal entries is indefinite ?

Anyway it would be great if scipy has a solver for symmetric indefinite
systems e.g. the LAPACK routine
http://www.netlib.org/lapack/double/dsysv.f

AFAIK, this subroutine is currently not available via linalg.flapack.

Nils

> ------------------------------------------------------------------------
>
> _______________________________________________
> SciPy-user mailing list
> SciPy-user@scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>   


 



More information about the SciPy-user mailing list