Intel® Fortran Compiler
Build applications that can scale for the future with optimized code designed for Intel® Xeon® and compatible processors.
29234 Discussions

invert SPD matrix

v-nijs
Beginner
‎08-01-2006 11:02 PM
286 Views
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
0 Kudos

Link Copied

0 Replies
Reply

Community support is provided Monday to Friday. Other contact methods are available here.

Intel does not verify all solutions, including but not limited to any file transfers that may appear in this community. Accordingly, Intel disclaims all express and implied warranties, including without limitation, the implied warranties of merchantability, fitness for a particular purpose, and non-infringement, as well as any warranty arising from course of performance, course of dealing, or usage in trade.

For more complete information about compiler optimizations, see our Optimization Notice.