Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.

MKL 7.0 Beta PARADISO: Cholesky factors

Intel_C_Intel
Empleados
714 Vistas
Is it possible to access cholesky factors of A (= C*C')?
0 kudos
5 Respuestas
Intel_C_Intel
Empleados
714 Vistas

I have gone back to the developers of sparse solver and they would like to understand why you want this particular quantity. Clearly it would be possible to add such functionality, but if they can understand the unlying motivation for this they may be able to address that.

Bruce

Intel_C_Intel
Empleados
714 Vistas
diagonal elements of "C" are needed for computing the determinant of the A matrix (along with some other things).
igino-pitacco
Principiante
714 Vistas

There are many situations where the factors of the Cholesky decompositon are needed. One of them, very important for many number smashers, regards the solution of the generalized eigenvalue problem:

find (L, X) (L real and X vector) such that

K X = L M X

where K,M are symmetric positive definite matrices.

The generalized problem could be transformed into a classical eigenvalue problem, retaining the symmetry and posiveness of the matrix operator A:

A X = 1/L X

where A depends on K,M. Using A Lanczos or a subspace iteration method for this problem (when coming from the generalized one) requires access to the factors of K (or M).

Please note that generalized eigenvalue analysis is a very common task in many engineering sciences.

Intel_C_Intel
Empleados
714 Vistas

This is exactly the information I need. I can't say for certain when I will have more information on this butI will forward this to the developers.

Bruce

mepad_si
Principiante
714 Vistas

has there been any development re. sparse cholesky ?

Thank you.

mepa

Responder