I have to solve a sequence (iteration) of systems AX = B with multiple RHS.
At each iteration, X and B do not change and A changes but retains its sparsity patterns.
So I think I can use iparm=62 but I know exactly how to use it.
Below is a summary of what I do :
Itration 1 :
Itration 2 :
Itration 3 :
I never release memory (phase=0 or -1) between iterations.
Is it the correct way to solve my problem ?