I want to find the inverse of a symmetric positive definite matrix. I am using lapack’s “dpotri” subroutine which uses Cholesky decomposition to find the inverse. However it returns the same matrix with upper or lower part modified to contain the inverse elements. Let us say I choose the UPLO=’U’ thus saving inverse elements to upper triangular part. I believe the lower part still contains the earlier elements.
My question is how to copy the upper part to lower part or vice-versa in the most efficient way? Is there a lapack/Blas routine to do so ? I am programming in C.