combining direct pardiso and Fgmres?

combining direct pardiso and Fgmres?

i am currently writting a nonlinear finite element code. so i have to solve a nonlinear system.
for each iteration my matrixes changes, i use direct solution for each matrix sparsity maybe change in each iteration so i cant use pardiso iterative solver at phase 13 and phase 23 as suggested here:

in this thread:

Alexander Kalinkin sugested we can use pardiso  LU decomposition of first matrix as preconditioner for Fgemres iterative solver.

i want to know how i can do that?

1 post / 0 new
For more complete information about compiler optimizations, see our Optimization Notice.