Sharp Stone <[email protected]> writes: > Hi all, > > I'm sorry I'm still confusing that if I have calculated the Jacobian > elements, how can I go ahead to solve the linear system with dof>1 at each > node? The example code ex19 solves the linear system with Guass-Seidel > method,
As we just said, the same code (ex19.c) can use many different methods. The default method is Newton-Krylov with block Jacobi preconditioning. Run with -snes_view to see details.
pgpxvqMFmYmzM.pgp
Description: PGP signature
