Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.
Announcements
Welcome to the Intel Community. If you get an answer you like, please mark it as an Accepted Solution to help others. Thank you!

Incomplete Cholesky decomposition

Robin_T_
Novice
279 Views

Hi,

my question relates to this previous thread: https://software.intel.com/en-us/forums/intel-math-kernel-library/topic/550882

I am trying to solve a large (N > 50e6) sparse (~1e-6) symmetric linear system of equations.

With PARDISO I am running into memory problems with a 576 GB RAM machine, even when utilizing the OOC capability.

So, I am trying to implement the RCI CG method. I need a suitable preconditioner because I am very quickly reaching large numbers of iterations with the CG method. How can I implement an incomplete Cholesky factorization? MKL only provides LU factorization apparently, which could be used in conjunction with GMRES. However, then I would have to provide the whole matrix (not just one triangle) in CSR format. It might be easier to apply incomplete Cholesky factorization. Is there any way to do so?

Thank you for any advice.

1 Reply
Royi
Novice
272 Views

I also wait for Incomplete Cholesky Decomposition support in Intel MK.

I hope they will support both Thresholding and Zero Filling with Modified Incomplete Cholesky Decomposition.

 

I also expects it to use modern techniques so it will be parallel and the fastest algorithms out there.

Reply