Currently very verbose when a degenerate shape is encountered.
Adding mesh cell inertia calc.
Keeping both functions to ensure that the octant ordering specified in treeBoundBox is obeyed at the octree level, irrespective of the implementation of the basic boundBox.
Use new "T" function to avoid duplicate registration of grad(U) object.