[SciPy-user] solve for symetric matrix ?

fred fredmfp at gmail.com
Thu Apr 26 06:52:42 EDT 2007


Ronan Perrussel a écrit :
> fred a écrit :
>> But in fact, he was confused. They are negative definite (kriging 
>> matrices).
>> Bad news :-(
> Then you can solve with -A instead of A?
I tried it, of course ;-)

But cho_factor() rejects -A too.
One of the definition of the positive definite (if I understand right) 
is that the eigenvalues must be positive.

So in my case, as A has negative egv, -A has positive ones.
Although -A has positive ones, cho_factor() raise error.

-- 
http://scipy.org/FredericPetit




More information about the SciPy-User mailing list