What is the difference between the fe and feAux fields in PetscFEM? Do they (and feBd) have to refer to distinct PetscFE objects, or can I reuse one set of finite element spaces?
In particular, based on snes ex12, it doesn't look like feBd contains the logic for what the boundary is. This information appears to be extracted by DMPlexComputeResidualFEM and friends. If so, does that mean PetscFEM doesn't have support for several different boundary routines evaluated over different parts of the boundary? I don't necessarily need this support at the moment, but I want to check whether I'm reading the code correctly. Thanks, Geoffrey
