<div dir="ltr">Dear PETSc developers,<div><br></div><div>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: </div><div><br></div><div>1.  Add my solver name into MatSolverType.</div><div>2.  Register my solver by calling MatSolverTypeRegister</div><div><br></div><div>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</div><div><br></div><div>1. Am I doing it the right way? And if so</div><div>2. How to set such parameters as MatFactorType, lufactorsymbolic, lufactornumeric.</div><div><br></div><div>Many thanks,</div><div><br></div><div>Yuan</div><div>Ph.D. in Solid Mechanics</div></div>