6 November  2012   9:51:38.919 PM      
 
SPARSE_INTERP_ND_PRB
  FORTRAN77 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.27E-01
   2         3  0.51E-03
   3         5  0.16E-05
   4         9  0.22E-11
   5        17  0.14E-16
   6        33  0.25E-16
   7        65  0.32E-16
   8       129  0.46E-16
   9       257  0.80E-16
  10       513  0.14E-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.74E-02
   3        25  0.59E-03
   4        67  0.59E-03
   5       161  0.22E-03
   6       371  0.94E-04
   7       837  0.31E-04
   8      1863  0.13E-04
   9      4105  0.45E-05
  10      8971  0.20E-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 =    3
  Maximum sparse grid level =    9
  Number of interpolation points is NI =    100

   L     Order    ApproxError

   1         1  0.15E-01
   2        10  0.56E-02
   3        52  0.22E-02
   4       195  0.25E-03
   5       609  0.47E-03
   6      1710  0.20E-03
   7      4502  0.98E-04
   8     11369  0.30E-04
   9     27887  0.15E-04
  10     66936  0.74E-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.14E-01
   2        13  0.61E-02
   3        87  0.24E-02
   4       411  0.12E-02
   5      1573  0.13E-03
   6      5257  0.41E-03
   7     16035  0.21E-03
   8     45879  0.16E-03
 
SPARSE_INTERP_ND_PRB
  Normal end of execution.
 
 6 November  2012   9:51:44.248 PM