Hi all I'm interested in using the Sacado package in Trilinos to numerically evaluate the global Newton matrix for a problem in which I have have hanging nodes handled via standard hanging-node constraints (i.e strongly). What I would like to know is how hard is this to do?
I have read the warning in tut-33... "because we would have to thread the automatic differentiation we use to compute the elements of the Newton matrix from the residual through the operations of the ConstraintMatrix class. This would be possible, but is not trivial, and so we choose this alternative approach." Many thanks Andrew
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
