Pardiso Access Violation in multiple calls

[WaterDog]

I'm using MKL versiion 10.2. I'm trying a program that calls Pardiso multiple times. Each time has multiple right hand sides. PARDISO works good for the first call and it gaves access violation error for second call. As the error happens, the program will freeze the computer and I have to push the power button to turn off the computer. I set PT()=0 after the first call. I follow the procedure in pardiso_sym_f.f to call Pardiso. Any help?


WaterDog

[WaterDog]

Update: the problem is solved. PT()=0 is actually needed after releasing memory (phase =-1).

[Alexander_K_Intel2]

Hi,

Sorry for delay. In MKL manual pt describe in next way: "ptis an integer array with 64 entries. It is very important that the pointer is initialized with zero at the first call of...". So you need to initialize pt by zero before call to pardiso, so if you call pardiso several times you need to set it by zero before each initialization step.

With best regards,

Alexander Kalinkin