15 September 2021   3:23:20.118 PM
 
SPARSE_INTERP_ND_TEST
  FORTRAN90 version.
  Test the SPARSE_INTERP_ND library.
  The R8LIB library is also required.

TEST01:
  Sparse interpolation for a function f(x) of M-dimensional argument.
  Use a sequence of sparse grids of levels 0 through SPARSE_MAX.
  Invoke a general Lagrange interpolant function to do this.

  Compare the exact function and the interpolants at a grid of points.

  The "order" is the sum of the orders of all the product grids
  used to make a particular sparse grid.

  Spatial dimension M =    1
  Maximum sparse grid level =    9
  Number of interpolation points is NI =    100

   L     Order    ApproxError

   1         1  0.24E-01
   2         3  0.49E-03
   3         5  0.18E-05
   4         9  0.21E-11
   5        17  0.17E-16
   6        33  0.35E-16
   7        65  0.49E-16
   8       129  0.56E-16
   9       257  0.82E-16
  10       513  0.12E-15

TEST01:
  Sparse interpolation for a function f(x) of M-dimensional argument.
  Use a sequence of sparse grids of levels 0 through SPARSE_MAX.
  Invoke a general Lagrange interpolant function to do this.

  Compare the exact function and the interpolants at a grid of points.

  The "order" is the sum of the orders of all the product grids
  used to make a particular sparse grid.

  Spatial dimension M =    2
  Maximum sparse grid level =    9
  Number of interpolation points is NI =    100

   L     Order    ApproxError

   1         1  0.21E-01
   2         7  0.61E-02
   3        25  0.52E-03
   4        67  0.59E-03
   5       161  0.22E-03
   6       371  0.81E-04
   7       837  0.22E-04
   8      1863  0.12E-04
   9      4105  0.16E-05
  10      8971  0.69E-06

TEST01:
  Sparse interpolation for a function f(x) of M-dimensional argument.
  Use a sequence of sparse grids of levels 0 through SPARSE_MAX.
  Invoke a general Lagrange interpolant function to do this.

  Compare the exact function and the interpolants at a grid of points.

  The "order" is the sum of the orders of all the product grids
  used to make a particular sparse grid.

  Spatial dimension M =    3
  Maximum sparse grid level =    9
  Number of interpolation points is NI =    100

   L     Order    ApproxError

   1         1  0.15E-01
   2        10  0.51E-02
   3        52  0.16E-02
   4       195  0.23E-03
   5       609  0.35E-03
   6      1710  0.27E-03
   7      4502  0.15E-03
   8     11369  0.50E-04
   9     27887  0.17E-04
  10     66936  0.81E-05

TEST01:
  Sparse interpolation for a function f(x) of M-dimensional argument.
  Use a sequence of sparse grids of levels 0 through SPARSE_MAX.
  Invoke a general Lagrange interpolant function to do this.

  Compare the exact function and the interpolants at a grid of points.

  The "order" is the sum of the orders of all the product grids
  used to make a particular sparse grid.

  Spatial dimension M =    4
  Maximum sparse grid level =    7
  Number of interpolation points is NI =    100

   L     Order    ApproxError

   1         1  0.12E-01
   2        13  0.52E-02
   3        87  0.16E-02
   4       411  0.59E-03
   5      1573  0.14E-03
   6      5257  0.16E-03
   7     16035  0.17E-03
   8     45879  0.13E-03
 
SPARSE_INTERP_ND_TEST
  Normal end of execution.
 
15 September 2021   3:23:21.147 PM