[petsc-users] Solving with SNES

[behzad baghapour]

So, I am solving the Euler equations (inviscid compressible flow) using
discontinuous Galerkin method:

M (dQ/dt) = R,

R = nonlinear residual,
M = mass matrix,

Implicit procedure:
F = M*(Q^(n+1)-Q^n)/DT - R,
dF/dQ=J=M/DT-dR/dQ,
Newton => Q^(n+1) = Q^n - (dF/dQ)^(-1) F ( with proper preconditioning )

The nonlinear residual and Jacobian matrix are evaluated with a face-based
method. Two integrals are involved (the effect of element internal
connections and face flux connections). All of my calculations of flow
states are solved and stored in arrays of objects called elements.So my
residual and Derivative of residual routines are developed based on the
element (and face) objects. Then it is better for me to keep them when
using petsc as a nonlinear solver.

Hope these make any help.
Thanks,
Behzad
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111112/81738a30/attachment.htm>