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?



1 post / 0 new
For more complete information about compiler optimizations, see our Optimization Notice.