> 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).
I don't know. I think what it takes the most is for someone to really think in detail about - what the ConstraintMatrix::condense() function does on vectors - how this interacts with automatic differentiation I think you could probably write all this down explicitly on a piece of paper for two cells, one of which is refined. I am not completely sure anything needs to be done at all, but I can't say I really thought about it in sufficient detail. Best W. -- ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
