topic Getting the content of *pt in Pardiso in Intel® oneAPI Math Kernel Library

Getting the content of *pt in Pardiso

Orlando_R_ — Thu, 24 Sep 2015 07:49:01 GMT

is there a way to see the Content of Workspace for pardiso ( moslty call void PT[64] ), It seems that the Permutation Matrix of PAx=Py preserves quite well the sparsity in LU, or LL* , in my app I need only L (for Hermitian mtx) . I would like to get L and P from pt

Any pointer will be gretly appreciated !

Thanks

Hi Orlando,

Alexander_K_Intel2 — Thu, 24 Sep 2015 07:54:53 GMT

Hi Orlando,

Currently MKL pardiso have not support external access to elements of L and U. But there different way how to handle it - can i ask the reason of request to access of L and U?

Thanks,

Alex

Hi Alex,

Orlando_R_ — Wed, 30 Sep 2015 08:28:01 GMT

Hi Alex,

thanks for your answer. In my application I am not solving a System directly Ax=b , What I have is res=L*c where c is the sum of other vectors, for now I care only about res. My Memory usage is quite high having L and U and P in pt. I can use phase 1 and 2, store L and the Permutaion matirx P , release pardiso , discarding U.

Regards

Orlando