Hello,
I would like to use ISS from MKL for solving linear systems of equations, originating from discretization of partial differential equations used in CFD.
I found out that CG and FGMRES are available, but not everything is clear from the documentation:
- Can ISS be threaded? I am using Pardiso and it can be made parallel easily.
- Is there any example besides those in \example\solverf\source that would demonstrate how the matrix of the linear system is constructed and used by ISS? I might have overlooked information about this and the matrix format (CSR?) in the documentation.
- Is there any example that would show how preconditioning should be done? There is a general description of this in the documentation, but it does not feel sufficient to make sure the implementation is correct, e.g. in case of "apply the preconditioner C inverse to tmp".
By the way, can Pardiso as a direct solver work also in an iterative "mode", meaning that it would not calculate the solution, but just a user-defined number of iterations?
Thank you.