Journal of Computational Physics 400 (2020) 108979
Contents lists available at ScienceDirect
Journal of Computational Physics
www.elsevier.com/locate/jcp
Generation of nested quadrature rules for generic weight functions via numerical optimization: Application to sparse grids ∗ Vahid Keshavarzzadeh a, , Robert M. Kirby a,b, Akil Narayan a,c a Scientific Computing and Imaging Institute, University of Utah, United States b School of Computing, University of Utah, United States c Department of Mathematics, University of Utah, United States a r t i c l e i n f o a b s t r a c t
Article history: We present a numerical framework for computing nested quadrature rules for various Received 3 August 2018 weight functions. The well-known Kronrod method extends the Gauss-Legendre quadrature Received in revised form 30 August 2019 by adding new optimal nodes to the existing Gauss nodes for integration of higher Accepted 17 September 2019 order polynomials. Our numerical method generalizes the Kronrod rule for any continuous Available online 27 September 2019 probability density function on real line with finite moments. We develop a bi-level Keywords: optimization scheme to solve moment-matching conditions for two levels of main and Nested quadrature nested rule and use a penalty method to enforce the constraints on the limits of the nodes Optimal quadrature and weights. We demonstrate our nested quadrature rule for probability measures on Sparse grids finite/infinite and symmetric/asymmetric supports. We generate Gauss-Kronrod-Patterson Numerical integration rules by slightly modifying our algorithm and present results associated with Chebyshev Polynomial approximation polynomials which are not reported elsewhere. We finally show the application of our nested rules in construction of sparse grids where we validate the accuracy and efficiency of such nested quadrature-based sparse grids on parameterized boundary and initial value problems in multiple dimensions. © 2019 Elsevier Inc. All rights reserved.
1. Introduction
A quadrature formula for integration takes the form