naginterfaces.library.matop.real_​symm_​posdef_​inv_​noref

naginterfaces.library.matop.real_symm_posdef_inv_noref(a)

real_symm_posdef_inv_noref calculates the approximate inverse of a real symmetric positive definite matrix, using a Cholesky factorization.

For full information please refer to the NAG Library document for f01ad

https://support.nag.com/numeric/nl/nagdoc_30.3/flhtml/f01/f01adf.html

Parameters

afloat, array-like, shape $(n+1,n)$

The upper triangle of the $n\times n$ positive definite symmetric matrix $A$. The elements of the array below the diagonal need not be set.

Returns

afloat, ndarray, shape $(n+1,n)$

The lower triangle of the inverse matrix $X$ is stored in the elements of the array below the diagonal, in rows $2$ to $n+1$; $x_{ij}$ is stored in $a[i,j-1]$ for $i\ge j$. The upper triangle of the original matrix is unchanged.

Raises

NagValueError

(errno $1$)

The matrix $A$ is not positive definite, possibly due to rounding errors.

(errno $2$)

On entry, $n=⟨value⟩$.

Constraint: $n\ge 0$.

Notes

No equivalent traditional C interface for this routine exists in the NAG Library.

To compute the inverse $X$ of a real symmetric positive definite matrix $A$, real_symm_posdef_inv_noref first computes a Cholesky factorization of $A$ as $A=L{L}^T$, where $L$ is lower triangular. It then computes $L^{-1}$ and finally forms $X$ as the product $L^{-T}{L}^{-1}$.

References

Wilkinson, J H and Reinsch, C, 1971, Handbook for Automatic Computation II, Linear Algebra, Springer–Verlag