Getting the content of *pt in Pardiso

[Orlando_R_]

Hi

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

[Alexander_K_Intel2]

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

[Orlando_R_]

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