This version increases integration performance on meshes by neglecting to loop over elements. However, we still have a (small) loop over the elements and integration points on the reference triangle.
1.0.0
Changes - updated readme - changed signature of main integration routines to be more user-friendly, i.e. you only need to pass the `CubatureRuleEnum` to both `integrate_on_mesh` and `integrate_on_triangle`