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

FGMRES preconditioner applied to?

simon932
Beginner
520 Views

I am using the MKL's preconditioned FGMRES solver and I am trying to understand what exactly is the vector that FGMRES is asking to apply the preconditioner to. From the reference for the solver, Saad's Iterative methods for sparse linear systems, the left-preconditioned GMRES iteration (I'm assuming that FGMRES does left-preconditioning, please correct me if I'm wrong) involves computing at each step M^-1 A v_j.  That is, FGMRES first asks to compute the matrix vector product A v_j, and then I would assume that FGMRES would ask to apply the preconditioner on that result, i.e. compute M^-1 A v_j. Only, when I compute the squared 2 norm of the vectors involved, I get that the vector on which FGMRES asks to apply the preconditioner (which I would assume to be A v_j) always has a unit norm, regardless of the norm of A v_j. What is this unit norm vector that FGMRES is asking the user to apply the preconditioner to?

0 Kudos
3 Replies
simon932
Beginner
520 Views

Bump, I'm not an expert in iterative methods, just a user seeking more information, hopefully this question does not expose a gross misunderstanding of mine about the GMRES algorithm.

 

Thanks!

0 Kudos
Zhen_Z_Intel
Employee
520 Views

Hi,

MKL uses flexible GMRES which allows the preconditioner to vary from step to step. In the FGMRES, preconditioner will be changed at each step, and according to Saad's doc, the flexible seems do changes bases on right precondition that the action of AM^−1 j on a vector v do not in the span of Vm+1. MKL provides Incomplete LU factorization for preconditioners, you could use ILU0 & ILUT to apply it for FGMRES solver, please refer this document for detailed implement steps:

https://software.intel.com/en-us/node/521725
http://www.hpc.ut.ee/dokumendid/ips_xe_2015/composerxe/Documentation/en_US/mkl/mklman/GUID-A9A97100-09F8-46ED-A055-3DB4EF10DE0B.htm

Best regards,
Fiona

 

0 Kudos
simon932
Beginner
520 Views

Thank you Fiona, I had somehow missed the fact that F in FGMRES stands for Flexible, and it does say so in the documentation, my bad! Nevertheless I appreciate you clarifying this for me!

0 Kudos
Reply