LU factorization of a SPD matrx

Post here if you have a question about LAPACK or ScaLAPACK algorithm or data format

LU factorization of a SPD matrx

Postby gsamaras » Fri Apr 24, 2015 9:24 am

I have some code (given to me) that that uses DGETRF for computing an LU factorization of a general M*N matrix A. Then I use DGETRI to compute the inverse of A.

The problem is that A is huge, making this operation the bottleneck of the program. For that reason, I am looking for alterntives, in order to do better. I see that A matrix is symmetric positive definite (SPD) and not a general matrix, so could I take advantage of this property somehow?

For example, I saw a function that performs Cholesky factorization of a SPD matrix. However, I did not find any such function for LU. Is there any?

Now, I am refreshing my linear algebra knowledge on whether I can discard LU and go for Cholesky.
gsamaras
 
Posts: 1
Joined: Fri Apr 24, 2015 9:09 am

Return to Algorithm / Data

Who is online

Users browsing this forum: No registered users and 2 guests