Tue Oct 19 12:01:07 2021 newton_interp_1d_test(): Python version: Test newton_interp_1d(). p00_data_test() tests p00_data p00_data() returns the actual (MxN) data for any problem. Problem 1 DATA_num = 18 DIM_num = 2 Data array: Row: 0 1 Col 0 : 0 4 1 : 1 5 2 : 2 6 3 : 4 6 4 : 5 5 5 : 6 3 6 : 7 1 7 : 8 1 8 : 9 1 9 : 10 3 10 : 11 4 11 : 12 4 12 : 13 3 13 : 14 3 14 : 15 4 15 : 16 4 16 : 17 3 17 : 18 0 Problem 2 DATA_num = 18 DIM_num = 2 Data array: Row: 0 1 Col 0 : 0 0 1 : 1.34 5 2 : 5 8.66 3 : 10 10 4 : 10.6 10.4 5 : 10.7 12 6 : 10.705 28.6 7 : 10.8 30.2 8 : 11.4 30.6 9 : 19.6 30.6 10 : 20.2 30.2 11 : 20.295 28.6 12 : 20.3 12 13 : 20.4 10.4 14 : 21 10 15 : 26 8.66 16 : 29.66 5 17 : 31 0 Problem 3 DATA_num = 11 DIM_num = 2 Data array: Row: 0 1 Col 0 : 0 0 1 : 2 10 2 : 3 10 3 : 5 10 4 : 6 10 5 : 8 10 6 : 9 10.5 7 : 11 15 8 : 12 50 9 : 14 60 10 : 15 85 Problem 4 DATA_num = 8 DIM_num = 2 Data array: Row: 0 1 Col 0 : 0 0 1 : 0.05 0.7 2 : 0.1 1 3 : 0.2 1 4 : 0.8 0.3 5 : 0.85 0.05 6 : 0.9 0.1 7 : 1 1 Problem 5 DATA_num = 9 DIM_num = 2 Data array: Row: 0 1 Col 0 : 0 0 1 : 0.1 0.9 2 : 0.2 0.95 3 : 0.3 0.9 4 : 0.4 0.1 5 : 0.5 0.05 6 : 0.6 0.05 7 : 0.8 0.2 8 : 1 1 Problem 6 DATA_num = 49 DIM_num = 2 Data array: Row: 0 1 Col 0 : 595 0.644 1 : 605 0.622 2 : 615 0.638 3 : 625 0.649 4 : 635 0.652 5 : 645 0.639 6 : 655 0.646 7 : 665 0.657 8 : 675 0.652 9 : 685 0.655 10 : 695 0.644 11 : 705 0.663 12 : 715 0.663 13 : 725 0.668 14 : 735 0.676 15 : 745 0.676 16 : 755 0.686 17 : 765 0.679 18 : 775 0.678 19 : 785 0.683 20 : 795 0.694 21 : 805 0.699 22 : 815 0.71 23 : 825 0.73 24 : 835 0.763 25 : 845 0.812 26 : 855 0.907 27 : 865 1.044 28 : 875 1.336 29 : 885 1.881 30 : 895 2.169 31 : 905 2.075 32 : 915 1.598 33 : 925 1.211 34 : 935 0.916 35 : 945 0.746 36 : 955 0.672 37 : 965 0.627 38 : 975 0.615 39 : 985 0.607 40 : 995 0.606 41 : 1005 0.609 42 : 1015 0.603 43 : 1025 0.601 44 : 1035 0.603 45 : 1045 0.601 46 : 1055 0.611 47 : 1065 0.601 48 : 1075 0.608 Problem 7 DATA_num = 4 DIM_num = 2 Data array: Row: 0 1 Col 0 : 0 1 1 : 1 2 2 : 4 2 3 : 5 1 Problem 8 DATA_num = 12 DIM_num = 2 Data array: Row: 0 1 Col 0 : -1 1 1 : -0.8 0.64 2 : -0.6 0.36 3 : -0.4 0.16 4 : -0.2 0.04 5 : 0 0 6 : 0.2 0.04 7 : 0.20001 0.05 8 : 0.4 0.16 9 : 0.6 0.36 10 : 0.8 0.64 11 : 1 1 p00_data_test: Normal end of execution. p00_data_num_test p00_data_num() returns the number of data points for any problem. Problem Data Num 1 18 2 18 3 11 4 8 5 9 6 49 7 4 8 12 p00_data_num_test: Normal end of execution. p00_dim_num_test p00_dim_num() returns the spatial dimension for any problem. Problem Dimension 1 2 2 2 3 2 4 2 5 2 6 2 7 2 8 2 p00_dim_num_test: Normal end of execution. p00_prob_num_test p00_prob_num() returns the number of test problems. TEST_interp includes 8 test problems. p00_prob_num_test: Normal end of execution. newton_coef_1d_test newton_coef_1d sets the coefficients for a 1D Newton interpolation. Interpolation data: 0: 0 24 1: 1 0 2: 2 0 3: 3 0 4: 4 0 Newton interpolant coefficients: 0: 24 1: -24 2: 12 3: -4 4: 1 newton_coef_1d_test: Normal end of execution. newton_value_1d_test newton_value_1d evaluates a Newton 1d interpolant. The Newton interpolant data: 0: 0 24 1: 1 -24 2: 2 12 3: 3 -4 4: 4 1 Newton interpolant values: 0: 0 24 1: 0.333333 10.8642 2: 0.666667 3.45679 3: 1 0 4: 1.33333 -0.987654 5: 1.66667 -0.691358 6: 2 0 7: 2.33333 0.493827 8: 2.66667 0.493827 9: 3 0 10: 3.33333 -0.691358 11: 3.66667 -0.987654 12: 4 0 13: 4.33333 3.45679 14: 4.66667 10.8642 15: 5 24 newton_value_1d_test: Normal end of execution. newton_interp_1d_test01: Interpolate data from TEST_interp problem #1. Number of data points = 18 L2 interpolation error averaged per interpolant node = 1.48478e-11 Normalized length of piecewise linear interpolant = 3.13553 Normalized length of Shepard interpolant = 23.6452 Created plot file "p01_data.png" Created plot file "p01_newton.png". newton_interp_1d_test01: Normal end of execution. newton_interp_1d_test01: Interpolate data from TEST_interp problem #2. Number of data points = 18 L2 interpolation error averaged per interpolant node = 0.00141825 Normalized length of piecewise linear interpolant = 2.66086 Normalized length of Shepard interpolant = 86277.7 Created plot file "p02_data.png" Created plot file "p02_newton.png". newton_interp_1d_test01: Normal end of execution. newton_interp_1d_test01: Interpolate data from TEST_interp problem #3. Number of data points = 11 L2 interpolation error averaged per interpolant node = 2.23156e-14 Normalized length of piecewise linear interpolant = 1.68472 Normalized length of Shepard interpolant = 2.07523 Created plot file "p03_data.png" Created plot file "p03_newton.png". newton_interp_1d_test01: Normal end of execution. newton_interp_1d_test01: Interpolate data from TEST_interp problem #4. Number of data points = 8 Data array: 0: 0 0 1: 0.05 0.7 2: 0.1 1 3: 0.2 1 4: 0.8 0.3 5: 0.85 0.05 6: 0.9 0.1 7: 1 1 L2 interpolation error averaged per interpolant node = 1.76905e-15 Normalized length of piecewise linear interpolant = 3.25908 Normalized length of Shepard interpolant = 2.31988 Created plot file "p04_data.png" Created plot file "p04_newton.png". newton_interp_1d_test01: Normal end of execution. newton_interp_1d_test01: Interpolate data from TEST_interp problem #5. Number of data points = 9 Data array: 0: 0 0 1: 0.1 0.9 2: 0.2 0.95 3: 0.3 0.9 4: 0.4 0.1 5: 0.5 0.05 6: 0.6 0.05 7: 0.8 0.2 8: 1 1 L2 interpolation error averaged per interpolant node = 6.19464e-15 Normalized length of piecewise linear interpolant = 3.2218 Normalized length of Shepard interpolant = 6.09886 Created plot file "p05_data.png" Created plot file "p05_newton.png". newton_interp_1d_test01: Normal end of execution. newton_interp_1d_test01: Interpolate data from TEST_interp problem #6. Number of data points = 49 L2 interpolation error averaged per interpolant node = 20.2923 Normalized length of piecewise linear interpolant = 2.68209 Normalized length of Shepard interpolant = 74265.4 Created plot file "p06_data.png" Created plot file "p06_newton.png". newton_interp_1d_test01: Normal end of execution. newton_interp_1d_test01: Interpolate data from TEST_interp problem #7. Number of data points = 4 Data array: 0: 0 1 1: 1 2 2: 4 2 3: 5 1 L2 interpolation error averaged per interpolant node = 0 Normalized length of piecewise linear interpolant = 2.63961 Normalized length of Shepard interpolant = 1.83502 Created plot file "p07_data.png" Created plot file "p07_newton.png". newton_interp_1d_test01: Normal end of execution. newton_interp_1d_test01: Interpolate data from TEST_interp problem #8. Number of data points = 12 L2 interpolation error averaged per interpolant node = 2.44769e-12 Normalized length of piecewise linear interpolant = 2.31939 Normalized length of Shepard interpolant = 152.207 Created plot file "p08_data.png" Created plot file "p08_newton.png". newton_interp_1d_test01: Normal end of execution. newton_interp_1d_test: Normal end of execution. Tue Oct 19 12:01:09 2021