Use constexpr integration points number in for loop.

This allows loop unrolling.
Add ShapeMatrixType defining types.

The matrix types are dependent on given element's dimension and number of nodal
functions.

The particular version is for fixed size eigen matrices called
EigenFixedSizeShapeMatrices.
In X struct use shape matrix types from template.