Add a factory to construct Gaussian quadrature objects
Add reduced integration scheme