For tets, wedges, and pyramids, using trilinear shape functions for value interpolation may not be linearly exact because the trilinear integration point location may not lie at the geometric centroid of the face on which the integration point is located. In such situations, interpolation accuracy may be improved by using geometric shape functions. The default behavior is to use geometric shape functions for buoyant problems and the standard parametric shape functions otherwise.