On Wed, 3 Aug 2011, Truman Ellis wrote: > I would like to have discontinuous elements other than monomial. In > order to add a discontinuous hierarchic or lagrange element would I just > modify their n_dofs_at_node to always return 0 and n_dofs_per_elem to > return the total number of dofs?
Copying this to libmesh-devel in case anyone sees something I'm missing, but yes, that ought to be 90% of the work. The other 10% is that adding any new FE types requires manual entries in all libMesh's many switch-statements-masquerading-as-vtables. Basically grep for HIERARCHIC to find all those spots and add a new DISC_HIERARCHIC and DISC_LAGRANGE entry as appropriate. > I am hoping it is that easy... Me too. It would be very nice to have a discontinuous basis that wasn't so poorly conditioned for high p. It would also be nice practice for you for adding new FE types; for DPG problems you'll need another new type to handle flux degrees of freedom (in which type dofs will only be defined on face nodes), right? --- Roy ------------------------------------------------------------------------------ BlackBerry® DevCon Americas, Oct. 18-20, San Francisco, CA The must-attend event for mobile developers. Connect with experts. Get tools for creating Super Apps. See the latest technologies. Sessions, hands-on labs, demos & much more. Register early & save! http://p.sf.net/sfu/rim-blackberry-1 _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
