We need to pass the state _with_ boundary data to the patch residual
evaluator, but compute a residual _without_ boundary data. We therefore
need all of the *withArtificial datastructures in PCPATCH

We need to pass the state _with_ boundary data to the patch residual
evaluator, but compute a residual _without_ boundary data. We therefore
need all of the *withArtificial datastructures in PCPATCH if we're
running in nonlinear mode. This requires passing two dofmaps to
the residual code, also.

show more ...

Make SNESPATCH a fully featured SNES (line searches, checking convergence criteria,
monitors, logging, etc)

Start refactoring for multiplicative updates. Multiplicative + SNESPATCH
doesn't actually work yet because we need to pass the state _with BC values_
to the residual evaluation, but at the moment we'

Start refactoring for multiplicative updates. Multiplicative + SNESPATCH
doesn't actually work yet because we need to pass the state _with BC values_
to the residual evaluation, but at the moment we're only passing the
state _for the dofs we're to solve for_.

show more ...

Remove the composition type on the SNESPATCH (it will all be
handled by the composition type on the underlying PC, although
it isn't yet)

Look for "pc_patch_construct_type" in PCPATCH, and
"snes_patch_construct_type" in SNESPATCH.

Small fixes

A PC only deals with two vectors, an RHS and an update.
A SNES deals with 4: residual calculation workspace, state vector, state update vector, optional RHS vector.

Rename things in PCPATCH to make

A PC only deals with two vectors, an RHS and an update.
A SNES deals with 4: residual calculation workspace, state vector, state update vector, optional RHS vector.

Rename things in PCPATCH to make clear what is RHS and what is update.
Scatter the global state onto the patch solve before calling SNESSolve.
Give each SNES its own separate residual calculation workspace.

show more ...

Some logic fixes. With a no-op residual evaluation the code now runs
through the solves on each patch, doing nothing each time (as the
residual is zero).

Set options prefix on inner PC before calling SetFromOptions

Pass the SNES' DM on to the patch PC

Various functions for setting data on a SNESPATCH.

Merge branch 'master' into knepley/feature-snes-patch

Not sure if this has done the right thing because there are no tests
yet. Also, the fortran stubs break (it doesn't know what to do with
SNESCo

Merge branch 'master' into knepley/feature-snes-patch

Not sure if this has done the right thing because there are no tests
yet. Also, the fortran stubs break (it doesn't know what to do with
SNESCompositeType), and I don't know how to fix it.

show more ...

SNES: Adde type Patch with trial implementation