invert SPD matrix

invert SPD matrix

Ritratto di v-nijs@kellogg.northwestern.edu

For the following matrix:

3.0 2.0
2.0 3.0

i get an incorrect answer for the inverse matrix when using the routines for a symmetric positive-definite matrix:

call potrf(M)
call potri(M)

The result i get is:

0.6 2.0
-0.4 0.6

When i use

call getrf(M,ipiv)
call getri(M,ipiv)

i do get the correct answer:

0.6 -0.4
-0.4 0.6

Any ideas why this is happening?

Thanks,

Vincent

1 contenuto / 0 new
Per informazioni complete sulle ottimizzazioni del compilatore, consultare l'Avviso sull'ottimizzazione