**spquad**,
a MATLAB code which
can compute a sparse grid quadrature rule to estimate a multidimensional integral.
The quadrature rules created by the program are based on the Clenshaw-Curtis rule.
The code is by Greg von Winckel.

The original version of this software is available at the MATLAB Central File Exchange: https://www.mathworks.com/MATLABcentral/fileexchange/19063-sparse-grid-quadrature/content/spquad.m.

Note that the **spquad()** function, by default, assumes that the integration region
is [-1,+1]^d, where **d** is the spatial dimension. The function includes an optional
argument **bpt**, which is an array allowing the user to specify the lower and upper
endpoints of integration. However, in the experience of the installer, this extra argument
does not seem to have the correct effect, resulting in quadrature rules that are not
correctly shifted to the new domain. A note has been sent to the author asking for
clarification. In the meantime, the user can compute the rule for [-1,+1]^d and use a
simple linear transformation to adjust the points and weights to a nondefault integration region.

No licensing information was given with the text of the file at the MATLAB Central File Exchange.

**spquad** is available in
a MATLAB version

sparse_grid_cc, a MATLAB code which can define a multidimensional sparse grid based on a 1d clenshaw curtis rule.

sparse_grid_gl, a MATLAB code which creates sparse grids based on gauss-legendre rules.

sparse_grid_hermite, a MATLAB code which creates sparse grids based on gauss-hermite rules.

sparse_grid_hw, a MATLAB code which creates sparse grids based on gauss-legendre, gauss-hermite, gauss-patterson, or a nested variation of gauss-hermite rules, by florian heiss and viktor winschel.

sparse_grid_laguerre, a MATLAB code which creates sparse grids based on gauss-laguerre rules.

spinterp, a MATLAB code which carries out piecewise multilinear hierarchical sparse grid interpolation; an earlier version of this software is ACM TOMS Algorithm 847, by Andreas Klimke;

Greg von Winckel

- spquad.m returns the points and weights of a sparse grid for dimension D and level K, based on the Clenshaw Curtis rule.