I want to call pardiso with phase=33 to compute a solution to a given rhs. I'm sure I need to pass it the pt(:) array from the factoring phase, but do I also need to pass it any other arrays used in the factoring phase, like iparm, perm, a, ia and ja. In my calls to perform the factoring, I'm using iparm(1)=0 to let pardiso use all defaults.
I intend to pass the pt(:) array to different routines that need to compute solutions to different rhs's. I need to know if must also pass these other arrays.