Testing was performed by comparing the results given by the dasl and the optimized algorithm C. The paragraphs to follow show the results of two different input files, one with `large' rms deformation ( p6.spl ) and one with `small' rms deformation ( p6_low_31_72.SPL ). The results agree to within roundoff errors, this was done by running the following command:
The file ./test/p6.spl with the following characteristics: # spline deformation amplitudes # z_order 4 # theta_order 4 # qz 421 # num_z_knots 417 # num_theta_knots 144 # zmin zmax -1.000000e+00 1.000000e+00 # tmin tmax 0.000000e+00 6.283185e+00 # rms amplitude 129188.64688685712463 # yield the following results: futile-90: objs/sun4u-SunOS-5/evalbspline num=100 test=5 \ | column -v -a deltas N - | column -v -a deltadt N \ | column -v -a deltadz N \ | compute deltas = \( dasl - algoc \) \/ dasl \ | compute deltadz = \( ddasldz - dalgocdz \) \/ ddasldz \ | compute deltadt = \( ddasldt - dalgocdt \) \/ ddasldt \ | rdbstats deltas deltadz deltadt | tbl2lst deltas_n | N deltas_sum | N deltas_ave | N deltas_dev | N deltas_min | N deltas_max | N deltadz_n | N deltadz_sum | N deltadz_ave | N deltadz_dev | N deltadz_min | N deltadz_max | N deltadt_n | N deltadt_sum | N deltadt_ave | N deltadt_dev | N deltadt_min | N deltadt_max | N deltas_n | 10000 deltas_sum | -1.29628283548978e-13 deltas_ave | -1.29628283548978e-17 deltas_dev | 1.10245039134638e-16 deltas_min | -5.83761910253855e-16 deltas_max | 5.3919301269975e-16 deltadz_n | 10000 deltadz_sum | -9.28799965664604e-14 deltadz_ave | -9.28799965664602e-18 deltadz_dev | 1.33699222283282e-16 deltadz_min | -3.90195209692813e-15 deltadz_max | 3.74917737998133e-15 deltadt_n | 10000 deltadt_sum | -1.46179231627943e-13 deltadt_ave | -1.46179231627943e-17 deltadt_dev | 1.07162272062163e-16 deltadt_min | -6.82814047720486e-16 deltadt_max | 7.42207044201643e-16 The file ./test/p6_low_31_72.SPL with the following characteristics: # spline deformation amplitudes # z_order 4 # theta_order 4 # qz 33 # num_z_knots 29 # num_theta_knots 71 # zmin zmax -1.000000e+00 1.000000e+00 # tmin tmax 0.000000e+00 6.283185e+00 # rms amplitude 0.00101284258284 yield the following results: futile-93: objs/sun4u-SunOS-5/evalbspline num=100 test=5 \ input=p6_low_31_72.SPL \ | column -v -a deltas N - | column -v -a deltadt N \ | column -v -a deltadz N \ | compute deltas = \( dasl - algoc \) \/ dasl \ | compute deltadz = \( ddasldz - dalgocdz \) \/ ddasldz \ | compute deltadt = \( ddasldt - dalgocdt \) \/ ddasldt \ | rdbstats deltas deltadz deltadt | tbl2lst deltas_n | N deltas_sum | N deltas_ave | N deltas_dev | N deltas_min | N deltas_max | N deltadz_n | N deltadz_sum | N deltadz_ave | N deltadz_dev | N deltadz_min | N deltadz_max | N deltadt_n | N deltadt_sum | N deltadt_ave | N deltadt_dev | N deltadt_min | N deltadt_max | N deltas_n | 10000 deltas_sum | -0 deltas_ave | 0 deltas_dev | 0 deltas_min | -0 deltas_max | -0 deltadz_n | 10000 deltadz_sum | 0 deltadz_ave | 0 deltadz_dev | 0 deltadz_min | 0 deltadz_max | 0 deltadt_n | 10000 deltadt_sum | 0 deltadt_ave | 0 deltadt_dev | 0 deltadt_min | -0 deltadt_max | -0