| 0904074f | 15-Dec-2018 |
Patrick Farrell <patrick@pefarrell.org> |
When solving a nonlinear problem, we need _all_ the information
about the state on the patch, including the values of dofs that
feature in global (as well as artificial) boundary conditions.
So we mu
When solving a nonlinear problem, we need _all_ the information
about the state on the patch, including the values of dofs that
feature in global (as well as artificial) boundary conditions.
So we must build this function space on each patch, too, and
the scatters from the smallest patch to this one.
show more ...
|
| 39fd2e8a | 14-Dec-2018 |
Patrick Farrell <patrick@pefarrell.org> |
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 ...
|
| 636c62a6 | 14-Dec-2018 |
Patrick Farrell <patrick@pefarrell.org> |
Make SNESPATCH a fully featured SNES (line searches, checking convergence criteria,
monitors, logging, etc) |
| 6c9c532d | 13-Dec-2018 |
Patrick Farrell <patrick@pefarrell.org> |
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 ...
|
| 1202d238 | 12-Dec-2018 |
Patrick Farrell <patrick@pefarrell.org> |
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 ...
|