Remove out-commented code
Add composite quadrature rules for simplex elements

This short file implements a general mechanism for constructing
composite quadrature rules.  It takes some quadrature rule
and an integer refinement level as the input.  The reference
simplex is the refined uniformly as prescribed by the desired
level.  The given input quadrature rule is then used on
each subelement.

Currently this only works for simplices.  I'll try to implement
it also for other elements, but there are some fixes in the
virtual refinement code needed first.
Use 'coorddim' instead of 'dim' for coordinates, even when the two are the same
Bugfix: jacobian(Inverse)Transpose may crash for 0d geometries
Test 0d-in-3d geometry constructed with the 3d-in-3d constructor.

This is not as stupid as it sounds: the given lower left corner
should then be taken as coordinate vector.
Implement geometries for static-refinement vertices

The refinement mechanism allows to iterative over the vertices
of a refinement.  However, unlike the elements of a refinement,
the vertices didn't have a 'geometry' method (only a method
coord()).  This patch adds the missing geometry method for
the static refinement.
Test vertex geometries of static refinement
Test the composite quadrature rule for simplices