Generation of nested quadrature rules for generic weight functions via numerical optimization: application to sparse grids

DOI10.1016/J.JCP.2019.108979zbMATH Open1453.65053arXiv1808.03707OpenAlexW2975137534WikidataQ127191188 ScholiaQ127191188MaRDI QIDQ2222671FDOQ2222671

Authors: V. Keshavarzzadeh, Robert M. Kirby, Akil C. Narayan 

Publication date: 27 January 2021

Published in: Journal of Computational Physics (Search for Journal in Brave)

Abstract: We present a numerical framework for computing nested quadrature rules for various weight functions. The well-known Kronrod method extends the Gauss-Legendre quadrature by adding new optimal nodes to the existing Gauss nodes for integration of higher order polynomials. Our numerical method generalizes the Kronrod rule for any continuous probability density function on real line with finite moments. We develop a bi-level optimization scheme to solve moment-matching conditions for two levels of main and nested rule and use a penalty method to enforce the constraints on the limits of the nodes and weights. We demonstrate our nested quadrature rule for probability measures on finite/infinite and symmetric/asymmetric supports. We generate Gauss-Kronrod-Patterson rules by slightly modifying our algorithm and present results associated with Chebyshev 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.

Full work available at URL: https://arxiv.org/abs/1808.03707

Recommendations

• Algorithms for the construction of high-order Kronrod rule extensions with application to sparse-grid integration
• Generating nested quadrature rules with positive weights based on arbitrary sample sets
• Stratified nested and related quadrature rules
• Numerical integration in multiple dimensions with designed quadrature
• Numerical integration using sparse grids

zbMATH Keywords

numerical integrationpolynomial approximationsparse gridsoptimal quadraturenested quadrature

Mathematics Subject Classification ID

Numerical optimization and variational techniques (65K10) Numerical quadrature and cubature formulas (65D32) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60)

Cites Work

• Fully symmetric interpolatory rules for multiple integrals over infinite regions with Gaussian weight
• Sparse Optimization with Least-Squares Constraints
• Probing the Pareto frontier for basis pursuit solutions
• Title not available (Why is that?)
• Title not available (Why is that?)
• The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems
• Sparse grids
• Likelihood approximation by numerical integration on sparse grids
• Numerical integration using sparse grids
• Title not available (Why is that?)
• Title not available (Why is that?)
• Calculation of Gauss Quadrature Rules
• High-Order Collocation Methods for Differential Equations with Random Inputs
• Adaptive Leja sparse grid constructions for stochastic collocation and high-dimensional approximation
• Title not available (Why is that?)
• The Optimum Addition of Points to Quadrature Formulae
• Explicit cost bounds of algorithms for multivariate tensor product problems
• Adaptive sparse grid algorithms with applications to electromagnetic scattering under uncertainty
• Construction of Gauss-Christoffel Quadrature Formulas
• The circle theorem and related theorems for Gauss-type quadrature rules
• Numerical integration in multiple dimensions with designed quadrature

Cited In (7)

• A geometrical interpretation of the addition of nodes to an interpolatory quadrature rule while preserving positive weights
• A new multivariate quadrature rule for calculating statistical moments of stochastic response
• Algorithms for the construction of high-order Kronrod rule extensions with application to sparse-grid integration
• Efficient quadratures for high-dimensional Bayesian data assimilation
• Non-intrusive uncertainty quantification using reduced cubature rules
• Numerical integration in multiple dimensions with designed quadrature
• Generating nested quadrature rules with positive weights based on arbitrary sample sets

Uses Software

• SPGL1
• nwSpGr
• PATSYM
• HRMSYM

This page was built for publication: Generation of nested quadrature rules for generic weight functions via numerical optimization: application to sparse grids

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2222671)