Added composite hexahedral rules based on the Gauss-Legendre quadrature rule
Added a generic stiffness matrix routine