Testing CURVE_HERMITE_INTERPOLATOR Here are stored abscissas, ordinates and computed derivatives: (Compare with derivatives given as input in NAG test e01bgc) Derivatives = BASIC_VECTOR <9|, real: Index # Values 1 0.0000e+000 5.5251e-004 3.3587e-001 3.4944e-001 5.9696e-001 6 6.0326e-002 8.9833e-004 2.9405e-005 0.0000e+000 Abscissas = BASIC_VECTOR <9|, real: Index # Values 1 7.9900e+000 8.0900e+000 8.1900e+000 8.7000e+000 9.2000e+000 6 1.0000e+001 1.2000e+001 1.5000e+001 2.0000e+001 Ordinates = BASIC_VECTOR <9|, real: Index # Values 1 0.0000e+000 2.7643e-005 4.3750e-002 1.6918e-001 4.6943e-001 6 9.4374e-001 9.9864e-001 9.9992e-001 9.9999e-001 Testing 'integral' for four limits - see NAG e01bhc a = 7.9900e+000 b = 2.0000e+001-->integral = 1.0765e+001 a = 1.0000e+001 b = 1.2000e+001-->integral = 1.9622e+000 a = 1.2000e+001 b = 1.0000e+001-->integral = -1.9622e+000 a = 1.5000e+001 b = 1.5000e+001-->integral = 0.0000e+000 Interposed abscissas, ordinates and derivatives (Compare with output in NAG test e01bgc) interposed_abscissas = BASIC_VECTOR <11|, real: Index # Values 1 7.9900e+000 9.1910e+000 1.0392e+001 1.1593e+001 1.2794e+001 6 1.3995e+001 1.5196e+001 1.6397e+001 1.7598e+001 1.8799e+001 11 2.0000e+001 interposed_derivatives = BASIC_VECTOR <11|, real: Index # Values 1 0.0000e+000 6.0601e-001 4.5688e-002 9.9166e-003 6.2491e-004 6 2.7077e-004 2.8094e-005 2.0341e-005 1.3074e-005 6.2938e-006 11 0.0000e+000 interposed_ordinates = BASIC_VECTOR <11|, real: Index # Values 1 0.0000e+000 4.6402e-001 9.6447e-001 9.9649e-001 9.9924e-001 6 9.9977e-001 9.9993e-001 9.9995e-001 9.9997e-001 9.9999e-001 11 9.9999e-001 End of tests for CURVE_HERMITE_INTERPOLATOR class Test of CURVE_SPLINE_INTERPOLATOR (eo1bac) x = BASIC_VECTOR <7|, real: Index # Values 1 0.0000e+000 2.0000e-001 4.0000e-001 6.0000e-001 7.5000e-001 6 9.0000e-001 1.0000e+000 y = BASIC_VECTOR <7|, real: Index # Values 1 1.0000e+000 1.2214e+000 1.4918e+000 1.8221e+000 2.1170e+000 6 2.4596e+000 2.7182e+000 At x = 0.0000e+000 value = 1.0000e+000 At x = 1.0000e-001 value = 1.1052e+000 At x = 2.0000e-001 value = 1.2214e+000 At x = 3.0000e-001 value = 1.3498e+000 At x = 4.0000e-001 value = 1.4918e+000 At x = 5.0000e-001 value = 1.6487e+000 At x = 6.0000e-001 value = 1.8221e+000 At x = 6.7500e-001 value = 1.9640e+000 At x = 7.5000e-001 value = 2.1170e+000 At x = 8.2500e-001 value = 2.2819e+000 At x = 9.0000e-001 value = 2.4596e+000 At x = 9.5000e-001 value = 2.5857e+000 At x = 1.0000e+000 value = 2.7182e+000 Test of value_with_derivatives At x = 0 [y, y', y'', y'''] = BASIC_VECTOR <4|, real: Index # Values 1 1.0000e+000 1.0017e+000 9.6661e-001 1.2944e+000 At x = 0.10000000000000001 [y, y', y'', y'''] = BASIC_VECTOR <4|, real: Index # Values 1 1.1052e+000 1.1049e+000 1.0960e+000 1.2944e+000 At x = 0.20000000000000001 [y, y', y'', y'''] = BASIC_VECTOR <4|, real: Index # Values 1 1.2214e+000 1.2209e+000 1.2255e+000 1.2944e+000 At x = 0.29999999999999999 [y, y', y'', y'''] = BASIC_VECTOR <4|, real: Index # Values 1 1.3498e+000 1.3500e+000 1.3549e+000 1.2944e+000 At x = 0.40000000000000002 [y, y', y'', y'''] = BASIC_VECTOR <4|, real: Index # Values 1 1.4918e+000 1.4919e+000 1.4843e+000 1.6681e+000 At x = 0.5 [y, y', y'', y'''] = BASIC_VECTOR <4|, real: Index # Values 1 1.6487e+000 1.6487e+000 1.6512e+000 1.6681e+000 At x = 0.59999999999999998 [y, y', y'', y'''] = BASIC_VECTOR <4|, real: Index # Values 1 1.8221e+000 1.8221e+000 1.8180e+000 1.9684e+000 At x = 0.67500000000000004 [y, y', y'', y'''] = BASIC_VECTOR <4|, real: Index # Values 1 1.9640e+000 1.9640e+000 1.9656e+000 1.9684e+000 At x = 0.75 [y, y', y'', y'''] = BASIC_VECTOR <4|, real: Index # Values 1 2.1170e+000 2.1170e+000 2.1132e+000 2.2786e+000 At x = 0.82499999999999996 [y, y', y'', y'''] = BASIC_VECTOR <4|, real: Index # Values 1 2.2819e+000 2.2819e+000 2.2841e+000 2.2786e+000 At x = 0.90000000000000002 [y, y', y'', y'''] = BASIC_VECTOR <4|, real: Index # Values 1 2.4596e+000 2.4596e+000 2.4550e+000 2.2786e+000 At x = 0.94999999999999996 [y, y', y'', y'''] = BASIC_VECTOR <4|, real: Index # Values 1 2.5857e+000 2.5852e+000 2.5689e+000 2.2786e+000 At x = 1 [y, y', y'', y'''] = BASIC_VECTOR <4|, real: Index # Values 1 2.7182e+000 2.7165e+000 2.6829e+000 2.2786e+000 Integral (should be e -1.) = 1.7182835320344108 End of tests for CURVE_SPLINE_INTERPOLATOR class Testing CURVE_SPLINE_APPROXIMATOR - NAG e02bec For smoothness = = 1 Coefficients = BASIC_VECTOR <5|, real: Index # Values 1 -1.3201e+000 1.3542e+000 5.5510e+000 4.7031e+000 8.2277e+000 Closeness = 1.0003358092819732 Knots = BASIC_VECTOR <9|, real: Index # Values 1 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 4.0000e+000 6 8.0000e+000 8.0000e+000 8.0000e+000 8.0000e+000 For smoothness = = 0.5 Coefficients = BASIC_VECTOR <9|, real: Index # Values 1 -1.1072e+000 -6.5715e-001 4.3500e-001 2.8061e+000 4.6824e+000 6 4.6416e+000 5.1976e+000 6.9008e+000 7.9979e+000 Closeness = 0.50010095149758516 Knots = BASIC_VECTOR <13|, real: Index # Values 1 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 1.0000e+000 6 2.0000e+000 4.0000e+000 5.0000e+000 6.0000e+000 8.0000e+000 11 8.0000e+000 8.0000e+000 8.0000e+000 For smoothness = = 0.10000000000000001 Coefficients = BASIC_VECTOR <12|, real: Index # Values 1 -1.0900e+000 -6.4217e-001 3.6930e-002 1.6353e+000 2.1274e+000 6 4.5526e+000 4.2225e+000 4.9108e+000 4.4159e+000 5.4794e+000 11 6.8308e+000 7.9935e+000 Closeness = 0.10000016449613151 Knots = BASIC_VECTOR <16|, real: Index # Values 1 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 1.0000e+000 6 1.5000e+000 2.0000e+000 3.0000e+000 4.0000e+000 4.5000e+000 11 5.0000e+000 6.0000e+000 8.0000e+000 8.0000e+000 8.0000e+000 16 8.0000e+000 For smoothness = = 0 Coefficients = BASIC_VECTOR <15|, real: Index # Values 1 -1.1000e+000 -4.0554e-001 -3.4241e-001 1.9989e+000 1.9016e+000 6 3.0546e+000 4.9548e+000 3.8575e+000 5.0606e+000 4.4592e+000 11 4.7625e+000 5.3869e+000 6.1556e+000 7.5481e+000 7.9700e+000 Closeness = 0 Knots = BASIC_VECTOR <19|, real: Index # Values 1 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 1.0000e+000 6 1.5000e+000 2.0000e+000 2.5000e+000 3.0000e+000 4.0000e+000 11 4.5000e+000 5.0000e+000 5.5000e+000 6.0000e+000 7.0000e+000 16 8.0000e+000 8.0000e+000 8.0000e+000 8.0000e+000 -------------------------------------------------- Comparing the results of evaluation 'value' (e02bbc) for smoothnesses: 1.0, 0.5, 0.1, 0.0 Refer to results of NAG test e02bec - shown as three plots there. x = 0.0000e+000 -> -1.3201e+000 -1.1072e+000 -1.0900e+000 -1.1000e+000 x = 2.5000e-001 -> -8.2561e-001 -7.5278e-001 -7.4968e-001 -6.7953e-001 x = 5.0000e-001 -> -3.4558e-001 -3.6918e-001 -3.8632e-001 -3.7200e-001 x = 7.5000e-001 -> 1.1893e-001 3.6799e-002 2.1268e-002 -4.5722e-002 x = 1.0000e+000 -> 5.6693e-001 4.5835e-001 4.9427e-001 4.3100e-001 x = 1.2500e+000 -> 9.9742e-001 8.8879e-001 1.0259e+000 1.1058e+000 x = 1.5000e+000 -> 1.4094e+000 1.3219e+000 1.4970e+000 1.6900e+000 x = 1.7500e+000 -> 1.8018e+000 1.7517e+000 1.8177e+000 1.9377e+000 x = 2.0000e+000 -> 2.1738e+000 2.1721e+000 2.1251e+000 2.1100e+000 x = 2.2500e+000 -> 2.5242e+000 2.5771e+000 2.5650e+000 2.5099e+000 x = 2.5000e+000 -> 2.8520e+000 2.9611e+000 3.0889e+000 3.1000e+000 x = 2.7500e+000 -> 3.1564e+000 3.3185e+000 3.6001e+000 3.7417e+000 x = 3.0000e+000 -> 3.4361e+000 3.6438e+000 4.0015e+000 4.2300e+000 x = 3.5000e+000 -> 3.9181e+000 4.1757e+000 4.3254e+000 4.3771e+000 x = 4.0000e+000 -> 4.2898e+000 4.5124e+000 4.4740e+000 4.3500e+000 x = 4.2500e+000 -> 4.4333e+000 4.5983e+000 4.6441e+000 4.6021e+000 x = 4.5000e+000 -> 4.5548e+000 4.6462e+000 4.7629e+000 4.8100e+000 x = 4.7500e+000 -> 4.6621e+000 4.6782e+000 4.7182e+000 4.7474e+000 x = 5.0000e+000 -> 4.7627e+000 4.7162e+000 4.6315e+000 4.6100e+000 x = 5.2500e+000 -> 4.8642e+000 4.7795e+000 4.6477e+000 4.6331e+000 x = 5.5000e+000 -> 4.9744e+000 4.8773e+000 4.7718e+000 4.7900e+000 x = 5.7500e+000 -> 5.1008e+000 5.0158e+000 4.9740e+000 5.0000e+000 x = 6.0000e+000 -> 5.2510e+000 5.2015e+000 5.2244e+000 5.2300e+000 x = 6.5000e+000 -> 5.6536e+000 5.7220e+000 5.7912e+000 5.7183e+000 x = 7.0000e+000 -> 6.2432e+000 6.3997e+000 6.4375e+000 6.3500e+000 x = 7.2500e+000 -> 6.6273e+000 6.7796e+000 6.7924e+000 6.7555e+000 x = 7.5000e+000 -> 7.0809e+000 7.1774e+000 7.1695e+000 7.1900e+000 x = 7.7500e+000 -> 7.6119e+000 7.5859e+000 7.5697e+000 7.6095e+000 x = 8.0000e+000 -> 8.2277e+000 7.9979e+000 7.9935e+000 7.9700e+000 End of tests for CURVE_SPLINE_APPROXIMATOR class. Making object KNOTTED_CURVE_SPLINE_APPROXIMATOR - NAG e02bac Knots input = BASIC_VECTOR <1|, real: Index # Values 1 4.0000e+000 Coefficients = BASIC_VECTOR <5|, real: Index # Values 1 -1.2472e+000 1.1886e+000 5.8674e+000 4.4985e+000 8.3139e+000 Knots used = BASIC_VECTOR <9|, real: Index # Values 1 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 4.0000e+000 6 8.0000e+000 8.0000e+000 8.0000e+000 8.0000e+000 Closeness = 0.94936292171380576 Coefficients from CURVE_SPLINE_APPROXIMATOR, same knots = BASIC_VECTOR <5|, real: Index # Values 1 -1.3201e+000 1.3542e+000 5.5510e+000 4.7031e+000 8.2277e+000 Making object KNOTTED_CURVE_SPLINE_APPROXIMATOR - NAG e02bac Knots input = BASIC_VECTOR <5|, real: Index # Values 1 1.0000e+000 2.0000e+000 4.0000e+000 5.0000e+000 6.0000e+000 Coefficients = BASIC_VECTOR <9|, real: Index # Values 1 -1.0921e+000 -7.6519e-001 6.1154e-001 2.6027e+000 5.1276e+000 6 4.4505e+000 5.3499e+000 6.8745e+000 7.9863e+000 Knots used = BASIC_VECTOR <13|, real: Index # Values 1 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 1.0000e+000 6 2.0000e+000 4.0000e+000 5.0000e+000 6.0000e+000 8.0000e+000 11 8.0000e+000 8.0000e+000 8.0000e+000 Closeness = 0.40430594730850983 Coefficients from CURVE_SPLINE_APPROXIMATOR, same knots = BASIC_VECTOR <9|, real: Index # Values 1 -1.1072e+000 -6.5715e-001 4.3500e-001 2.8061e+000 4.6824e+000 6 4.6416e+000 5.1976e+000 6.9008e+000 7.9979e+000 Testing SURFACE_SPLINE_INTERPOLATOR (e01dac) Knots_x = BASIC_VECTOR <11|, real: Index # Values 1 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 1.3000e+000 6 1.5000e+000 1.6000e+000 2.0000e+000 2.0000e+000 2.0000e+000 11 2.0000e+000 Knots_y = BASIC_VECTOR <10|, real: Index # Values 1 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 4.0000e-001 6 7.0000e-001 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 B-spline coefficients = BASIC_MATRIX [7 x 6], real: Row/Col 1 2 3 4 5 1 1.0000e+000 1.1333e+000 1.3667e+000 1.7000e+000 1.9000e+000 2 1.2000e+000 1.3333e+000 1.5667e+000 1.9000e+000 2.1000e+000 3 1.5833e+000 1.7167e+000 1.9500e+000 2.2833e+000 2.4833e+000 4 2.1433e+000 2.2767e+000 2.5100e+000 2.8433e+000 3.0433e+000 5 2.8667e+000 3.0000e+000 3.2333e+000 3.5667e+000 3.7667e+000 6 3.4667e+000 3.6000e+000 3.8333e+000 4.1667e+000 4.3667e+000 7 4.0000e+000 4.1333e+000 4.3667e+000 4.7000e+000 4.9000e+000 Row/Col 6 1 2.0000e+000 2 2.2000e+000 3 2.5833e+000 4 3.1433e+000 5 3.8667e+000 6 4.4667e+000 7 5.0000e+000 Computing 'values_on_grid' on test grid values_on_grid = BASIC_MATRIX [6 x 6], real: Row/Col 1 2 3 4 5 1 1.0000e+000 1.2000e+000 1.4000e+000 1.6000e+000 1.8000e+000 2 1.4400e+000 1.6400e+000 1.8400e+000 2.0400e+000 2.2400e+000 3 1.9600e+000 2.1600e+000 2.3600e+000 2.5600e+000 2.7600e+000 4 2.5600e+000 2.7600e+000 2.9600e+000 3.1600e+000 3.3600e+000 5 3.2400e+000 3.4400e+000 3.6400e+000 3.8400e+000 4.0400e+000 6 4.0000e+000 4.2000e+000 4.4000e+000 4.6000e+000 4.8000e+000 Row/Col 6 1 2.0000e+000 2 2.4400e+000 3 2.9600e+000 4 3.5600e+000 5 4.2400e+000 6 5.0000e+000 x1 = BASIC_VECTOR <6|, real: Index # Values 1 1.0000e+000 1.2000e+000 1.4000e+000 1.6000e+000 1.8000e+000 6 2.0000e+000 y1 = BASIC_VECTOR <6|, real: Index # Values 1 0.0000e+000 2.0000e-001 4.0000e-001 6.0000e-001 8.0000e-001 6 1.0000e+000 values_at_points(x1, y1) = BASIC_VECTOR <6|, real: Index # Values 1 1.0000e+000 1.6400e+000 2.3600e+000 3.1600e+000 4.0400e+000 6 5.0000e+000 End of tests for SURFACE_SPLINE_INTERPOLATOR class Test of SURFACE_SPLINE_APPROXIMATOR Making object SURFACE_SPLINE_APPROXIMATOR - NAG e02dcc with assumed smoothness = 0.1 Closeness (sum of squared residuals) = = 0.10004011185035931 Distinct knots_x - see e02dcc test results 0.000e+000 1.500e+000 2.500e+000 5.000e+000 Distinct knots_y - see e02dcc test results 0.000e+000 1.000e+000 2.000e+000 2.500e+000 3.000e+000 3.500e+000 4.000e+000 B-spline coefficients - see e02dcc test results = BASIC_MATRIX [6 x 9], real: Row/Col 1 2 3 4 5 1 9.9180e-001 1.0546e+000 6.0982e-001 -2.9148e-001 -8.4761e-001 2 1.5381e+000 1.5270e+000 9.5574e-001 -4.1989e-001 -1.3296e+000 3 2.3913e+000 2.2441e+000 1.5587e+000 -7.3994e-001 -1.8521e+000 4 3.9845e+000 4.2217e+000 2.3458e+000 -1.1763e+000 -3.3468e+000 5 5.2138e+000 5.0860e+000 3.3860e+000 -1.5527e+000 -4.3628e+000 6 5.9965e+000 6.0821e+000 3.7716e+000 -1.7775e+000 -5.0085e+000 Row/Col 6 7 8 9 1 -1.0168e+000 -9.5292e-001 -7.7107e-001 -6.4757e-001 2 -1.5952e+000 -1.3381e+000 -1.0914e+000 -1.0373e+000 3 -2.4022e+000 -2.2844e+000 -1.8488e+000 -1.5936e+000 4 -3.9390e+000 -3.9559e+000 -3.2549e+000 -2.5887e+000 5 -5.4680e+000 -5.0032e+000 -3.9444e+000 -3.3485e+000 6 -6.1656e+000 -5.8709e+000 -4.7297e+000 -3.9330e+000 Values on grid (e01dcc) = BASIC_MATRIX [6 x 5], real: Row/Col 1 2 3 4 5 1 9.9180e-001 5.4074e-001 -4.1677e-001 -9.7795e-001 -6.4757e-001 2 2.0427e+000 1.0881e+000 -8.2926e-001 -1.9748e+000 -1.3627e+000 3 3.0286e+000 1.6071e+000 -1.2410e+000 -2.9145e+000 -1.9911e+000 4 4.0135e+000 2.1422e+000 -1.6647e+000 -3.9133e+000 -2.6055e+000 5 5.0213e+000 2.7054e+000 -2.0831e+000 -4.9652e+000 -3.2510e+000 6 5.9965e+000 3.2394e+000 -2.4846e+000 -5.9236e+000 -3.9330e+000 End of tests for SURFACE_SPLINE_APPROXIMATOR class e01sac Example Program Data x = BASIC_VECTOR <30|, real: Index # Values 1 1.1160e+001 1.2850e+001 1.9850e+001 1.9720e+001 1.5910e+001 6 0.0000e+000 2.0870e+001 3.4500e+000 1.4260e+001 1.7430e+001 11 2.2800e+001 7.5800e+000 2.5000e+001 0.0000e+000 9.6600e+000 16 5.2200e+000 1.7250e+001 2.5000e+001 1.2130e+001 2.2230e+001 21 1.1520e+001 1.5200e+001 7.5400e+000 1.7320e+001 2.1400e+000 26 5.1000e-001 2.2690e+001 5.4700e+000 2.1670e+001 3.3100e+000 y = BASIC_VECTOR <30|, real: Index # Values 1 1.2400e+000 3.0600e+000 1.0720e+001 1.3900e+000 7.7400e+000 6 2.0000e+001 2.0000e+001 1.2780e+001 1.7870e+001 3.4600e+000 11 1.2390e+001 1.9800e+000 1.1870e+001 0.0000e+000 2.0000e+001 16 1.4660e+001 1.9570e+001 3.8700e+000 1.0790e+001 6.2100e+000 21 8.5300e+000 0.0000e+000 1.0690e+001 1.3780e+001 1.5030e+001 26 8.3700e+000 1.9630e+001 1.7130e+001 1.4360e+001 3.3000e-001 f = BASIC_VECTOR <30|, real: Index # Values 1 2.2150e+001 2.2110e+001 7.9700e+000 1.6830e+001 1.5300e+001 6 3.4600e+001 5.7400e+000 4.1240e+001 1.0740e+001 1.8600e+001 11 5.4700e+000 2.9870e+001 4.4000e+000 5.8200e+001 4.7300e+000 16 4.0360e+001 6.4300e+000 8.7400e+000 1.3710e+001 1.0250e+001 21 1.5740e+001 2.1600e+001 1.9310e+001 1.2110e+001 5.3100e+001 26 4.9430e+001 3.2500e+000 2.8630e+001 5.5200e+000 4.4080e+001 Renka-Cline method values_on_grid = BASIC_MATRIX [7 x 6], real: Row/Col 1 2 3 4 5 1 4.3519e+001 4.0485e+001 3.7900e+001 3.8552e+001 4.7609e+001 2 3.3905e+001 2.9256e+001 2.3966e+001 2.5253e+001 3.6658e+001 3 2.6589e+001 2.2511e+001 1.6789e+001 1.6724e+001 2.2869e+001 4 2.2226e+001 2.0715e+001 1.6427e+001 1.3831e+001 1.4018e+001 5 2.1154e+001 1.9299e+001 1.5461e+001 1.3079e+001 1.3443e+001 6 1.8670e+001 1.6724e+001 1.3023e+001 1.0715e+001 1.1195e+001 7 1.4880e+001 1.2869e+001 9.2982e+000 6.8805e+000 6.4607e+000 Row/Col 6 1 4.1245e+001 2 2.7620e+001 3 1.8035e+001 4 1.2287e+001 5 1.1684e+001 6 9.0934e+000 7 5.3679e+000 rnw_used = 9.487744680835549 rnq_used = 13.417697203970826 minnq = 7 values_on_grid = BASIC_MATRIX [7 x 6], real: Row/Col 1 2 3 4 5 1 4.4583e+001 4.1251e+001 3.7502e+001 3.9419e+001 4.6958e+001 2 3.4353e+001 3.1761e+001 2.2360e+001 2.5417e+001 3.7369e+001 3 2.6470e+001 2.4736e+001 1.8569e+001 1.6318e+001 2.3742e+001 4 2.2272e+001 2.1169e+001 1.5629e+001 1.3776e+001 1.4670e+001 5 2.0981e+001 1.8934e+001 1.5548e+001 1.2602e+001 1.3245e+001 6 1.8690e+001 1.6830e+001 1.3049e+001 1.0391e+001 1.1293e+001 7 1.5057e+001 1.2651e+001 9.6903e+000 7.0334e+000 6.2588e+000 Row/Col 6 1 4.0233e+001 2 2.7720e+001 3 2.1232e+001 4 1.4595e+001 5 1.1998e+001 6 9.4282e+000 7 5.4643e+000 End of tests for SCATTERED_SURFACE_INTERPOLATOR class e02ddc Example Program Data x = BASIC_VECTOR <30|, real: Index # Values 1 1.1160e+001 1.2850e+001 1.9850e+001 1.9720e+001 1.5910e+001 6 0.0000e+000 2.0870e+001 3.4500e+000 1.4260e+001 1.7430e+001 11 2.2800e+001 7.5800e+000 2.5000e+001 0.0000e+000 9.6600e+000 16 5.2200e+000 1.7250e+001 2.5000e+001 1.2130e+001 2.2230e+001 21 1.1520e+001 1.5200e+001 7.5400e+000 1.7320e+001 2.1400e+000 26 5.1000e-001 2.2690e+001 5.4700e+000 2.1670e+001 3.3100e+000 y = BASIC_VECTOR <30|, real: Index # Values 1 1.2400e+000 3.0600e+000 1.0720e+001 1.3900e+000 7.7400e+000 6 2.0000e+001 2.0000e+001 1.2780e+001 1.7870e+001 3.4600e+000 11 1.2390e+001 1.9800e+000 1.1870e+001 0.0000e+000 2.0000e+001 16 1.4660e+001 1.9570e+001 3.8700e+000 1.0790e+001 6.2100e+000 21 8.5300e+000 0.0000e+000 1.0690e+001 1.3780e+001 1.5030e+001 26 8.3700e+000 1.9630e+001 1.7130e+001 1.4360e+001 3.3000e-001 f = BASIC_VECTOR <30|, real: Index # Values 1 2.2150e+001 2.2110e+001 7.9700e+000 1.6830e+001 1.5300e+001 6 3.4600e+001 5.7400e+000 4.1240e+001 1.0740e+001 1.8600e+001 11 5.4700e+000 2.9870e+001 4.4000e+000 5.8200e+001 4.7300e+000 16 4.0360e+001 6.4300e+000 8.7400e+000 1.3710e+001 1.0250e+001 21 1.5740e+001 2.1600e+001 1.9310e+001 1.2110e+001 5.3100e+001 26 4.9430e+001 3.2500e+000 2.8630e+001 5.5200e+000 4.4080e+001 w = BASIC_VECTOR <30|, real: Index # Values 1 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 6 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 11 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 16 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 21 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 26 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 1.0000e+000 s = 10 rank = 30 closeness = 10.002108871019857 values_on_grid = BASIC_MATRIX [7 x 6], real: Row/Col 1 2 3 4 5 1 4.5039e+001 3.6639e+001 3.0251e+001 3.7265e+001 4.8338e+001 2 3.3696e+001 2.6748e+001 1.9659e+001 2.4456e+001 3.3966e+001 3 2.6247e+001 2.3069e+001 1.6904e+001 1.7209e+001 2.1556e+001 4 2.2876e+001 2.1134e+001 1.6285e+001 1.4137e+001 1.4711e+001 5 2.1618e+001 1.8970e+001 1.5211e+001 1.3025e+001 1.2318e+001 6 1.9395e+001 1.5897e+001 1.2708e+001 1.1234e+001 1.0816e+001 7 1.3405e+001 1.1982e+001 8.9851e+000 7.2938e+000 7.1497e+000 Row/Col 6 1 4.0743e+001 2 2.8624e+001 3 1.9842e+001 4 1.4294e+001 5 1.1208e+001 6 9.4621e+000 7 7.0894e+000 rank = 30 closeness = 9.0000426592967191 values_on_grid = BASIC_MATRIX [7 x 6], real: Row/Col 1 2 3 4 5 1 4.5770e+001 3.7285e+001 3.0152e+001 3.6985e+001 4.8373e+001 2 3.3916e+001 2.6942e+001 1.9544e+001 2.4346e+001 3.4069e+001 3 2.6246e+001 2.3113e+001 1.6898e+001 1.7167e+001 2.1508e+001 4 2.2834e+001 2.1145e+001 1.6332e+001 1.4118e+001 1.4595e+001 5 2.1617e+001 1.8979e+001 1.5232e+001 1.3012e+001 1.2259e+001 6 1.9419e+001 1.5897e+001 1.2674e+001 1.1215e+001 1.0844e+001 7 1.3347e+001 1.1957e+001 8.9442e+000 7.2767e+000 7.1913e+000 Row/Col 6 1 4.0683e+001 2 2.8581e+001 3 1.9770e+001 4 1.4221e+001 5 1.1176e+001 6 9.5019e+000 7 7.1632e+000 End of tests for SCATTERED_SURFACE_SPLINE_APPROXIMATOR class