SPARSE_GRID_CC_DATASET
Create Clenshaw-Curtis Sparse Grid Dataset


SPARSE_GRID_CC_DATASET is a FORTRAN77 program which creates a Clenshaw-Curtis sparse grid dataset.

Usage:

sparse_grid_cc_dataset dim_num level_max
where

Each sparse grid is stored using the "quadrature rule" format, that is, as three files:

The abscissas are ordered to respect the natural nesting of the sparse grids by level. That is, the file of level 3 points begins by listing the points in the level 2 grid.

Licensing:

The code described and made available on this web page is distributed under the GNU LGPL license.

Languages:

SPARSE_GRID_CC_DATASET is available in a C version and a C++ version and a FORTRAN77 version and a FORTRAN90 version and a MATLAB version.

Related Data and Programs:

MONTE_CARLO_RULE, a FORTRAN90 program which generates a dataset of N random M-dimensional points, regards it as a quadrature rule for the unit hypercube, and writes out three files of information.

SGMGA, a FORTRAN90 library which creates sparse grids based on a mixture of 1D quadrature rules, allowing anisotropic weights for each dimension.

SMOLPACK, a C library which implements Novak and Ritter's method for estimating the integral of a function over a multidimensional hypercube using sparse grids.

SPARSE_GRID_CC, a FORTRAN77 library which creates a Clenshaw-Curtis sparse grid.

SPARSE_GRID_CC, a dataset directory which contains examples of Clenshaw-Curtis sparse grids.

SPARSE_GRID_GL_DATASET, a FORTRAN90 program which creates a Gauss-Legendre sparse grid and write the data to three files.

SPARSE_GRID_HW, a FORTRAN77 library 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_DATASET, a FORTRAN90 program which creates a Gauss-Laguerre sparse grid and write the data to three files.

SPARSE_GRID_MIXED_DATASET, a FORTRAN90 program which creates a sparse grid dataset based on a mixture of 1D rules.

SPARSE_GRID_OPEN_DATASET, a FORTRAN90 program which creates a sparse grid dataset based on open rules (Fejer 2, Gauss-Patterson, Newton-Cotes-Open).

TOMS847, a MATLAB program which uses sparse grids to carry out multilinear hierarchical interpolation. It is commonly known as SPINTERP, and is by Andreas Klimke.

Reference:

  1. Volker Barthelmann, Erich Novak, Klaus Ritter,
    High Dimensional Polynomial Interpolation on Sparse Grids,
    Advances in Computational Mathematics,
    Volume 12, Number 4, 2000, pages 273-288.
  2. Thomas Gerstner, Michael Griebel,
    Numerical Integration Using Sparse Grids,
    Numerical Algorithms,
    Volume 18, Number 3-4, 1998, pages 209-232.
  3. Albert Nijenhuis, Herbert Wilf,
    Combinatorial Algorithms for Computers and Calculators,
    Second Edition,
    Academic Press, 1978,
    ISBN: 0-12-519260-6,
    LC: QA164.N54.
  4. Fabio Nobile, Raul Tempone, Clayton Webster,
    A Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data,
    SIAM Journal on Numerical Analysis,
    Volume 46, Number 5, 2008, pages 2309-2345.
  5. Sergey Smolyak,
    Quadrature and Interpolation Formulas for Tensor Products of Certain Classes of Functions,
    Doklady Akademii Nauk SSSR,
    Volume 4, 1963, pages 240-243.
  6. Dennis Stanton, Dennis White,
    Constructive Combinatorics,
    Springer, 1986,
    ISBN: 0387963472,
    LC: QA164.S79.

Source Code:

Examples and Tests:

Here are the three quadrature files created by the command

sparse_grid_cc_dataset 2 3

List of Routines:

You can go up one level to the FORTRAN77 source codes.


Last revised on 16 March 2013.