Gaussian quadrature rules for C1 quintic splines with uniform knot vectors

We provide explicit quadrature rules for spaces of C1 quintic splines with uniform knot sequences over finite domains. The quadrature nodes and weights are derived via an explicit recursion that avoids numerical solvers. Each rule is optimal, that is, requires the minimal number of nodes, for a given function space. For each of n subintervals, generically, only two nodes are required which reduces the evaluation cost by 2/3 when compared to the classical Gaussian quadrature for polynomials over each knot span. Numerical experiments show fast convergence, as n grows, to the “two-third” quadrature rule of Hughes et al. (2010) for infinite domains.