15 September 2021 3:23:18.932 PM SHEPARD_INTERP_1D_TEST: FORTRAN90 version Test the SHEPARD_INTERP_1D library. The R8LIB library is needed. This test needs the TEST_INTERP library as well. SHEPARD_BASIS_1D_TEST: SHEPARD_BASIS_1D evaluates the Shepard 1D basis functions. Using power P = 2.00000 The Shepard basis functions: Col 1 2 3 4 Row 1: 1.00000 0.00000 0.00000 0.00000 2: 0.887921 0.986578E-01 0.109620E-01 0.245961E-02 3: 0.481964 0.481964 0.301227E-01 0.595017E-02 4: 0.978970E-01 0.881073 0.179811E-01 0.304870E-02 5: 0.00000 1.00000 0.00000 0.00000 6: 0.368852E-01 0.922131 0.368852E-01 0.409836E-02 7: 0.804268E-01 0.723841 0.180960 0.147723E-01 8: 0.821097E-01 0.447042 0.447042 0.238070E-01 9: 0.466321E-01 0.186528 0.746114 0.207254E-01 10: 0.116402E-01 0.377141E-01 0.942854 0.779218E-02 11: 0.00000 0.00000 1.00000 0.00000 12: 0.793883E-02 0.196040E-01 0.960598 0.118592E-01 13: 0.240964E-01 0.542169E-01 0.867470 0.542169E-01 14: 0.395051E-01 0.824243E-01 0.741818 0.136252 15: 0.484157E-01 0.948948E-01 0.593093 0.263597 16: 0.479398E-01 0.891442E-01 0.431458 0.431458 17: 0.386266E-01 0.686695E-01 0.274678 0.618026 18: 0.245585E-01 0.419965E-01 0.144845 0.788600 19: 0.112720E-01 0.186333E-01 0.570644E-01 0.913030 20: 0.271694E-02 0.435918E-02 0.121088E-01 0.980815 21: 0.00000 0.00000 0.00000 1.00000 SHEPARD_VALUE_1D_TEST: SHEPARD_VALUE_1D evaluates a Shepard 1D interpolant. Using power P = 2.00000 Table of interpolant values: 1 0.00000 -28.0000 2 0.500000 -23.5644 3 1.00000 -7.95240 4 1.50000 6.41966 5 2.00000 10.0000 6 2.50000 8.55738 7 3.00000 5.93188 8 3.50000 2.45174 9 4.00000 -2.05181 10 4.50000 -6.22928 11 5.00000 -8.00000 12 5.50000 -5.78983 13 6.00000 1.71084 14 6.50000 15.8564 15 7.00000 37.5512 16 7.50000 65.9937 17 8.00000 97.5279 18 8.50000 126.327 19 9.00000 147.325 20 9.50000 158.763 21 10.0000 162.000 TEST01: Interpolate data from TEST_INTERP problem # 1 using Shepard interpolation with P = 0.00000 Number of data points = 18 Data array: Row 1 2 Col 1: 0.00000 4.00000 2: 1.00000 5.00000 3: 2.00000 6.00000 4: 4.00000 6.00000 5: 5.00000 5.00000 6: 6.00000 3.00000 7: 7.00000 1.00000 8: 8.00000 1.00000 9: 9.00000 1.00000 10: 10.0000 3.00000 11: 11.0000 4.00000 12: 12.0000 4.00000 13: 13.0000 3.00000 14: 14.0000 3.00000 15: 15.0000 4.00000 16: 16.0000 4.00000 17: 17.0000 3.00000 18: 18.0000 0.00000 L2 interpolation error averaged per interpolant node = 0.392837 Normalized length of piecewise linear interpolant = 3.13553 Normalized length of Shepard interpolant = 1.00000 TEST01: Interpolate data from TEST_INTERP problem # 1 using Shepard interpolation with P = 1.00000 Number of data points = 18 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 3.13553 Normalized length of Shepard interpolant = 4.81280 TEST01: Interpolate data from TEST_INTERP problem # 1 using Shepard interpolation with P = 2.00000 Number of data points = 18 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 3.13553 Normalized length of Shepard interpolant = 3.59344 TEST01: Interpolate data from TEST_INTERP problem # 1 using Shepard interpolation with P = 4.00000 Number of data points = 18 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 3.13553 Normalized length of Shepard interpolant = 3.36283 TEST01: Interpolate data from TEST_INTERP problem # 1 using Shepard interpolation with P = 8.00000 Number of data points = 18 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 3.13553 Normalized length of Shepard interpolant = 3.47367 TEST01: Interpolate data from TEST_INTERP problem # 2 using Shepard interpolation with P = 0.00000 Number of data points = 18 Data array: Row 1 2 Col 1: 0.00000 0.00000 2: 1.34000 5.00000 3: 5.00000 8.66000 4: 10.0000 10.0000 5: 10.6000 10.4000 6: 10.7000 12.0000 7: 10.7050 28.6000 8: 10.8000 30.2000 9: 11.4000 30.6000 10: 19.6000 30.6000 11: 20.2000 30.2000 12: 20.2950 28.6000 13: 20.3000 12.0000 14: 20.4000 10.4000 15: 21.0000 10.0000 16: 26.0000 8.66000 17: 29.6600 5.00000 18: 31.0000 0.00000 L2 interpolation error averaged per interpolant node = 2.58162 Normalized length of piecewise linear interpolant = 2.66086 Normalized length of Shepard interpolant = 1.00000 TEST01: Interpolate data from TEST_INTERP problem # 2 using Shepard interpolation with P = 1.00000 Number of data points = 18 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.66086 Normalized length of Shepard interpolant = 5.98166 TEST01: Interpolate data from TEST_INTERP problem # 2 using Shepard interpolation with P = 2.00000 Number of data points = 18 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.66086 Normalized length of Shepard interpolant = 5.44204 TEST01: Interpolate data from TEST_INTERP problem # 2 using Shepard interpolation with P = 4.00000 Number of data points = 18 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.66086 Normalized length of Shepard interpolant = 4.62358 TEST01: Interpolate data from TEST_INTERP problem # 2 using Shepard interpolation with P = 8.00000 Number of data points = 18 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.66086 Normalized length of Shepard interpolant = 3.57979 TEST01: Interpolate data from TEST_INTERP problem # 3 using Shepard interpolation with P = 0.00000 Number of data points = 11 Data array: Row 1 2 Col 1: 0.00000 0.00000 2: 2.00000 10.0000 3: 3.00000 10.0000 4: 5.00000 10.0000 5: 6.00000 10.0000 6: 8.00000 10.0000 7: 9.00000 10.5000 8: 11.0000 15.0000 9: 12.0000 50.0000 10: 14.0000 60.0000 11: 15.0000 85.0000 L2 interpolation error averaged per interpolant node = 7.87737 Normalized length of piecewise linear interpolant = 1.68472 Normalized length of Shepard interpolant = 1.00000 TEST01: Interpolate data from TEST_INTERP problem # 3 using Shepard interpolation with P = 1.00000 Number of data points = 11 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 1.68472 Normalized length of Shepard interpolant = 2.41306 TEST01: Interpolate data from TEST_INTERP problem # 3 using Shepard interpolation with P = 2.00000 Number of data points = 11 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 1.68472 Normalized length of Shepard interpolant = 1.84529 TEST01: Interpolate data from TEST_INTERP problem # 3 using Shepard interpolation with P = 4.00000 Number of data points = 11 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 1.68472 Normalized length of Shepard interpolant = 1.81905 TEST01: Interpolate data from TEST_INTERP problem # 3 using Shepard interpolation with P = 8.00000 Number of data points = 11 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 1.68472 Normalized length of Shepard interpolant = 1.87204 TEST01: Interpolate data from TEST_INTERP problem # 4 using Shepard interpolation with P = 0.00000 Number of data points = 8 Data array: Row 1 2 Col 1: 0.00000 0.00000 2: 0.500000E-01 0.700000 3: 0.100000 1.00000 4: 0.200000 1.00000 5: 0.800000 0.300000 6: 0.850000 0.500000E-01 7: 0.900000 0.100000 8: 1.00000 1.00000 L2 interpolation error averaged per interpolant node = 0.149984 Normalized length of piecewise linear interpolant = 3.25908 Normalized length of Shepard interpolant = 1.00000 TEST01: Interpolate data from TEST_INTERP problem # 4 using Shepard interpolation with P = 1.00000 Number of data points = 8 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 3.25908 Normalized length of Shepard interpolant = 3.92156 TEST01: Interpolate data from TEST_INTERP problem # 4 using Shepard interpolation with P = 2.00000 Number of data points = 8 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 3.25908 Normalized length of Shepard interpolant = 3.48948 TEST01: Interpolate data from TEST_INTERP problem # 4 using Shepard interpolation with P = 4.00000 Number of data points = 8 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 3.25908 Normalized length of Shepard interpolant = 3.54049 TEST01: Interpolate data from TEST_INTERP problem # 4 using Shepard interpolation with P = 8.00000 Number of data points = 8 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 3.25908 Normalized length of Shepard interpolant = 3.69223 TEST01: Interpolate data from TEST_INTERP problem # 5 using Shepard interpolation with P = 0.00000 Number of data points = 9 Data array: Row 1 2 Col 1: 0.00000 0.00000 2: 0.100000 0.900000 3: 0.200000 0.950000 4: 0.300000 0.900000 5: 0.400000 0.100000 6: 0.500000 0.500000E-01 7: 0.600000 0.500000E-01 8: 0.800000 0.200000 9: 1.00000 1.00000 L2 interpolation error averaged per interpolant node = 0.143324 Normalized length of piecewise linear interpolant = 3.22180 Normalized length of Shepard interpolant = 1.00000 TEST01: Interpolate data from TEST_INTERP problem # 5 using Shepard interpolation with P = 1.00000 Number of data points = 9 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 3.22180 Normalized length of Shepard interpolant = 4.81294 TEST01: Interpolate data from TEST_INTERP problem # 5 using Shepard interpolation with P = 2.00000 Number of data points = 9 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 3.22180 Normalized length of Shepard interpolant = 3.58093 TEST01: Interpolate data from TEST_INTERP problem # 5 using Shepard interpolation with P = 4.00000 Number of data points = 9 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 3.22180 Normalized length of Shepard interpolant = 3.44206 TEST01: Interpolate data from TEST_INTERP problem # 5 using Shepard interpolation with P = 8.00000 Number of data points = 9 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 3.22180 Normalized length of Shepard interpolant = 3.57007 TEST01: Interpolate data from TEST_INTERP problem # 6 using Shepard interpolation with P = 0.00000 Number of data points = 49 Data array: Row 1 2 Col 1: 595.000 0.644000 2: 605.000 0.622000 3: 615.000 0.638000 4: 625.000 0.649000 5: 635.000 0.652000 6: 645.000 0.639000 7: 655.000 0.646000 8: 665.000 0.657000 9: 675.000 0.652000 10: 685.000 0.655000 11: 695.000 0.644000 12: 705.000 0.663000 13: 715.000 0.663000 14: 725.000 0.668000 15: 735.000 0.676000 16: 745.000 0.676000 17: 755.000 0.686000 18: 765.000 0.679000 19: 775.000 0.678000 20: 785.000 0.683000 21: 795.000 0.694000 22: 805.000 0.699000 23: 815.000 0.710000 24: 825.000 0.730000 25: 835.000 0.763000 26: 845.000 0.812000 27: 855.000 0.907000 28: 865.000 1.04400 29: 875.000 1.33600 30: 885.000 1.88100 31: 895.000 2.16900 32: 905.000 2.07500 33: 915.000 1.59800 34: 925.000 1.21100 35: 935.000 0.916000 36: 945.000 0.746000 37: 955.000 0.672000 38: 965.000 0.627000 39: 975.000 0.615000 40: 985.000 0.607000 41: 995.000 0.606000 42: 1005.00 0.609000 43: 1015.00 0.603000 44: 1025.00 0.601000 45: 1035.00 0.603000 46: 1045.00 0.601000 47: 1055.00 0.611000 48: 1065.00 0.601000 49: 1075.00 0.608000 L2 interpolation error averaged per interpolant node = 0.530251E-01 Normalized length of piecewise linear interpolant = 2.68209 Normalized length of Shepard interpolant = 1.00000 TEST01: Interpolate data from TEST_INTERP problem # 6 using Shepard interpolation with P = 1.00000 Number of data points = 49 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.68209 Normalized length of Shepard interpolant = 5.91812 TEST01: Interpolate data from TEST_INTERP problem # 6 using Shepard interpolation with P = 2.00000 Number of data points = 49 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.68209 Normalized length of Shepard interpolant = 3.00366 TEST01: Interpolate data from TEST_INTERP problem # 6 using Shepard interpolation with P = 4.00000 Number of data points = 49 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.68209 Normalized length of Shepard interpolant = 2.80798 TEST01: Interpolate data from TEST_INTERP problem # 6 using Shepard interpolation with P = 8.00000 Number of data points = 49 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.68209 Normalized length of Shepard interpolant = 2.88659 TEST01: Interpolate data from TEST_INTERP problem # 7 using Shepard interpolation with P = 0.00000 Number of data points = 4 Data array: Row 1 2 Col 1: 0.00000 1.00000 2: 1.00000 2.00000 3: 4.00000 2.00000 4: 5.00000 1.00000 L2 interpolation error averaged per interpolant node = 0.250000 Normalized length of piecewise linear interpolant = 2.63961 Normalized length of Shepard interpolant = 1.00000 TEST01: Interpolate data from TEST_INTERP problem # 7 using Shepard interpolation with P = 1.00000 Number of data points = 4 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.63961 Normalized length of Shepard interpolant = 3.10373 TEST01: Interpolate data from TEST_INTERP problem # 7 using Shepard interpolation with P = 2.00000 Number of data points = 4 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.63961 Normalized length of Shepard interpolant = 2.91001 TEST01: Interpolate data from TEST_INTERP problem # 7 using Shepard interpolation with P = 4.00000 Number of data points = 4 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.63961 Normalized length of Shepard interpolant = 2.84058 TEST01: Interpolate data from TEST_INTERP problem # 7 using Shepard interpolation with P = 8.00000 Number of data points = 4 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.63961 Normalized length of Shepard interpolant = 2.87137 TEST01: Interpolate data from TEST_INTERP problem # 8 using Shepard interpolation with P = 0.00000 Number of data points = 12 Data array: Row 1 2 Col 1: -1.00000 1.00000 2: -0.800000 0.640000 3: -0.600000 0.360000 4: -0.400000 0.160000 5: -0.200000 0.400000E-01 6: 0.00000 0.00000 7: 0.200000 0.400000E-01 8: 0.200010 0.500000E-01 9: 0.400000 0.160000 10: 0.600000 0.360000 11: 0.800000 0.640000 12: 1.00000 1.00000 L2 interpolation error averaged per interpolant node = 0.101554 Normalized length of piecewise linear interpolant = 2.31939 Normalized length of Shepard interpolant = 1.00000 TEST01: Interpolate data from TEST_INTERP problem # 8 using Shepard interpolation with P = 1.00000 Number of data points = 12 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.31939 Normalized length of Shepard interpolant = 3.26227 TEST01: Interpolate data from TEST_INTERP problem # 8 using Shepard interpolation with P = 2.00000 Number of data points = 12 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.31939 Normalized length of Shepard interpolant = 2.48289 TEST01: Interpolate data from TEST_INTERP problem # 8 using Shepard interpolation with P = 4.00000 Number of data points = 12 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.31939 Normalized length of Shepard interpolant = 2.60246 TEST01: Interpolate data from TEST_INTERP problem # 8 using Shepard interpolation with P = 8.00000 Number of data points = 12 L2 interpolation error averaged per interpolant node = 0.00000 Normalized length of piecewise linear interpolant = 2.31939 Normalized length of Shepard interpolant = 2.75121 SHEPARD_INTERP_1D_TEST: Normal end of execution. 15 September 2021 3:23:18.942 PM