I try to use PARDISO as an iterative solver not a direct solver. I have problem for defining IPARMs and DPARMs. I want an example to show hove can i use them and which phases i need to call. I have read the manual and also the example for pardiso_symm_iter.f but in that one, the auther has used default iparm and dparm. In iterative solver do i need to call other phases that phase=3? I want iterative method with conjugate gradient preconditioner . I appreciate if some one can help me.
iparm(4) = 0 phase = 13 PARDISO(phase, A_0,…) Do I = 1, number_of_matrix iparm(4) = 61 phase=23 PARDISO(phase, A_i,…) EnddoPlease let us know if this solves your problem. Thanks, Zhang