Re: [Libmesh-users] p-refinement for Lagrange variables?

[Roy S. <roy...@ic...>]

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