Boundary elements for HDivConstant spaces

[jefersondossa]

There is an issue when handling with HDivConstant spaces: The boundary elements need to be oriented in PZ standard. In my external project, I used to delete all boundary elements coming from the geometric mesh and then recreate them with TPZGeoElBC to avoid this problem. As it is, TPZHDivApproxCreator is not working for HdivConstant spaces (mostly in 3D), unless the mesh boundary is properly oriented. Any suggestions to improve robustness or we should just incorporate the code to delete and recreate the boundary elements in TPZHDivApproxCreator when we have HDivConstant spaces?

[philippedevloo]

I believe the verification of the orientation is necessary. In the "standard" H(div) spaces, the value of the shape functions can be altered by the variable fSideOrient. This is not possible (I believe) for spaces derived from H(curl)