is there a way to see the Content of Workspace for pardiso ( moslty call void PT ), 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 !
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 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.