I used PARDISO to solve large sparse matrix with phase = 33 many times, I found that the memory usage increased everytime the PARDISO is called with phase=33. So I checked the sample program 'pardiso_sym_c.c' with Intel Inspector. The conculsion is as follows:
1) When call the PARDISO 100 times with phase = 33, the object size by memory leak is 28800;
2) When call the PARDISO 1000 times with phase = 33, the object size by memory leak is 288000.
I also attached the whole solution for you to reproduce the problem.
Could you please help me to take a look at it?