[petsc-users] How to introduce external iterative solver into PETSc

[袁煕]

Dear PETSc developers,

I have my own linear solver and am trying to put it into PETSc as an
external solver.  Following the implementation of mumps, mkl_cpardiso,
supelu etc, I think I should do the follow:

1.  Add my solver name into MatSolverType.
2.  Register my solver by calling MatSolverTypeRegister

to let petsc record the existence of a new solver. The problem is that the
above external solvers are all direct solvers, and a MatFactorType
parameter should be set to indicate its factorization type, such as LU, QR,
or Cholesky. But my solver is an iterative one, that means I cannot specify
its MatFactorType. I wish to understand

1. Am I doing it the right way? And if so
2. How to set such parameters as MatFactorType,
lufactorsymbolic, lufactornumeric.

Many thanks,

Yuan
Ph.D. in Solid Mechanics
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221205/2ee8f23d/attachment.html>