On Thu, 19 Nov 2009, Jed Brown wrote: > 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).
Good point. We're wary of AD (some nasty constitutive equations require calls to external codes which can't easily be retrofit with Sacado). But if memory does become a problem, our current "backup plan" is to only build the adjoint matrix one element submatrix at a time, use those for local submatrix-subvector multiplies, then assemble the resulting subvector products and throw away the submatrices. --- Roy ------------------------------------------------------------------------------ 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
