Summation-by-parts operators and high-order quadrature

Summation-by-parts (SBP) operators are finite-difference operators that mimic integration by parts. The SBP operator definition includes a weight matrix that is used formally for discrete integration; however, the accuracy of the weight matrix as a quadrature rule is not explicitly part of the SBP definition. We show that SBP weight matrices are related to trapezoid rules with end corrections whose accuracy matches the corresponding difference operator at internal nodes. For diagonal weight matrices, the accuracy of SBP quadrature extends to curvilinear domains provided the Jacobian is approximated with the same SBP operator used for the quadrature. This quadrature has significant implications for SBP-based discretizations; in particular, the diagonal norm accurately approximates the L2 norm for functions, and multi-dimensional SBP discretizations accurately approximate the divergence theorem.