09 May 2025   9:16:52.793 PM

lagrange_interp_nd_test():
  FORTRAN90 version
  Test lagrange_interp_nd().
  The R8LIB library is needed.

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

         Zinterp          Zexact      Error

    0.442095        0.442086        0.94E-05
    0.710794        0.710808        0.14E-04
    0.209516        0.209501        0.15E-04
    0.791311        0.791301        0.10E-04
    0.835530        0.835526        0.42E-05

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

         Zinterp          Zexact      Error

    0.840769        0.840926        0.16E-03
    0.763250        0.762932        0.32E-03
    0.580150        0.577254        0.29E-02
    0.961127        0.961119        0.84E-05
    0.702283        0.699640        0.26E-02

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

         Zinterp          Zexact      Error

    0.939020        0.939008        0.12E-04
    0.983632        0.983410        0.22E-03
    0.538514        0.538273        0.24E-03
    0.333019        0.338486        0.55E-02
    0.598913        0.598045        0.87E-03

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

  Level     Order   Average Error

      0         1    0.46E-01
      1        27    0.78E-02
      2       125    0.36E-03
      3       729    0.11E-04
      4      4913    0.17E-06
      5     35937    0.17E-08

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

         Zinterp          Zexact      Error

    0.286094        0.286069        0.26E-04
    0.479454        0.479454        0.72E-08
    0.195155        0.195143        0.12E-04
    0.667483        0.667506        0.23E-04
    0.430694        0.430682        0.12E-04

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

         Zinterp          Zexact      Error

    0.190367        0.188491        0.19E-02
    0.564128        0.564789        0.66E-03
    0.854787        0.854902        0.11E-03
    0.894211        0.894698        0.49E-03
    0.686730        0.683307        0.34E-02

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

         Zinterp          Zexact      Error

    0.932648        0.932771        0.12E-03
    0.697544        0.696996        0.55E-03
    0.661050        0.658493        0.26E-02
    0.906009        0.906491        0.48E-03
    0.876303        0.876352        0.49E-04

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

  Level     Order   Average Error

      0         1    0.50E-01
      1        27    0.77E-02
      2       125    0.41E-03
      3       729    0.21E-04
      4      4913    0.11E-05
      5     35937    0.42E-07

TEST09:
  Interpolate in 3D, using orders.
  Use a sequence of increasing orders.
  Use anisotropic resolution.
  The interpolant 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    0.11E+00
      1      2      1      1    0.65E-01
      2      1      2      1    0.99E-01
      3      1      1      2    0.10E+00
      4      4      2      2    0.49E-01
      5      2      4      2    0.66E-01
      6      2      2      4    0.91E-01
      7      2      2      4    0.91E-01
      8      6      4      4    0.15E-02
      9      4      6      4    0.14E-15
     10      4      4      6    0.15E-02

lagrange_interp_nd_test():
  Normal end of execution.

09 May 2025   9:16:52.812 PM