Possible bug in TreatmentOfEssentialEnv ?
TreatmentOfEssentialEnv
in xAlgorithm.h
contains the following statement:
xFilteredRegion<xIter, xAcceptUnion> dir_reg(all.begin(1), all.end(1), filter_dirichlet);
Meaning that the filter is applied only on edges, which is fine in 2D, but not in 3D. I would rather use:
xFilteredRegion<xIter, xAcceptUnion> dir_reg(all.begin(mesh->getDim()-1), all.end(mesh->getDim()-1), filter_dirichlet);
What do you think ?