Roy Stogner wrote:

> Block variable ordering with "--node_major_dofs" is perhaps the least
> well-documented option in libMesh, which is saying a lot, but we're
> already using it.

Neat.

> We'll have to construct the true Jacobian for the adjoint problem, but
> might try to set up a preconditioner based on reduced physics or a
> coarser grid.

You still don't need to assemble the adjoint, you just need to be able
to multiply with the adjoint.  Finite differencing can't do this for
you, but AD can and it's easy to do yourself when you have a Galerkin
method (basically just swap test and trial functions in the right
place).

Jed

Attachment: signature.asc
Description: OpenPGP digital signature

------------------------------------------------------------------------------
Let Crystal Reports handle the reporting - Free Crystal Reports 2008 30-Day 
trial. Simplify your report design, integration and deployment - and focus on 
what you do best, core application coding. Discover what's new with
Crystal Reports now.  http://p.sf.net/sfu/bobj-july
_______________________________________________
Libmesh-devel mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-devel

Reply via email to