07-Jan-2022 22:39:23

lagrange_interp_nd_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Test lagrange_interp_nd().
  The R8LIB library is needed.

LAGRANGE_INTERP_ND_TEST01:
  Interpolate in 1D, using orders.
  LAGRANGE_INTERP_ND_GRID sets the interpolant.
  LAGRANGE_INTERP_ND_VALUE evaluates it.

         Zinterp          Zexact      Error

        0.727527        0.727536    8.66e-06
        0.786924        0.786914    9.55e-06
        0.126641        0.126646    4.37e-06
        0.791582        0.791571    1.04e-05
        0.591026        0.591049    2.37e-05

LAGRANGE_INTERP_ND_TEST02:
  Interpolate in 2D, using orders.
  LAGRANGE_INTERP_ND_GRID sets the interpolant.
  LAGRANGE_INTERP_ND_VALUE evaluates it.

         Zinterp          Zexact      Error

        0.285976        0.290822    4.85e-03
        0.892463        0.892419    4.46e-05
        0.828141        0.829201    1.06e-03
        0.978622        0.978522    1.00e-04
        0.805267        0.805174    9.31e-05

LAGRANGE_INTERP_ND_TEST03:
  Interpolate in 3D, using orders.
  LAGRANGE_INTERP_ND_GRID sets the interpolant.
  LAGRANGE_INTERP_ND_VALUE evaluates it.

         Zinterp          Zexact      Error

        0.851145        0.851127    1.83e-05
         0.98641        0.986532    1.22e-04
        0.952994         0.95294    5.40e-05
        0.951515        0.952029    5.14e-04
        0.706125        0.705277    8.48e-04

LAGRANGE_INTERP_ND_TEST04:
  Interpolate in 3D, using orders.
  Use a sequence of increasing orders.

  Level     Order   Average Error

      0         1    3.32e-02
      1        27    4.95e-03
      2       125    1.91e-04
      3       729    4.72e-06
      4      4913    2.70e-08
      5     35937    6.09e-10

LAGRANGE_INTERP_ND_TEST05:
  Repeat test #1, using levels.
  LAGRANGE_INTERP_ND_GRID2 sets the interpolant.
  LAGRANGE_INTERP_ND_VALUE2 evaluates it.

         Zinterp          Zexact      Error

       0.0757693       0.0757816    1.22e-05
       0.0539113        0.053924    1.26e-05
        0.506214        0.506221    6.88e-06
        0.702671        0.702687    1.59e-05
        0.804023        0.804011    1.17e-05

LAGRANGE_INTERP_ND_TEST06:
  Repeat test #2, using levels.
  LAGRANGE_INTERP_ND_GRID2 sets the interpolant.
  LAGRANGE_INTERP_ND_VALUE2 evaluates it.

         Zinterp          Zexact      Error

         0.55308        0.550922    2.16e-03
        0.452125        0.452477    3.52e-04
        0.359459         0.36544    5.98e-03
        0.754842        0.753309    1.53e-03
         0.52719        0.525994    1.20e-03

LAGRANGE_INTERP_ND_TEST07:
  Repeat test #3,  using levels.
  LAGRANGE_INTERP_ND_GRID2 sets the interpolant.
  LAGRANGE_INTERP_ND_VALUE2 evaluates it.

         Zinterp          Zexact      Error

        0.799501        0.799831    3.31e-04
         0.89643        0.896823    3.93e-04
        0.809899        0.810688    7.89e-04
        0.839979        0.840117    1.38e-04
        0.887949         0.88795    4.23e-07

LAGRANGE_INTERP_ND_TEST08:
  Interpolate in 3D, using levels.
  Use a sequence of increasing levels.

  Level     Order   Average Error

      0         1    4.76e-02
      1        27    8.30e-03
      2       125    3.88e-04
      3       729    1.16e-05
      4      4913    4.75e-08
      5     35937    2.63e-10

LAGRANGE_INTERP_ND_TEST09:
  Interpolate in 3D, using orders.
  Use a sequence of increasing orders.
  Use anisotropic resolution.
  The interpoland is a polynomial of degrees 3, 5, 2
  so our orders need to be at least 4, 6, 3 to match it.

  Level     Orders   Average Error

      0   1   1   1    6.10e-02
      1   2   1   1    5.88e-02
      2   1   2   1    6.21e-02
      3   1   1   2    6.18e-02
      4   4   2   2    2.52e-02
      5   2   4   2    6.71e-02
      6   2   2   4    8.86e-02
      7   2   2   4    8.86e-02
      8   6   4   4    3.41e-04
      9   4   6   4    8.38e-17
     10   4   4   6    3.41e-04

lagrange_interp_nd_test():
  Normal end of execution.

07-Jan-2022 22:39:25