This PR fixes a long standing issue in the boundary area calculation. Before using the FaceElementTransformations object, we have to set the integration point by calling SetIntPoint, which we had failed to do. Since we're only using the FaceElementTransformations::Weight function, which won't actually depend on the quadrature point for linear elements, it seems we were getting away with it... until #168.
After b5ea27a, I have been unable to reproduce #168, so this PR resolves #168.
This PR fixes a long standing issue in the boundary area calculation. Before using the
FaceElementTransformationsobject, we have to set the integration point by callingSetIntPoint, which we had failed to do. Since we're only using theFaceElementTransformations::Weightfunction, which won't actually depend on the quadrature point for linear elements, it seems we were getting away with it... until #168.After b5ea27a, I have been unable to reproduce #168, so this PR resolves #168.