"De Groof, Vincent Frans Maria" <[email protected]> writes:
> I want to do some experiments with my own iterative solvers. Is there > any alternative for a user-defined KSP similar to PCSHELL? There is not a KSPSHELL, but you can create and register new implementations from user code. I recommend copying an existing implementation most like what you want to implement, or start from something simple like src/ksp/ksp/impls/cr/cr.c. Namespace the implementation, compile it, and call KSPRegister() somewhere in your application. If it turns out to work better than other options for some problem of interest, please consider contributing to PETSc. > And a question on the side. Are the orthogonalization methods ( (M)GS > ) readily available? Or are they private for GMRES? The GMRES implementation is incremental. If you have all the vectors at once, it's just a couple lines of code using VecMDot and VecMAXPY. We may add a general-purpose "TSQR" in the near future.
pgpsErcWTqsFX.pgp
Description: PGP signature
