invert SPD matrix

invert SPD matrix

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 contribution / 0 nouveau(x)
Reportez-vous à notre Notice d'optimisation pour plus d'informations sur les choix et l'optimisation des performances dans les produits logiciels Intel.