Remove unused import.
Remove unnecessary mut.
Use optimal order in Gaussian quadrature.
Make aux integration functions standalone.  Unify tags.
Implement extra subdivisions in tagged regions of input mesh.