Add find_bounding_box
Distinguish between interpolatory quadrature and quadrature element groups