14 September 2021 7:36:19.470 PM PWL_INTERP_2D_SCATTERED_TEST: FORTRAN90 version Test PWL_INTERP_2D_SCATTERED. The R8LIB library is needed. This test also needs the TEST_INTERP_2D library. TEST01 R8TRIS2 computes the Delaunay triangulation of a set of nodes in 2D. TRIANGULATION_ORDER3_PRINT Information defining an order3 triangulation. The number of nodes is 9 Node coordinates Row 1 2 Col 1: 0.00000 0.00000 2: 0.00000 1.00000 3: 0.200000 0.500000 4: 0.300000 0.600000 5: 0.400000 0.500000 6: 0.600000 0.400000 7: 0.600000 0.500000 8: 1.00000 0.00000 9: 1.00000 1.00000 The number of triangles is 12 Sets of three nodes are used as vertices of the triangles. For each triangle, the nodes are listed in counterclockwise order. Triangle nodes: Row 1 2 3 Col 1 2 1 3 2 3 1 5 3 2 3 4 4 4 3 5 5 6 7 5 6 5 1 6 7 7 4 5 8 9 4 7 9 6 1 8 10 7 6 8 11 7 8 9 12 2 4 9 On each side of a given triangle, there is either another triangle, or a piece of the convex hull. For each triangle, we list the indices of the three neighbors, or (if negative) the codes of the segments of the convex hull. Triangle neighbors Row 1 2 3 Col 1 -28 2 3 2 1 6 4 3 1 4 12 4 3 2 7 5 10 7 6 6 2 9 5 7 8 4 5 8 12 7 11 9 6 -34 10 10 5 9 11 11 10 -38 8 12 3 8 -3 The number of boundary points is 4 The segments that make up the convex hull can be determined from the negative entries of the triangle neighbor list. # Tri Side N1 N2 1 9 2 1 8 2 11 2 8 9 3 12 3 9 2 4 1 1 2 1 TEST02 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.0000 0.0000 0.0000 0.0000 2 0.0000 0.2500 0.5000 0.5000 3 0.0000 0.5000 1.0000 1.0000 4 0.0000 0.7500 1.5000 1.5000 5 0.0000 1.0000 2.0000 2.0000 6 0.2500 0.0000 0.2500 0.2500 7 0.2500 0.2500 0.7500 0.7500 8 0.2500 0.5000 1.2500 1.2500 9 0.2500 0.7500 1.7500 1.7500 10 0.2500 1.0000 2.2500 2.2500 11 0.5000 0.0000 0.5000 0.5000 12 0.5000 0.2500 1.0000 1.0000 13 0.5000 0.5000 1.5000 1.5000 14 0.5000 0.7500 2.0000 2.0000 15 0.5000 1.0000 2.5000 2.5000 16 0.7500 0.0000 0.7500 0.7500 17 0.7500 0.2500 1.2500 1.2500 18 0.7500 0.5000 1.7500 1.7500 19 0.7500 0.7500 2.2500 2.2500 20 0.7500 1.0000 2.7500 2.7500 21 1.0000 0.0000 1.0000 1.0000 22 1.0000 0.2500 1.5000 1.5000 23 1.0000 0.5000 2.0000 2.0000 24 1.0000 0.7500 2.5000 2.5000 25 1.0000 1.0000 3.0000 3.0000 TEST03 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. Here, we use grid number 2 with 33 scattered points in the unit square on problem 1 RMS error is 0.65E-01 K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.1000 0.1000 0.9857 0.9857 2 0.1000 0.3000 0.8957 0.9716 3 0.1000 0.5000 0.5653 0.5182 4 0.1000 0.7000 0.3714 0.3412 5 0.1000 0.9000 0.2618 0.2805 6 0.3000 0.1000 0.9258 0.9610 7 0.3000 0.3000 0.9385 0.9836 8 0.3000 0.5000 0.6485 0.4690 9 0.3000 0.7000 0.3716 0.2576 10 0.3000 0.9000 0.1794 0.2174 11 0.5000 0.1000 0.5476 0.4855 12 0.5000 0.3000 0.6090 0.5210 13 0.5000 0.5000 0.4525 0.3258 14 0.5000 0.7000 0.2178 0.1080 15 0.5000 0.9000 0.1380 0.1166 16 0.7000 0.1000 0.3886 0.3614 17 0.7000 0.3000 0.5806 0.6136 18 0.7000 0.5000 0.4127 0.3994 19 0.7000 0.7000 0.1619 0.1503 20 0.7000 0.9000 0.1021 0.1021 21 0.9000 0.1000 0.2848 0.2372 22 0.9000 0.3000 0.4360 0.4569 23 0.9000 0.5000 0.2846 0.2904 24 0.9000 0.7000 0.1079 0.0910 25 0.9000 0.9000 0.0583 0.0563 TEST03 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. Here, we use grid number 2 with 33 scattered points in the unit square on problem 2 RMS error is 0.21E-01 K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.1000 0.1000 0.1147 0.1111 2 0.1000 0.3000 0.1961 0.2163 3 0.1000 0.5000 0.2121 0.2221 4 0.1000 0.7000 0.2202 0.2222 5 0.1000 0.9000 0.2222 0.2222 6 0.3000 0.1000 0.0319 0.0059 7 0.3000 0.3000 0.1220 0.1111 8 0.3000 0.5000 0.1781 0.2163 9 0.3000 0.7000 0.2048 0.2221 10 0.3000 0.9000 0.2222 0.2222 11 0.5000 0.1000 0.0079 0.0002 12 0.5000 0.3000 0.0595 0.0059 13 0.5000 0.5000 0.1383 0.1111 14 0.5000 0.7000 0.1794 0.2163 15 0.5000 0.9000 0.2216 0.2221 16 0.7000 0.1000 0.0000 0.0000 17 0.7000 0.3000 0.0002 0.0002 18 0.7000 0.5000 0.0556 0.0059 19 0.7000 0.7000 0.1083 0.1111 20 0.7000 0.9000 0.2163 0.2163 21 0.9000 0.1000 0.0000 0.0000 22 0.9000 0.3000 0.0000 0.0000 23 0.9000 0.5000 0.0029 0.0002 24 0.9000 0.7000 0.0206 0.0059 25 0.9000 0.9000 0.0946 0.1111 TEST03 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. Here, we use grid number 2 with 33 scattered points in the unit square on problem 3 RMS error is 0.34E-01 K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.1000 0.1000 0.2311 0.2358 2 0.1000 0.3000 0.1394 0.1343 3 0.1000 0.5000 0.0695 0.0387 4 0.1000 0.7000 0.0707 0.0500 5 0.1000 0.9000 0.1396 0.1563 6 0.3000 0.1000 0.2926 0.3478 7 0.3000 0.3000 0.1889 0.1982 8 0.3000 0.5000 0.1575 0.0571 9 0.3000 0.7000 0.1616 0.0738 10 0.3000 0.9000 0.1899 0.2305 11 0.5000 0.1000 0.2507 0.2810 12 0.5000 0.3000 0.1461 0.1601 13 0.5000 0.5000 0.0876 0.0461 14 0.5000 0.7000 0.0908 0.0596 15 0.5000 0.9000 0.1762 0.1863 16 0.7000 0.1000 0.1630 0.1590 17 0.7000 0.3000 0.0960 0.0906 18 0.7000 0.5000 0.0445 0.0261 19 0.7000 0.7000 0.0369 0.0337 20 0.7000 0.9000 0.1054 0.1054 21 0.9000 0.1000 0.0871 0.0903 22 0.9000 0.3000 0.0491 0.0514 23 0.9000 0.5000 0.0211 0.0148 24 0.9000 0.7000 0.0272 0.0191 25 0.9000 0.9000 0.0632 0.0599 TEST03 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. Here, we use grid number 2 with 33 scattered points in the unit square on problem 4 RMS error is 0.24E-01 K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.1000 0.1000 0.0700 0.0660 2 0.1000 0.3000 0.1148 0.1211 3 0.1000 0.5000 0.1354 0.1483 4 0.1000 0.7000 0.1135 0.1211 5 0.1000 0.9000 0.0766 0.0660 6 0.3000 0.1000 0.1112 0.1211 7 0.3000 0.3000 0.2157 0.2223 8 0.3000 0.5000 0.2138 0.2722 9 0.3000 0.7000 0.1749 0.2223 10 0.3000 0.9000 0.1245 0.1211 11 0.5000 0.1000 0.1428 0.1483 12 0.5000 0.3000 0.2371 0.2722 13 0.5000 0.5000 0.2630 0.3333 14 0.5000 0.7000 0.2396 0.2722 15 0.5000 0.9000 0.1411 0.1483 16 0.7000 0.1000 0.1165 0.1211 17 0.7000 0.3000 0.2145 0.2223 18 0.7000 0.5000 0.2398 0.2722 19 0.7000 0.7000 0.2205 0.2223 20 0.7000 0.9000 0.1211 0.1211 21 0.9000 0.1000 0.0766 0.0660 22 0.9000 0.3000 0.1224 0.1211 23 0.9000 0.5000 0.1450 0.1483 24 0.9000 0.7000 0.1150 0.1211 25 0.9000 0.9000 0.0668 0.0660 TEST03 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. Here, we use grid number 2 with 33 scattered points in the unit square on problem 5 RMS error is 0.47E-01 K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.1000 0.1000 0.0011 0.0005 2 0.1000 0.3000 0.0063 0.0058 3 0.1000 0.5000 0.0188 0.0131 4 0.1000 0.7000 0.0067 0.0058 5 0.1000 0.9000 0.0047 0.0005 6 0.3000 0.1000 0.0058 0.0058 7 0.3000 0.3000 0.0695 0.0660 8 0.3000 0.5000 0.0693 0.1483 9 0.3000 0.7000 0.0363 0.0660 10 0.3000 0.9000 0.0121 0.0058 11 0.5000 0.1000 0.0263 0.0131 12 0.5000 0.3000 0.0865 0.1483 13 0.5000 0.5000 0.1344 0.3333 14 0.5000 0.7000 0.1166 0.1483 15 0.5000 0.9000 0.0125 0.0131 16 0.7000 0.1000 0.0068 0.0058 17 0.7000 0.3000 0.0594 0.0660 18 0.7000 0.5000 0.0965 0.1483 19 0.7000 0.7000 0.0768 0.0660 20 0.7000 0.9000 0.0058 0.0058 21 0.9000 0.1000 0.0033 0.0005 22 0.9000 0.3000 0.0075 0.0058 23 0.9000 0.5000 0.0211 0.0131 24 0.9000 0.7000 0.0092 0.0058 25 0.9000 0.9000 0.0013 0.0005 TEST03 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. Here, we use grid number 2 with 33 scattered points in the unit square on problem 6 RMS error is 0.18E-01 K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.1000 0.1000 0.1746 0.1857 2 0.1000 0.3000 0.2565 0.2682 3 0.1000 0.5000 0.2762 0.2938 4 0.1000 0.7000 0.2580 0.2682 5 0.1000 0.9000 0.1526 0.1857 6 0.3000 0.1000 0.2531 0.2682 7 0.3000 0.3000 0.3364 0.3427 8 0.3000 0.5000 0.3343 0.3661 9 0.3000 0.7000 0.3094 0.3427 10 0.3000 0.9000 0.2412 0.2682 11 0.5000 0.1000 0.2766 0.2938 12 0.5000 0.3000 0.3501 0.3661 13 0.5000 0.5000 0.3617 0.3889 14 0.5000 0.7000 0.3470 0.3661 15 0.5000 0.9000 0.2856 0.2938 16 0.7000 0.1000 0.2602 0.2682 17 0.7000 0.3000 0.3379 0.3427 18 0.7000 0.5000 0.3507 0.3661 19 0.7000 0.7000 0.3380 0.3427 20 0.7000 0.9000 0.2682 0.2682 21 0.9000 0.1000 0.1594 0.1857 22 0.9000 0.3000 0.2573 0.2682 23 0.9000 0.5000 0.2813 0.2938 24 0.9000 0.7000 0.2545 0.2682 25 0.9000 0.9000 0.1768 0.1857 TEST03 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. Here, we use grid number 2 with 33 scattered points in the unit square on problem 7 RMS error is 0.79E+00 K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.1000 0.1000 0.4882 1.0091 2 0.1000 0.3000 0.0632 0.4480 3 0.1000 0.5000 -0.3180 -0.5568 4 0.1000 0.7000 1.2926 1.3542 5 0.1000 0.9000 -0.3702 1.2287 6 0.3000 0.1000 -0.4943 -1.3706 7 0.3000 0.3000 0.8294 0.5039 8 0.3000 0.5000 0.8474 2.8962 9 0.3000 0.7000 -0.1054 -0.4376 10 0.3000 0.9000 -1.2731 -0.3886 11 0.5000 0.1000 0.5320 0.9568 12 0.5000 0.3000 1.8509 1.0776 13 0.5000 0.5000 0.6402 0.0545 14 0.5000 0.7000 -0.5243 0.0219 15 0.5000 0.9000 -0.7664 -0.7437 16 0.7000 0.1000 1.2361 1.9130 17 0.7000 0.3000 1.0933 1.0760 18 0.7000 0.5000 -0.1089 -1.7967 19 0.7000 0.7000 -0.2092 0.0082 20 0.7000 0.9000 0.6382 0.6382 21 0.9000 0.1000 0.3073 -0.7501 22 0.9000 0.3000 0.4824 0.1702 23 0.9000 0.5000 0.1845 0.7699 24 0.9000 0.7000 -0.6264 -1.1804 25 0.9000 0.9000 0.0840 0.2189 TEST03 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. Here, we use grid number 2 with 33 scattered points in the unit square on problem 8 RMS error is 0.43E+00 K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.1000 0.1000 0.0038 0.0006 2 0.1000 0.3000 0.2559 0.1019 3 0.1000 0.5000 0.4853 0.7506 4 0.1000 0.7000 0.1237 0.1019 5 0.1000 0.9000 0.0768 0.0006 6 0.3000 0.1000 0.3056 0.1356 7 0.3000 0.3000 0.3164 0.2506 8 0.3000 0.5000 0.3669 0.9868 9 0.3000 0.7000 0.3070 0.2506 10 0.3000 0.9000 0.3512 0.1356 11 0.5000 0.1000 0.7820 1.0005 12 0.5000 0.3000 0.5910 1.2030 13 0.5000 0.5000 0.7601 2.5000 14 0.5000 0.7000 0.7415 1.2030 15 0.5000 0.9000 0.6612 1.0005 16 0.7000 0.1000 0.2203 0.1356 17 0.7000 0.3000 0.3795 0.2506 18 0.7000 0.5000 0.5014 0.9868 19 0.7000 0.7000 0.3480 0.2506 20 0.7000 0.9000 0.1356 0.1356 21 0.9000 0.1000 0.0204 0.0006 22 0.9000 0.3000 0.1838 0.1019 23 0.9000 0.5000 0.5678 0.7506 24 0.9000 0.7000 0.2287 0.1019 25 0.9000 0.9000 0.0274 0.0006 TEST03 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. Here, we use grid number 2 with 33 scattered points in the unit square on problem 9 RMS error is 0.17E+02 K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.1000 0.1000 1.6034 1.0265 2 0.1000 0.3000 4.6348 5.7630 3 0.1000 0.5000 5.9660 0.0000 4 0.1000 0.7000 -4.1937 -8.5418 5 0.1000 0.9000 -6.1631 -1.8651 6 0.3000 0.1000 4.5365 5.7630 7 0.3000 0.3000 26.6637 32.3542 8 0.3000 0.5000 20.7600 0.0000 9 0.3000 0.7000 -2.4982 -47.9552 10 0.3000 0.9000 -13.5956 -10.4712 11 0.5000 0.1000 -8.6655 0.0000 12 0.5000 0.3000 -5.0148 0.0000 13 0.5000 0.5000 48.0306 0.0000 14 0.5000 0.7000 39.8545 -0.0000 15 0.5000 0.9000 1.1601 -0.0000 16 0.7000 0.1000 -8.4011 -8.5418 17 0.7000 0.3000 -38.3621 -47.9552 18 0.7000 0.5000 15.2016 -0.0000 19 0.7000 0.7000 65.8221 71.0789 20 0.7000 0.9000 15.5203 15.5203 21 0.9000 0.1000 -5.5072 -1.8651 22 0.9000 0.3000 -11.3125 -10.4712 23 0.9000 0.5000 -0.4242 -0.0000 24 0.9000 0.7000 14.3110 15.5203 25 0.9000 0.9000 4.6781 3.3889 TEST03 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. Here, we use grid number 2 with 33 scattered points in the unit square on problem 10 RMS error is 0.29E+00 K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.1000 0.1000 0.0972 0.0857 2 0.1000 0.3000 0.1262 0.1643 3 0.1000 0.5000 0.0718 0.0243 4 0.1000 0.7000 0.1675 0.1643 5 0.1000 0.9000 0.0279 0.0857 6 0.3000 0.1000 0.1405 0.1926 7 0.3000 0.3000 -0.2736 -0.3401 8 0.3000 0.5000 -0.2637 -0.3888 9 0.3000 0.7000 -0.0383 -0.3401 10 0.3000 0.9000 0.0614 0.1926 11 0.5000 0.1000 -0.0103 0.1504 12 0.5000 0.3000 -0.3775 -0.4401 13 0.5000 0.5000 -0.3918 1.0000 14 0.5000 0.7000 -0.2250 -0.4401 15 0.5000 0.9000 0.1317 0.1504 16 0.7000 0.1000 0.1520 0.1926 17 0.7000 0.3000 -0.2929 -0.3401 18 0.7000 0.5000 -0.3715 -0.3888 19 0.7000 0.7000 -0.3087 -0.3401 20 0.7000 0.9000 0.1926 0.1926 21 0.9000 0.1000 0.0479 0.0857 22 0.9000 0.3000 0.1060 0.1643 23 0.9000 0.5000 -0.0217 0.0243 24 0.9000 0.7000 0.1066 0.1643 25 0.9000 0.9000 0.0610 0.0857 TEST03 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. Here, we use grid number 2 with 33 scattered points in the unit square on problem 11 RMS error is 0.50E-02 K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.1000 0.1000 0.1125 0.1100 2 0.1000 0.3000 0.1275 0.1300 3 0.1000 0.5000 0.1469 0.1500 4 0.1000 0.7000 0.1694 0.1700 5 0.1000 0.9000 0.1822 0.1900 6 0.3000 0.1000 0.3229 0.3300 7 0.3000 0.3000 0.3908 0.3900 8 0.3000 0.5000 0.4512 0.4500 9 0.3000 0.7000 0.5129 0.5100 10 0.3000 0.9000 0.5650 0.5700 11 0.5000 0.1000 0.5515 0.5500 12 0.5000 0.3000 0.6450 0.6500 13 0.5000 0.5000 0.7650 0.7500 14 0.5000 0.7000 0.8433 0.8500 15 0.5000 0.9000 0.9525 0.9500 16 0.7000 0.1000 0.7733 0.7700 17 0.7000 0.3000 0.9129 0.9100 18 0.7000 0.5000 1.0425 1.0500 19 0.7000 0.7000 1.1915 1.1900 20 0.7000 0.9000 1.3300 1.3300 21 0.9000 0.1000 0.9829 0.9900 22 0.9000 0.3000 1.1677 1.1700 23 0.9000 0.5000 1.3500 1.3500 24 0.9000 0.7000 1.5267 1.5300 25 0.9000 0.9000 1.7100 1.7100 TEST03 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. Here, we use grid number 2 with 33 scattered points in the unit square on problem 12 RMS error is 0.49E-01 K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.1000 0.1000 0.4167 0.4511 2 0.1000 0.3000 0.7905 0.8997 3 0.1000 0.5000 0.7571 0.8084 4 0.1000 0.7000 0.5331 0.5179 5 0.1000 0.9000 0.6090 0.5351 6 0.3000 0.1000 0.3824 0.3660 7 0.3000 0.3000 0.7416 0.7790 8 0.3000 0.5000 0.7244 0.6979 9 0.3000 0.7000 0.6278 0.4959 10 0.3000 0.9000 0.6843 0.6435 11 0.5000 0.1000 0.1755 0.2100 12 0.5000 0.3000 0.4733 0.4633 13 0.5000 0.5000 0.5586 0.4604 14 0.5000 0.7000 0.5280 0.4830 15 0.5000 0.9000 0.8826 0.8459 16 0.7000 0.1000 0.0877 0.0990 17 0.7000 0.3000 0.2633 0.2495 18 0.7000 0.5000 0.3861 0.3691 19 0.7000 0.7000 0.6628 0.6349 20 0.7000 0.9000 1.1678 1.1678 21 0.9000 0.1000 0.0798 0.0810 22 0.9000 0.3000 0.2659 0.2618 23 0.9000 0.5000 0.5606 0.5274 24 0.9000 0.7000 0.9797 0.9556 25 0.9000 0.9000 1.5049 1.5179 TEST03 PWL_INTERP_2D_SCATTERED_VALUE evaluates a piecewise linear interpolant to scattered data. Here, we use grid number 2 with 33 scattered points in the unit square on problem 13 RMS error is 0.16E+00 K Xi(K) Yi(K) Zi(K) Z(X,Y) 1 0.1000 0.1000 0.0316 0.0303 2 0.1000 0.3000 0.0462 0.0476 3 0.1000 0.5000 0.0600 0.0588 4 0.1000 0.7000 0.0464 0.0476 5 0.1000 0.9000 0.0354 0.0303 6 0.3000 0.1000 0.0451 0.0476 7 0.3000 0.3000 0.1142 0.1111 8 0.3000 0.5000 0.1136 0.2000 9 0.3000 0.7000 0.0813 0.1111 10 0.3000 0.9000 0.0521 0.0476 11 0.5000 0.1000 0.0663 0.0588 12 0.5000 0.3000 0.1315 0.2000 13 0.5000 0.5000 0.1886 1.0000 14 0.5000 0.7000 0.1730 0.2000 15 0.5000 0.9000 0.0566 0.0588 16 0.7000 0.1000 0.0470 0.0476 17 0.7000 0.3000 0.1050 0.1111 18 0.7000 0.5000 0.1459 0.2000 19 0.7000 0.7000 0.1206 0.1111 20 0.7000 0.9000 0.0476 0.0476 21 0.9000 0.1000 0.0345 0.0303 22 0.9000 0.3000 0.0489 0.0476 23 0.9000 0.5000 0.0628 0.0588 24 0.9000 0.7000 0.0481 0.0476 25 0.9000 0.9000 0.0310 0.0303 PWL_INTERP_2D_SCATTERED_TEST: Normal end of execution. 14 September 2021 7:36:19.474 PM