combining direct pardiso and Fgmres?

Intel® oneAPI Math Kernel Library

Ask questions and share information with other developers who use Intel® Math Kernel Library.

combining direct pardiso and Fgmres?

200 Views

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?

Link Copied

0 Replies

Community support is provided during standard business hours (Monday to Friday 7AM - 5PM PST). Other contact methods are available here.

Intel does not verify all solutions, including but not limited to any file transfers that may appear in this community. Accordingly, Intel disclaims all express and implied warranties, including without limitation, the implied warranties of merchantability, fitness for a particular purpose, and non-infringement, as well as any warranty arising from course of performance, course of dealing, or usage in trade.

For more complete information about compiler optimizations, see our Optimization Notice.