combining direct pardiso and Fgmres?

combining direct pardiso and Fgmres?

aslan666's picture

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 iteration.my matrix sparsity maybe change in each iteration so i cant use pardiso iterative solver at phase 13 and phase 23 as suggested here:

http://software.intel.com/en-us/forums/topic/326721

in this thread:

http://software.intel.com/en-us/forums/topic/389822

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.