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
