Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.
Announcements
This community is designed for sharing of public information. Please do not share Intel or third-party confidential information here.
6590 Discussions

Can I free original matrix after Pardiso factorization

lengfei_h_
Beginner
120 Views

Hi, 

I am wondering when can I release the original matrix memory (*a, *ia, *ja)? Since pardiso will keep the factor matrix internally, I think the original matrix memory can be release after factorization.

Any comments would be great.

Thanks

Lengfei

0 Kudos
3 Replies
Gennady_F_Intel
Moderator
120 Views

yes, that's correct.

lengfei_h_
Beginner
120 Views

Hi Gennady,

Thanks for the reply. I tried several matrices. Some of them work fine when I free the original matrix before going to solve phase. However, one case gave me zero Pivot error. If I free the original matrix after all solve steps, it will run without giving any error.

Thanks

Lengfei

Gennady_F_Intel
Moderator
120 Views

some more clarification: In the case if iparm[7] == 0 and iparm[20] == 2 || 3 ( in C notation),  You may delete the original matrix and call solving stage ( 33 ).This is applied to the all types of matrixes. In the case if you will see some sort of failure or hanging  at the Solving stage, then it may be the real issue. pls report this.

Reply