On Fri, 2 Oct 2009, John Peterson wrote: > In FEMSystem systems, it might be cool if there was an additional > virtual function, something like > > virtual std::vector<std::vector<Real> >& user_test_function(DiffContext&); > > which could be overloaded in derived classes to compute any "custom" > test function the user desired, given current values in a DiffContext. > This, combined with a user-provided forcing, could conceivably then > be used to add an "automatic" forcing. Also, this code could be > re-used in normal assembly as well so there would be no duplication.
This is exactly what I was thinking about, but I think it needs one more user-provided function: if the user also supplies derivatives of test functions with respect to components of the (fixed) solution, then stabilized formulations might get a bit simpler: in addition to reducing some of the complexity in time_derivative and constraint methods, you ought to be able to forgo overriding mass_residual entirely. --- Roy ------------------------------------------------------------------------------ Come build with us! The BlackBerry® Developer Conference in SF, CA is the only developer event you need to attend this year. Jumpstart your developing skills, take BlackBerry mobile applications to market and stay ahead of the curve. Join us from November 9-12, 2009. Register now! http://p.sf.net/sfu/devconf _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
