On Fri, 9 Mar 2018, David Knezevic wrote:
- If I were to look into adding a special case for non-uniform 1st to 2nd order refinement for LAGRANGE variables, do you think this would be of interest to include that in libMesh, or would it be too specific to include? (I'd like to know if it's potentially of broader interest before looking further into this.)
Uninteresting to me, but I try like mad to avoid writing LAGRANGE-specific code, so if I want C0 p refinement I just use HIERARCHIC. Others who've coded themselves too far into a LAGRANGE-only corner might disagree... but from what I've seen the very first way to screw up is to assume that every node has a dof for every variable, and any code like that will *still* be broken if those users have second-order geometric elements (so they can support p=2) but start with p=1.
- How complex do you think it would be to add that special case?
Not very. I personally don't think it's worth it when you'd just end up stuck restricted by p<=2 anyways, but if anyone else disagrees I'd still happily merge their work. ;-) --- Roy ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users