15 September 2021   8:27:05.584 AM

LAGRANGE_INTERP_ND_TEST:
  FORTRAN90 version
  Test the LAGRANGE_INTERP_ND library.
  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.147727        0.147726        0.43E-06
    0.292315        0.292289        0.26E-04
    0.680839        0.680860        0.21E-04
    0.695232        0.695250        0.18E-04
    0.463772        0.463768        0.40E-05

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

         Zinterp          Zexact      Error

    0.771167        0.770812        0.35E-03
    0.796560        0.797536        0.98E-03
    0.864320        0.864426        0.11E-03
    0.746707        0.745196        0.15E-02
    0.736863        0.735224        0.16E-02

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

         Zinterp          Zexact      Error

    0.849191        0.849590        0.40E-03
    0.827371        0.828498        0.11E-02
    0.971872        0.971844        0.28E-04
    0.954639        0.955111        0.47E-03
    0.895286        0.895968        0.68E-03

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

  Level     Order   Average Error

      0         1    0.36E-01
      1        27    0.51E-02
      2       125    0.20E-03
      3       729    0.54E-05
      4      4913    0.22E-07
      5     35937    0.43E-09

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

         Zinterp          Zexact      Error

    0.806153        0.806142        0.12E-04
    0.353698        0.353673        0.25E-04
    0.602272        0.602296        0.25E-04
    0.616962        0.616988        0.26E-04
    0.610117        0.610143        0.25E-04

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

         Zinterp          Zexact      Error

    0.453657        0.454673        0.10E-02
    0.422585        0.426353        0.38E-02
    0.878607        0.879120        0.51E-03
    0.690033        0.688910        0.11E-02
    0.174610        0.170691        0.39E-02

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

         Zinterp          Zexact      Error

    0.629679        0.626806        0.29E-02
    0.934364        0.934391        0.27E-04
    0.892032        0.892091        0.59E-04
    0.968055        0.968107        0.53E-04
    0.673496        0.673147        0.35E-03

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

  Level     Order   Average Error

      0         1    0.37E-01
      1        27    0.56E-02
      2       125    0.32E-03
      3       729    0.83E-05
      4      4913    0.14E-06
      5     35937    0.59E-09

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.91E-01
      1      2      1      1    0.73E-01
      2      1      2      1    0.83E-01
      3      1      1      2    0.93E-01
      4      4      2      2    0.36E-01
      5      2      4      2    0.63E-01
      6      2      2      4    0.80E-01
      7      2      2      4    0.80E-01
      8      6      4      4    0.78E-03
      9      4      6      4    0.19E-15
     10      4      4      6    0.78E-03

TEST:
  Normal end of execution.

15 September 2021   8:27:05.600 AM