[petsc-users] Augmented Lagrangian examples?

Matthew Knepley knepley at gmail.com
Mon Aug 10 10:01:07 CDT 2015

On Mon, Aug 10, 2015 at 9:56 AM, Justin Chang <jychang48 at gmail.com> wrote:

> Hi all,
> 1) I ran across this paper:
> http://web.stanford.edu/~egawlik/pdf/GaMuSaWi2012.pdf
> and was wondering if there are any current TAO examples that do this.
> also
> 2) If I integrate this into an FEM (from DMPlex) I will need to assemble
> an equality jacobian matrix and constraint vector. But the element-wise
> constraints that I need to compute (e.g., the divergence) needs all degrees
> of freedom within the element closure including the essential BCs DMPlex
> removes from the global matrix/vector. So how can I work around this and/or
> access said removed terms inside a TAO routine?

Those are not really dofs, they are boundary values, so they are not in the
global system. The local vectors have
the boundary values, so you can calculate the correct constraint eqns to
put in the global system.


> Thanks,
> Justin

What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150810/19286a5a/attachment.html>

More information about the petsc-users mailing list