topic Hi Gennady, in Intel® oneAPI Math Kernel Library

Can I free original matrix after Pardiso factorization

lengfei_h_ — Fri, 13 Jan 2017 23:54:56 GMT

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

yes, that's correct.

Gennady_F_Intel — Sun, 15 Jan 2017 05:11:09 GMT

yes, that's correct.

Hi Gennady,

lengfei_h_ — Sun, 15 Jan 2017 05:20:11 GMT

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

some more clarification: In

Gennady_F_Intel — Tue, 17 Jan 2017 06:42:07 GMT

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.