On Wed, 27 Jan 2016, Roy Stogner wrote: > That's nothing I'd be able to get to myself in the immediate future, > although obviously "subtle inaccuracy when mixing time-dependent > Dirichlet + AMR" ought to be a high priority to fix.
On the other hand, Paul Bauman just IM'ed me to say that he'd be interested in fixing this ASAP. :-) So here's more of an outline in what needs to be done: 1. DofMap::process_constraints() needs to be fixed to properly update adjoint and not just primal rhs'es when expanding constraint equations. I'm doing this myself right now, partly because I know even more people it could affect, partly because it's very straightforward to fix, partly because it's a prerequisite for the rest, and partly because it was just plain embarrassing to run across. 2. AdjointDofConstraintValues needs to be changed to index over strings (vector names), not just integers. Loops over qoi_index which directly look up adjoint constraint rhs via that index need to instead prepend the rest of the adjoint vector name first. 3. We need a new enforce_constraints_exactly() version which takes a vector name. 4. We need to use this function (probably with a new overloaded argument to figure out which) in project_vector(). --- Roy ------------------------------------------------------------------------------ Site24x7 APM Insight: Get Deep Visibility into Application Performance APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month Monitor end-to-end web transactions and take corrective actions now Troubleshoot faster and improve end-user experience. Signup Now! http://pubads.g.doubleclick.net/gampad/clk?id=267308311&iu=/4140 _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
