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

How to export the results of PARDISO numerical decomposition

HUANGK
Novice
618 Views

Hello!

I want to use PARDISO for matrix decomposition, and I want to export the results of the decomposition, for example, I want to export the L matrix after the LLT decomposition, how do I implement this method? If PARDISO can't implement this step, is there any way to do LLT decomposition for a large sparse matrix?

Thank you for your help

0 Kudos
1 Solution
VidyalathaB_Intel
Moderator
583 Views

Hi Huang,


Thanks for reaching out to us.


>>I want to export the L matrix after the LLT decomposition, how do I implement this method?


Could you please confirm if you are asking to know whether we could extract factors from pardiso?

If yes, this feature is already enabled for the distributed cluster sparse solver via cluster_sparse_solver_export functionality. So, as a workaround, we suggest you try it. You can use the cluster sparse solver with a single MPI process.

Reference Link: https://www.intel.com/content/www/us/en/develop/documentation/onemkl-developer-reference-c/top/sparse-solver-routines/parallel-direct-sp-solver-for-clusters-iface/cluster-sparse-solver-export.html

You can also refer to the example cl_solver_export_c.c under C:\Program Files (x86)\Intel\oneAPI\mkl\2022.0.2\examples\examples_cluster_c\c_mpi\cluster_sparse_solver\source directory.

And here is the thread with a similar discussion which you can refer to for more details.

https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-get-L-U-factors/m-p/1265720#M31046


Please let us know if it helps in answering your queries.


Regards,

Vidya.


View solution in original post

0 Kudos
2 Replies
VidyalathaB_Intel
Moderator
584 Views

Hi Huang,


Thanks for reaching out to us.


>>I want to export the L matrix after the LLT decomposition, how do I implement this method?


Could you please confirm if you are asking to know whether we could extract factors from pardiso?

If yes, this feature is already enabled for the distributed cluster sparse solver via cluster_sparse_solver_export functionality. So, as a workaround, we suggest you try it. You can use the cluster sparse solver with a single MPI process.

Reference Link: https://www.intel.com/content/www/us/en/develop/documentation/onemkl-developer-reference-c/top/sparse-solver-routines/parallel-direct-sp-solver-for-clusters-iface/cluster-sparse-solver-export.html

You can also refer to the example cl_solver_export_c.c under C:\Program Files (x86)\Intel\oneAPI\mkl\2022.0.2\examples\examples_cluster_c\c_mpi\cluster_sparse_solver\source directory.

And here is the thread with a similar discussion which you can refer to for more details.

https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/MKL-PARDISO-get-L-U-factors/m-p/1265720#M31046


Please let us know if it helps in answering your queries.


Regards,

Vidya.


0 Kudos
VidyalathaB_Intel
Moderator
561 Views

Hi Huang,


Thanks for accepting our solution.

As the issue is resolved we are closing this thread. Please post a new question if you need any additional information from Intel as this thread will longer be monitored.


Have a Great Day!


Regards,

Vidya.


0 Kudos
Reply