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.
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.
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.
