Community
cancel
Showing results for 
Search instead for 
Did you mean: 
lengfei_h_
Beginner
49 Views

Can I free original matrix after Pardiso factorization

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
49 Views

yes, that's correct.

lengfei_h_
Beginner
49 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
49 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