04 March 2022 02:13:07 PM FN_TEST: C++ version Test the FN library. I4_MACH_TEST: I4_MACH evaluates integer machine numbers. I4_MACH(1) = the standard input unit. I4_MACH(2) = the standard output unit. I4_MACH(3) = the standard punch unit. I4_MACH(4) = the standard error message unit. I4_MACH(5) = the number of bits per integer storage unit. I4_MACH(6) = the number of characters per integer storage unit. I4_MACH(7) = A, the base. I4_MACH(8) = S, the number of base A digits. I4_MACH(9) = A^S-1, the largest integer. I4_MACH(10) = B, the base. I4_MACH(11) = T, the number of single precision base B digits. I4_MACH(12) = EMIN, the smallest single precision exponent E. I4_MACH(13) = EMAX, the largest single precision exponent E. I4_MACH(14) = T, the number of double precision base B digits. I4_MACH(15) = EMIN, the smallest double precision exponent E. I4_MACH(16) = EMAX, the largest double precision exponent E. I I4_MACH(I) 1 5 2 6 3 7 4 6 5 32 6 4 7 2 8 31 9 2147483647 10 2 11 24 12 -125 13 128 14 53 15 -1021 16 1024 R8_ACOS_TEST: R8_ACOS evaluates the arccosine function. X ARCCOS(X) R8_ACOS(X) Diff -0.1 1.67096 1.67096 2.22045e-16 0 1.5708 1.5708 0 0.1 1.47063 1.47063 0 0.2 1.36944 1.36944 2.22045e-16 0.3 1.2661 1.2661 2.22045e-16 0.4 1.15928 1.15928 0 0.5 1.0472 1.0472 0 0.6 0.927295 0.927295 0 0.7 0.795399 0.795399 0 0.8 0.643501 0.643501 1.11022e-16 0.9 0.451027 0.451027 5.55112e-17 1 0 0 0 R8_ACOSH_TEST: R8_ACOSH evaluates the hyperbolic arccosine function X ARCCOSH(X) R8_ACOSH(X) Diff 1 0 0 0 1.01 0.141304 0.141304 8.32667e-17 1.1 0.443568 0.443568 1.11022e-16 1.2 0.622363 0.622363 1.11022e-16 1.3 0.756433 0.756433 0 1.4 0.867015 0.867015 1.11022e-16 1.5 0.962424 0.962424 0 2 1.31696 1.31696 2.22045e-16 3 1.76275 1.76275 0 3.14159 1.81153 1.81153 0 4 2.06344 2.06344 4.44089e-16 5 2.29243 2.29243 0 10 2.99322 2.99322 0 100 5.29829 5.29829 0 1000 7.6009 7.6009 0 R8_AI_TEST: Test R8_AI. X AIRY_AI(X) R8_AI(X) Diff 0 0.355028 0.355028 0 0.1 0.329203 0.329203 0 0.2 0.303703 0.303703 0 0.3 0.278806 0.278806 0 0.4 0.254742 0.254742 0 0.5 0.231694 0.231694 0 0.6 0.2098 0.2098 2.77556e-17 0.7 0.189162 0.189162 5.55112e-17 0.8 0.169846 0.169846 5.55112e-17 0.9 0.151887 0.151887 5.55112e-17 1 0.135292 0.135292 2.77556e-17 R8_AID_TEST: Test R8_AID. X AIRY_AID(X) R8_AID(X) Diff 0 -0.258819 -0.258819 0 0.1 -0.25713 -0.25713 5.55112e-17 0.2 -0.252405 -0.252405 5.55112e-17 0.3 -0.245146 -0.245146 0 0.4 -0.235832 -0.235832 2.77556e-17 0.5 -0.224911 -0.224911 0 0.6 -0.212793 -0.212793 2.77556e-17 0.7 -0.199851 -0.199851 2.77556e-17 0.8 -0.186413 -0.186413 0 0.9 -0.172764 -0.172764 0 1 -0.159147 -0.159147 0 R8_AINT_TEST: R8_AINT rounds an R8 towards 0. X AINT(X) R8_AINT(X) Diff -2.01 -2 -2 0 -1.99 -1 -1 0 -1.5 -1 -1 0 -1.1 -1 -1 0 -1.01 -1 -1 0 -1 -1 -1 0 -0.99 0 -0 0 -0.9 0 -0 0 -0.51 0 -0 0 -0.5 0 -0 0 -0.49 0 -0 0 -0.01 0 -0 0 0 0 0 0 0.01 0 0 0 0.49 0 0 0 0.5 0 0 0 0.51 0 0 0 0.9 0 0 0 0.99 0 0 0 1 1 1 0 1.01 1 1 0 1.1 1 1 0 1.5 1 1 0 1.99 1 1 0 2.01 2 2 0 R8_ASIN_TEST: Test R8_ASIN. X ARCSIN(X) R8_ASIN(X) Diff -0.1 -0.100167 -0.100167 0 0 0 0 0 0.1 0.100167 0.100167 0 0.2 0.201358 0.201358 0 0.3 0.304693 0.304693 0 0.4 0.411517 0.411517 5.55112e-17 0.5 0.523599 0.523599 1.11022e-16 0.6 0.643501 0.643501 0 0.7 0.775397 0.775397 1.11022e-16 0.8 0.927295 0.927295 1.11022e-16 0.9 1.11977 1.11977 0 1 1.5708 1.5708 0 R8_ASINH_TEST: Test R8_ASINH X ARCSINH(X) R8_ASINH(X) Diff -5 -2.31244 -2.31244 0 -1 -0.881374 -0.881374 0 0 0 0 0 0.1 0.0998341 0.0998341 1.38778e-17 0.2 0.19869 0.19869 0 0.3 0.295673 0.295673 0 0.4 0.390035 0.390035 5.55112e-17 0.5 0.481212 0.481212 0 0.6 0.568825 0.568825 0 0.7 0.652667 0.652667 0 0.8 0.732668 0.732668 1.11022e-16 0.9 0.808867 0.808867 0 1 0.881374 0.881374 0 2 1.44364 1.44364 0 3 1.81845 1.81845 0 4 2.09471 2.09471 0 5 2.31244 2.31244 0 10 2.99822 2.99822 0 100 5.29834 5.29834 0 1000 7.6009 7.6009 0 R8_ATAN_TEST: Test R8_ATAN. X ARCTAN(X) R8_ATAN(X) Diff 0 0 0 0 0.25 0.244979 0.244979 2.77556e-17 0.333333 0.321751 0.321751 0 0.5 0.463648 0.463648 5.55112e-17 1 0.785398 0.785398 0 2 1.10715 1.10715 2.22045e-16 3 1.24905 1.24905 0 4 1.32582 1.32582 0 5 1.3734 1.3734 0 10 1.47113 1.47113 0 20 1.52084 1.52084 0 R8_ATAN2_TEST: Test R8_ATAN2. X Y ARCTAN2(Y,X) R8_ATAN2(Y,X) Diff 0 -1 -1.5708 -1.5708 0 0.5 -0.866025 -1.0472 -1.0472 2.22045e-16 0.866025 -0.5 -0.523599 -0.523599 0 1 0 0 0 0 0.866025 0.5 0.523599 0.523599 0 0.5 0.866025 1.0472 1.0472 2.22045e-16 0 1 1.5708 1.5708 0 -0.5 0.866025 2.0944 2.0944 0 -0.866025 0.5 2.61799 2.61799 4.44089e-16 -1 0 3.14159 3.14159 0 -0.866025 -0.5 -2.61799 -2.61799 4.44089e-16 -0.5 -0.866025 -2.0944 -2.0944 0 0 -1 -1.5708 -1.5708 0 0.5 -0.866025 -1.0472 -1.0472 2.22045e-16 0.866025 -0.5 -0.523599 -0.523599 0 1 0 0 0 0 0.866025 0.5 0.523599 0.523599 0 0.5 0.866025 1.0472 1.0472 2.22045e-16 0 1 1.5708 1.5708 0 R8_ATANH_TEST: Test ARCTANH_VALUES, R4_ATANH, R8_ATANH X ARCTANH(X) R8_ATANH(X) Diff -0.5 -0.549306 -0.549306 0 0 0 0 0 0.001 0.001 0.001 0 0.1 0.100335 0.100335 0 0.2 0.202733 0.202733 2.77556e-17 0.3 0.30952 0.30952 0 0.4 0.423649 0.423649 5.55112e-17 0.5 0.549306 0.549306 0 0.6 0.693147 0.693147 0 0.7 0.867301 0.867301 0 0.8 1.09861 1.09861 0 0.9 1.47222 1.47222 2.22045e-16 0.99 2.64665 2.64665 4.44089e-16 0.999 3.8002 3.8002 0 0.999999 7.25433 7.25433 1.43778e-11 R8_BESI0_TEST: Test R8_BESI0 X BESI0(X) R8_BESI0(X) Diff 0 1 1 0 0.2 1.01003 1.01003 0 0.4 1.0404 1.0404 0 0.6 1.09205 1.09205 4.44089e-16 0.8 1.16651 1.16651 4.44089e-16 1 1.26607 1.26607 6.66134e-16 1.2 1.39373 1.39373 2.22045e-16 1.4 1.5534 1.5534 4.44089e-16 1.6 1.74998 1.74998 4.44089e-16 1.8 1.98956 1.98956 0 2 2.27959 2.27959 4.44089e-16 2.5 3.28984 3.28984 4.44089e-16 3 4.88079 4.88079 0 3.5 7.3782 7.3782 8.88178e-16 4 11.3019 11.3019 1.77636e-15 4.5 17.4812 17.4812 3.55271e-15 5 27.2399 27.2399 0 6 67.2344 67.2344 1.42109e-14 8 427.564 427.564 5.68434e-14 10 2815.72 2815.72 4.54747e-13 R8_BESI1_TEST: Test R8_BESI1 X BESI1(X) R8_BESI1(X) Diff 0 0 0 0 0.2 0.100501 0.100501 4.16334e-17 0.4 0.204027 0.204027 2.77556e-17 0.6 0.313704 0.313704 0 0.8 0.432865 0.432865 5.55112e-17 1 0.565159 0.565159 0 1.2 0.714678 0.714678 0 1.4 0.886092 0.886092 2.22045e-16 1.6 1.08481 1.08481 4.44089e-16 1.8 1.31717 1.31717 0 2 1.59064 1.59064 0 2.5 2.51672 2.51672 4.44089e-16 3 3.95337 3.95337 4.44089e-16 3.5 6.20583 6.20583 0 4 9.75947 9.75947 0 4.5 15.3892 15.3892 3.55271e-15 5 24.3356 24.3356 3.55271e-15 6 61.3419 61.3419 7.10543e-15 8 399.873 399.873 0 10 2670.99 2670.99 4.54747e-13 R8_BESJ0_TEST: Test R8_BESJ0 X BESJ0(X) R8_BESJ0(X) Diff -5 -0.177597 -0.177597 1.38778e-16 -4 -0.39715 -0.39715 5.55112e-17 -3 -0.260052 -0.260052 5.55112e-17 -2 0.223891 0.223891 8.32667e-17 -1 0.765198 0.765198 0 0 1 1 0 1 0.765198 0.765198 0 2 0.223891 0.223891 8.32667e-17 3 -0.260052 -0.260052 5.55112e-17 4 -0.39715 -0.39715 5.55112e-17 5 -0.177597 -0.177597 1.38778e-16 6 0.150645 0.150645 5.55112e-17 7 0.300079 0.300079 0 8 0.171651 0.171651 2.77556e-17 9 -0.0903336 -0.0903336 2.91434e-16 10 -0.245936 -0.245936 0 11 -0.17119 -0.17119 5.55112e-17 12 0.0476893 0.0476893 2.70617e-16 13 0.206926 0.206926 8.32667e-17 14 0.171073 0.171073 0 15 -0.0142245 -0.0142245 1.30104e-16 R8_BESJ1_TEST: Test R8_BESJ1 X BESJ1(X) R8_BESJ1(X) Diff -5 0.327579 0.327579 1.11022e-16 -4 0.0660433 0.0660433 1.52656e-16 -3 -0.339059 -0.339059 0 -2 -0.576725 -0.576725 0 -1 -0.440051 -0.440051 0 0 0 0 0 1 0.440051 0.440051 0 2 0.576725 0.576725 0 3 0.339059 0.339059 0 4 -0.0660433 -0.0660433 1.52656e-16 5 -0.327579 -0.327579 1.11022e-16 6 -0.276684 -0.276684 0 7 -0.00468282 -0.00468282 9.28077e-17 8 0.234636 0.234636 0 9 0.245312 0.245312 2.77556e-17 10 0.0434727 0.0434727 5.55112e-17 11 -0.176785 -0.176785 2.77556e-16 12 -0.223447 -0.223447 1.11022e-16 13 -0.0703181 -0.0703181 2.08167e-16 14 0.133375 0.133375 0 15 0.205104 0.205104 5.55112e-17 R8_BESK0_TEST: Test R8_BESK0 X BESK0(X) R8_BESK0(X) Diff 0.1 2.42707 2.42707 4.44089e-16 0.2 1.7527 1.7527 0 0.4 1.11453 1.11453 0 0.6 0.777522 0.777522 1.11022e-16 0.8 0.565347 0.565347 1.11022e-16 1 0.421024 0.421024 1.66533e-16 1.2 0.318508 0.318508 5.55112e-17 1.4 0.243655 0.243655 0 1.6 0.187955 0.187955 5.55112e-17 1.8 0.145931 0.145931 0 2 0.113894 0.113894 1.38778e-17 2.5 0.0623476 0.0623476 0 3 0.0347395 0.0347395 6.93889e-18 3.5 0.0195989 0.0195989 0 4 0.0111597 0.0111597 3.46945e-18 4.5 0.00639986 0.00639986 0 5 0.0036911 0.0036911 4.33681e-19 6 0.00124399 0.00124399 2.1684e-19 8 0.000146471 0.000146471 0 10 1.77801e-05 1.77801e-05 3.38813e-21 R8_BESK1_TEST: Test R8_BESK1 X BESK1(X) R8_BESK1(X) Diff 0.1 9.85384 9.85384 0 0.2 4.77597 4.77597 8.88178e-16 0.4 2.18435 2.18435 0 0.6 1.30283 1.30283 2.22045e-16 0.8 0.861782 0.861782 0 1 0.601907 0.601907 0 1.2 0.434592 0.434592 0 1.4 0.320836 0.320836 0 1.6 0.240634 0.240634 2.77556e-17 1.8 0.182623 0.182623 1.11022e-16 2 0.139866 0.139866 2.77556e-17 2.5 0.0738908 0.0738908 0 3 0.0401564 0.0401564 6.93889e-18 3.5 0.0222394 0.0222394 3.46945e-18 4 0.0124835 0.0124835 0 4.5 0.00707809 0.00707809 0 5 0.00404461 0.00404461 8.67362e-19 6 0.00134392 0.00134392 2.1684e-19 8 0.000155369 0.000155369 2.71051e-20 10 1.86488e-05 1.86488e-05 3.38813e-21 R8_BESY0_TEST: Test R8_BESY0 X BESY0(X) R8_BESY0(X) Diff 0.1 -1.53424 -1.53424 2.22045e-16 1 0.088257 0.088257 4.16334e-17 2 0.510376 0.510376 0 3 0.37685 0.37685 0 4 -0.0169407 -0.0169407 3.1225e-17 5 -0.308518 -0.308518 5.55112e-17 6 -0.288195 -0.288195 0 7 -0.0259497 -0.0259497 1.11022e-16 8 0.223521 0.223521 5.55112e-17 9 0.249937 0.249937 5.55112e-17 10 0.0556712 0.0556712 8.32667e-17 11 -0.168847 -0.168847 0 12 -0.225237 -0.225237 8.32667e-17 13 -0.0782079 -0.0782079 2.91434e-16 14 0.127193 0.127193 2.77556e-17 15 0.205464 0.205464 2.77556e-17 R8_BESY1_TEST: Test R8_BESY1 X BESY1(X) R8_BESY1(X) Diff 0.1 -6.45895 -6.45895 8.88178e-16 1 -0.781213 -0.781213 1.11022e-16 2 -0.107032 -0.107032 0 3 0.324674 0.324674 0 4 0.397926 0.397926 5.55112e-17 5 0.147863 0.147863 2.77556e-17 6 -0.17501 -0.17501 2.77556e-17 7 -0.302667 -0.302667 1.11022e-16 8 -0.15806 -0.15806 5.55112e-17 9 0.104315 0.104315 1.38778e-17 10 0.249015 0.249015 5.55112e-17 11 0.163706 0.163706 3.05311e-16 12 -0.0570992 -0.0570992 4.09395e-16 13 -0.210081 -0.210081 8.32667e-17 14 -0.166645 -0.166645 5.55112e-17 15 0.0210736 0.0210736 4.85723e-17 R8_BETA_TEST: Test R8_BETA. X BETA(A,B) R8_BETA(A,B) Diff 0.2 1 5 5 8.88178e-16 0.4 1 2.5 2.5 0 0.6 1 1.66667 1.66667 4.44089e-16 0.8 1 1.25 1.25 2.22045e-16 1 0.2 5 5 8.88178e-16 1 0.4 2.5 2.5 0 1 1 1 1 0 2 2 0.166667 0.166667 5.55112e-17 3 3 0.0333333 0.0333333 0 4 4 0.00714286 0.00714286 0 5 5 0.0015873 0.0015873 2.1684e-19 6 2 0.0238095 0.0238095 3.46945e-18 6 3 0.00595238 0.00595238 0 6 4 0.00198413 0.00198413 0 6 5 0.000793651 0.000793651 1.84314e-18 6 6 0.00036075 0.00036075 1.6263e-18 7 7 8.32501e-05 8.32501e-05 1.49078e-19 R8_BETAI_TEST: Test R8_BETAI. X BETA(A,B,X) R8_BETAI(A,B,X) Diff 0.5 0.5 0.01 0.0637686 0.0637686 0 0.5 0.5 0.1 0.204833 0.204833 5.55112e-17 0.5 0.5 1 1 1 0 1 0.5 0 0 0 0 1 0.5 0.01 0.00501256 0.00501256 3.46945e-18 1 0.5 0.1 0.0513167 0.0513167 1.38778e-17 1 0.5 0.5 0.292893 0.292893 1.11022e-16 1 1 0.5 0.5 0.5 0 2 2 0.1 0.028 0.028 1.04083e-17 2 2 0.2 0.104 0.104 1.38778e-17 2 2 0.3 0.216 0.216 8.32667e-17 2 2 0.4 0.352 0.352 5.55112e-17 2 2 0.5 0.5 0.5 0 2 2 0.6 0.648 0.648 1.11022e-16 2 2 0.7 0.784 0.784 1.11022e-16 2 2 0.8 0.896 0.896 0 2 2 0.9 0.972 0.972 0 5.5 5 0.5 0.436191 0.436191 7.77156e-16 10 0.5 0.9 0.151641 0.151641 1.44329e-15 10 5 0.5 0.0897827 0.0897827 2.63678e-16 10 5 1 1 1 0 10 10 0.5 0.5 0.5 1.11022e-16 20 5 0.8 0.459877 0.459877 1.11022e-16 20 10 0.6 0.214682 0.214682 1.08247e-15 20 10 0.8 0.950736 0.950736 2.22045e-16 20 20 0.5 0.5 0.5 0 20 20 0.6 0.897941 0.897941 2.55351e-15 30 10 0.7 0.22413 0.22413 4.16334e-16 30 10 0.8 0.758641 0.758641 3.33067e-16 40 20 0.7 0.700178 0.700178 5.55112e-16 1 0.5 0.1 0.0513167 0.0513167 1.38778e-17 1 0.5 0.2 0.105573 0.105573 1.38778e-17 1 0.5 0.3 0.16334 0.16334 0 1 0.5 0.4 0.225403 0.225403 8.32667e-17 1 2 0.2 0.36 0.36 0 1 3 0.2 0.488 0.488 1.11022e-16 1 4 0.2 0.5904 0.5904 1.11022e-16 1 5 0.2 0.67232 0.67232 1.11022e-16 2 2 0.3 0.216 0.216 8.32667e-17 3 2 0.3 0.0837 0.0837 5.55112e-17 4 2 0.3 0.03078 0.03078 1.73472e-17 5 2 0.3 0.010935 0.010935 3.46945e-18 1.30625 11.7562 0.225609 0.918885 0.918885 2.22045e-16 1.30625 11.7562 0.0335568 0.21053 0.21053 3.88578e-16 1.30625 11.7562 0.0295222 0.182413 0.182413 2.498e-16 R8_BI_TEST: Test R8_BI. X AIRY_BI(X) R8_BI(X) Diff 0 0.614927 0.614927 0 0.1 0.659862 0.659862 1.11022e-16 0.2 0.705464 0.705464 0 0.3 0.752486 0.752486 1.11022e-16 0.4 0.801773 0.801773 1.11022e-16 0.5 0.854277 0.854277 0 0.6 0.911063 0.911063 1.11022e-16 0.7 0.973329 0.973329 1.11022e-16 0.8 1.04242 1.04242 2.22045e-16 0.9 1.11987 1.11987 2.22045e-16 1 1.20742 1.20742 2.22045e-16 R8_BID_TEST: Test R8_BID. X AIRY_BID(X) R8_BID(X) Diff 0 0.448288 0.448288 5.55112e-17 0.1 0.451513 0.451513 0 0.2 0.461789 0.461789 5.55112e-17 0.3 0.480049 0.480049 5.55112e-17 0.4 0.507282 0.507282 1.11022e-16 0.5 0.544573 0.544573 0 0.6 0.593144 0.593144 0 0.7 0.654406 0.654406 0 0.8 0.730007 0.730007 1.11022e-16 0.9 0.821904 0.821904 1.11022e-16 1 0.932436 0.932436 0 R8_BINOM_TEST: Test R8_BINOM. X BINOM(A,B) R8_BINOM(A,B) Diff 1 0 1 1 0 6 1 6 6 0 6 3 20 20 0 6 5 6 6 0 15 1 15 15 0 15 3 455 455 0 15 5 3003 3003 0 15 7 6435 6435 0 15 9 5005 5005 0 15 11 1365 1365 0 15 13 105 105 0 25 1 25 25 0 25 3 2300 2300 0 25 5 53130 53130 0 25 7 480700 480700 0 25 9 2042975 2.04298e+06 0 25 11 4457400 4.4574e+06 0 25 13 5200300 5.2003e+06 0 25 15 3268760 3.26876e+06 0 25 17 1081575 1.08158e+06 0 R8_CBRT_TEST: Test R8_CBRT X CBRT(X) R8_CBRT(X) Diff 0 0 0 0 -8.1e-09 -0.0020083 -0.0020083 4.33681e-19 0.09 0.44814 0.44814 0 -0.1 -0.464159 -0.464159 5.55112e-17 0.4 0.736806 0.736806 1.11022e-16 -1 -1 -1 0 2 1.25992 1.25992 0 -3 -1.44225 -1.44225 2.22045e-16 3.14159 1.46459 1.46459 0 -19 -2.6684 -2.6684 4.44089e-16 29 3.07232 3.07232 4.44089e-16 -71 -4.14082 -4.14082 0 97 4.5947 4.5947 8.88178e-16 -1.23457e+08 -497.934 -497.934 0 R8_CHI_TEST: Test R8_CHI. X CHI(X) R8_CHI(X) Diff 0.5 -0.0527768 -0.0527768 6.93889e-18 0.6 0.157751 0.157751 8.32667e-17 0.7 0.345569 0.345569 1.11022e-16 0.8 0.5184 0.5184 1.11022e-16 0.9 0.681314 0.681314 1.11022e-16 1 0.837867 0.837867 1.11022e-16 1.2 1.14184 1.14184 6.66134e-16 1.4 1.44549 1.44549 2.22045e-16 1.6 1.75951 1.75951 0 1.8 2.09258 2.09258 8.88178e-16 2 2.45267 2.45267 4.44089e-16 2.5 3.52443 3.52443 4.44089e-16 3 4.96039 4.96039 8.88178e-16 3.5 6.95919 6.95919 0 4 9.81355 9.81355 3.55271e-15 4.5 13.9658 13.9658 3.55271e-15 R8_CHU_TEST: Test R8_CHU. A B X CHU(A,B,X) R8_CHU(A,B,X) Diff -2.5 3.3 0.25 -68.6936 -68.6936 5.68434e-14 -0.5 1.1 0.25 -0.00297106 -0.00297106 4.75748e-16 0.5 1.1 0.25 1.50086 1.50086 1.11022e-15 2.5 3.3 0.25 20.6147 20.6147 1.06581e-14 -2.5 3.3 1.55 7.45638 7.45638 1.77636e-15 -0.5 1.1 1.55 1.01558 1.01558 2.22045e-16 0.5 1.1 1.55 0.734465 0.734465 9.99201e-16 2.5 3.3 1.55 0.280464 0.280464 1.66533e-16 -2.5 3.3 2.85 3.45082 3.45082 5.77316e-15 -0.5 1.1 2.85 1.51566 1.51566 4.44089e-16 0.5 1.1 2.85 0.560421 0.560421 1.11022e-16 2.5 3.3 2.85 0.0648971 0.0648971 0 0.825 6.7 0.25 223432 223432 5.82077e-11 1.1 6.7 0.25 263079 263079 1.74623e-10 1.65 6.7 0.25 269803 269803 1.16415e-10 3.3 6.7 0.25 82809.3 82809.3 4.36557e-11 0.825 6.7 1.55 26.4657 26.4657 3.55271e-15 1.1 6.7 1.55 28.0935 28.0935 7.10543e-15 1.65 6.7 1.55 23.8892 23.8892 0 3.3 6.7 1.55 4.53388 4.53388 8.88178e-16 0.825 6.7 2.85 3.02245 3.02245 8.88178e-16 1.1 6.7 2.85 2.80407 2.80407 1.33227e-15 1.65 6.7 2.85 1.92626 1.92626 2.22045e-16 3.3 6.7 2.85 0.230205 0.230205 0 R8_CI_TEST: Test R8_CI. X CI(X) R8_CI(X) Diff 0.5 -0.177784 -0.177784 5.55112e-17 0.6 -0.0222707 -0.0222707 7.63278e-17 0.7 0.100515 0.100515 6.93889e-17 0.8 0.198279 0.198279 2.22045e-16 0.9 0.276068 0.276068 5.55112e-17 1 0.337404 0.337404 1.11022e-16 1.2 0.420459 0.420459 5.55112e-17 1.4 0.462007 0.462007 5.55112e-17 1.6 0.471733 0.471733 1.66533e-16 1.8 0.456811 0.456811 1.11022e-16 2 0.422981 0.422981 1.11022e-16 2.5 0.285871 0.285871 1.11022e-16 3 0.11963 0.11963 1.52656e-16 3.5 -0.0321285 -0.0321285 2.77556e-17 4 -0.140982 -0.140982 8.32667e-17 4.5 -0.193491 -0.193491 5.55112e-17 R8_CIN_TEST: Test R8_CIN. X CIN(X) R8_CIN(X) Diff 0.5 0.0618526 0.0618526 6.93889e-18 0.6 0.0886607 0.0886607 1.38778e-17 0.7 0.120026 0.120026 4.16334e-17 0.8 0.155793 0.155793 5.55112e-17 0.9 0.195787 0.195787 0 1 0.239812 0.239812 0 1.2 0.339078 0.339078 5.55112e-17 1.4 0.451681 0.451681 5.55112e-17 1.6 0.575487 0.575487 0 1.8 0.708191 0.708191 0 2 0.847382 0.847382 0 2.5 1.20764 1.20764 4.44089e-16 3 1.5562 1.5562 4.44089e-16 3.5 1.86211 1.86211 2.22045e-16 4 2.10449 2.10449 0 4.5 2.27478 2.27478 4.44089e-16 R8_CINH_TEST: Test R8_CINH. X CINH(X) R8_CINH(X) Diff 0 0 0 0 0.5 0.0631547 0.0631547 1.38778e-17 0.6 0.0913609 0.0913609 1.38778e-17 0.7 0.125028 0.125028 0 0.8 0.164328 0.164328 1.11022e-16 0.9 0.209459 0.209459 2.77556e-17 1 0.260651 0.260651 0 1.2 0.382305 0.382305 1.11022e-16 1.4 0.531806 0.531806 1.11022e-16 1.6 0.712287 0.712287 2.22045e-16 1.8 0.927575 0.927575 0 2 1.1823 1.1823 2.22045e-16 2.5 2.03092 2.03092 4.44089e-16 3 3.28456 3.28456 4.44089e-16 3.5 5.12921 5.12921 8.88178e-16 4 7.85004 7.85004 3.55271e-15 4.5 11.8845 11.8845 5.32907e-15 R8_COS_TEST: Test R8_COS. X COS(X) R8_COS(X) Diff 0 1 1 0 0.261799 0.965926 0.965926 1.11022e-16 0.5 0.877583 0.877583 1.11022e-16 0.523599 0.866025 0.866025 0 0.785398 0.707107 0.707107 0 1 0.540302 0.540302 0 1.0472 0.5 0.5 5.55112e-17 1.5708 0 6.12574e-17 6.12574e-17 2 -0.416147 -0.416147 0 3 -0.989992 -0.989992 0 3.14159 -1 -1 1.11022e-16 4 -0.653644 -0.653644 0 5 0.283662 0.283662 5.55112e-17 R8_COS_DEG_TEST: Test R8_COS_DEG. X COS_DEG(X) R8_COS_DEG(X) Diff -5 0.996195 0.996195 0 0 1 1 0 1 0.999848 0.999848 0 2 0.999391 0.999391 0 3 0.99863 0.99863 0 4 0.997564 0.997564 0 5 0.996195 0.996195 0 10 0.984808 0.984808 0 15 0.965926 0.965926 0 30 0.866025 0.866025 1.11022e-16 45 0.707107 0.707107 0 60 0.5 0.5 1.11022e-16 75 0.258819 0.258819 0 85 0.0871557 0.0871557 4.16334e-17 86 0.0697565 0.0697565 6.93889e-17 87 0.052336 0.052336 1.31839e-16 88 0.0348995 0.0348995 1.11022e-16 89 0.0174524 0.0174524 8.67362e-17 90 0 0 0 91 -0.0174524 -0.0174524 3.46945e-17 105 -0.258819 -0.258819 1.11022e-16 180 -1 -1 0 R8_COSH_TEST: Test R8_COSH X COSH(X) R8_COSH(X) Diff -5 74.2099 74.2099 0 -1 1.54308 1.54308 0 0 1 1 0 0.1 1.005 1.005 0 0.2 1.02007 1.02007 0 0.3 1.04534 1.04534 0 0.4 1.08107 1.08107 2.22045e-16 0.5 1.12763 1.12763 0 0.6 1.18547 1.18547 2.22045e-16 0.7 1.25517 1.25517 0 0.8 1.33743 1.33743 2.22045e-16 0.9 1.43309 1.43309 2.22045e-16 1 1.54308 1.54308 0 2 3.7622 3.7622 0 3 10.0677 10.0677 0 4 27.3082 27.3082 0 5 74.2099 74.2099 0 10 11013.2 11013.2 1.81899e-12 R8_COT_TEST: Test R8_COT. X COT(X) R8_COT(X) Diff 0.0833333 11.9722 11.9722 0 0.261799 3.73205 3.73205 0 0.5 1.83049 1.83049 2.22045e-16 0.523599 1.73205 1.73205 0 0.785398 1 1 0 1 0.642093 0.642093 1.11022e-16 1.0472 0.57735 0.57735 1.11022e-16 1.309 0.267949 0.267949 5.55112e-17 1.5708 0 -0 0 1.4399 0.131652 0.131652 8.32667e-17 1.50535 0.0655435 0.0655435 1.38778e-16 2 -0.457658 -0.457658 1.11022e-16 3 -7.01525 -7.01525 8.88178e-16 4 0.863691 0.863691 0 5 -0.295813 -0.295813 5.55112e-17 R8_CSEVL_TEST: R8_CSEVL evaluates a Chebyshev approximant of N terms at a point X. Here we use an approximant to the exponential function. and average the absolute error at 21 points. N error 1 0.605859 2 0.172657 3 0.0290247 4 0.00362759 5 0.000367891 6 3.07777e-05 7 2.19352e-06 8 1.37199e-07 9 7.67356e-09 10 3.5653e-10 11 1.62784e-11 12 6.53432e-13 R8_DAWSON_TEST: Test R8_DAWSON. X DAWSON(X) R8_DAWSON(X) Diff 0 0 0 0 0.1 0.099336 0.099336 1.38778e-17 0.2 0.194751 0.194751 5.55112e-17 0.3 0.282632 0.282632 5.55112e-17 0.4 0.359943 0.359943 0 0.5 0.424436 0.424436 0 0.6 0.474763 0.474763 5.55112e-17 0.7 0.510504 0.510504 0 0.8 0.532102 0.532102 0 0.9 0.540724 0.540724 1.11022e-16 1 0.53808 0.53808 0 1.1 0.526207 0.526207 1.11022e-16 1.2 0.507273 0.507273 2.22045e-16 1.3 0.483398 0.483398 1.11022e-16 1.4 0.456507 0.456507 0 1.5 0.428249 0.428249 5.55112e-17 1.6 0.39994 0.39994 0 1.7 0.372559 0.372559 5.55112e-17 1.8 0.346773 0.346773 0 1.9 0.322974 0.322974 1.11022e-16 2 0.30134 0.30134 5.55112e-17 R8_E1_TEST: Test R8_E1. X E1(X) R8_E1(X) Diff 0.5 0.559774 0.559774 0 0.6 0.45438 0.45438 0 0.7 0.373769 0.373769 5.55112e-17 0.8 0.310597 0.310597 0 0.9 0.260184 0.260184 5.55112e-17 1 0.219384 0.219384 5.55112e-17 1.1 0.185991 0.185991 5.55112e-17 1.2 0.158408 0.158408 2.77556e-17 1.3 0.135451 0.135451 2.77556e-17 1.4 0.116219 0.116219 1.38778e-17 1.5 0.10002 0.10002 4.16334e-17 1.6 0.0863083 0.0863083 2.77556e-17 1.7 0.0746546 0.0746546 0 1.8 0.0647131 0.0647131 1.38778e-17 1.9 0.0562044 0.0562044 0 2 0.0489005 0.0489005 6.93889e-18 R8_EI_TEST: Test R8_EI. X EI(X) R8_EI(X) Diff 0.5 0.45422 0.45422 0 0.6 0.769881 0.769881 0 0.7 1.06491 1.06491 6.66134e-16 0.8 1.3474 1.3474 0 0.9 1.62281 1.62281 4.44089e-16 1 1.89512 1.89512 2.22045e-16 1.1 2.16738 2.16738 1.33227e-15 1.2 2.44209 2.44209 4.44089e-16 1.3 2.7214 2.7214 1.77636e-15 1.4 3.00721 3.00721 0 1.5 3.30129 3.30129 1.33227e-15 1.6 3.60532 3.60532 8.88178e-16 1.7 3.92096 3.92096 8.88178e-16 1.8 4.24987 4.24987 8.88178e-16 1.9 4.59371 4.59371 8.88178e-16 2 4.95423 4.95423 8.88178e-16 R8_ERF_TEST: Test R8_ERF. X ERF(X) R8_ERF(X) Diff 0 0 0 0 0.1 0.112463 0.112463 0 0.2 0.222703 0.222703 2.77556e-17 0.3 0.328627 0.328627 0 0.4 0.428392 0.428392 5.55112e-17 0.5 0.5205 0.5205 0 0.6 0.603856 0.603856 0 0.7 0.677801 0.677801 1.11022e-16 0.8 0.742101 0.742101 0 0.9 0.796908 0.796908 1.11022e-16 1 0.842701 0.842701 1.11022e-16 1.1 0.880205 0.880205 0 1.2 0.910314 0.910314 1.11022e-16 1.3 0.934008 0.934008 0 1.4 0.952285 0.952285 0 1.5 0.966105 0.966105 1.11022e-16 1.6 0.976348 0.976348 0 1.7 0.98379 0.98379 1.11022e-16 1.8 0.989091 0.989091 1.11022e-16 1.9 0.99279 0.99279 0 2 0.995322 0.995322 0 R8_ERFC_TEST: Test R8_ERFC. X ERFC(X) R8_ERFC(X) Diff 0 1 1 0 0.2 0.777297 0.777297 0 0.4 0.571608 0.571608 1.11022e-16 0.6 0.396144 0.396144 0 0.8 0.257899 0.257899 0 1 0.157299 0.157299 1.11022e-16 1.2 0.089686 0.089686 0 1.4 0.0477149 0.0477149 1.38778e-17 1.6 0.0236516 0.0236516 6.93889e-18 1.8 0.0109095 0.0109095 6.93889e-18 2 0.00467773 0.00467773 8.67362e-19 2.2 0.00186285 0.00186285 1.0842e-18 2.4 0.000688514 0.000688514 2.1684e-19 2.6 0.000236034 0.000236034 1.6263e-19 2.8 7.50132e-05 7.50132e-05 6.77626e-20 3 2.20905e-05 2.20905e-05 6.77626e-21 3.2 6.02576e-06 6.02576e-06 1.35525e-20 3.4 1.52199e-06 1.52199e-06 2.11758e-21 3.6 3.55863e-07 3.55863e-07 3.70577e-22 3.8 7.70039e-08 7.70039e-08 5.29396e-23 4 1.54173e-08 1.54173e-08 3.30872e-24 R8_EXP_TEST: Test R8_EXP. X EXP(X) R8_EXP(X) Diff -10 4.53999e-05 4.53999e-05 0 -5 0.00673795 0.00673795 0 -1 0.367879 0.367879 0 0 1 1 0 1e-08 1 1 0 0.0001 1.0001 1.0001 0 0.001 1.001 1.001 0 0.01 1.01005 1.01005 0 0.1 1.10517 1.10517 0 0.2 1.2214 1.2214 0 0.3 1.34986 1.34986 0 0.4 1.49182 1.49182 0 0.5 1.64872 1.64872 0 0.6 1.82212 1.82212 0 0.7 2.01375 2.01375 4.44089e-16 0.8 2.22554 2.22554 4.44089e-16 0.9 2.4596 2.4596 4.44089e-16 1 2.71828 2.71828 0 2 7.38906 7.38906 0 3.14159 23.1407 23.1407 3.55271e-15 5 148.413 148.413 0 10 22026.5 22026.5 0 20 4.85165e+08 4.85165e+08 0 40 2.35385e+17 2.35385e+17 32 R8_FAC_TEST: R8_FAC evaluates the factorial function. N FAC(N) R8_FAC(N) Diff 0 1 1 0 1 1 1 0 2 2 2 0 3 6 6 0 4 24 24 0 5 120 120 0 6 720 720 0 7 5040 5040 0 8 40320 40320 0 9 362880 362880 0 10 3.6288e+06 3.6288e+06 0 11 3.99168e+07 3.99168e+07 0 12 4.79002e+08 4.79002e+08 0 13 6.22702e+09 6.22702e+09 0 14 8.71783e+10 8.71783e+10 0 15 1.30767e+12 1.30767e+12 0 16 2.09228e+13 2.09228e+13 0 17 3.55687e+14 3.55687e+14 0 18 6.40237e+15 6.40237e+15 0 19 1.21645e+17 1.21645e+17 0 20 2.4329e+18 2.4329e+18 0 25 1.55112e+25 1.55112e+25 4.29497e+09 50 3.04141e+64 3.04141e+64 3.62452e+50 100 9.33262e+157 9.33262e+157 1.90232e+144 150 5.71338e+262 5.71338e+262 4.99922e+249 R8_GAMIC_TEST: R8_GAMIC evaluates the incomplete gamma function. X GAMIC(A,X) R8_GAMIC(A,X) Diff 0.1 0.03 2.4903 2.4903 1.33227e-15 0.1 0.3 0.871837 0.871837 0 0.1 1.5 0.107921 0.107921 2.77556e-17 0.5 0.075 1.23812 1.23812 4.44089e-16 0.5 0.75 0.39113 0.39113 2.22045e-16 0.5 3.5 0.0144472 0.0144472 0 1 0.1 0.904837 0.904837 1.11022e-16 1 1 0.367879 0.367879 5.55112e-17 1 5 0.00673795 0.00673795 0 1.1 0.1 0.882797 0.882797 0 1.1 1 0.390833 0.390833 2.22045e-16 1.1 5 0.00805146 0.00805146 3.46945e-18 2 0.15 0.989814 0.989814 0 2 1.5 0.557825 0.557825 1.11022e-16 2 7 0.00729506 0.00729506 1.73472e-18 6 2.5 114.957 114.957 9.9476e-14 6 12 2.44092 2.44092 2.66454e-15 11 16 280855 280855 5.82077e-11 26 25 8.57648e+24 8.57648e+24 3.22123e+09 41 45 2.08503e+47 2.08503e+47 2.83954e+32 R8_GAMIT_TEST: R8_GAMIT evaluates Tricomi's incomplete Gamma function. X GAMIT(A,X) R8_GAMIT(A,X) Diff 0.1 0.03 1.04829 1.04829 4.44089e-16 0.1 0.3 1.02458 1.02458 2.22045e-16 0.1 1.5 0.949371 0.949371 0 0.5 0.075 1.10079 1.10079 2.22045e-16 0.5 0.75 0.899891 0.899891 1.11022e-16 0.5 3.5 0.530166 0.530166 0 1 0.1 0.951626 0.951626 0 1 1 0.632121 0.632121 1.11022e-16 1 5 0.198652 0.198652 2.77556e-17 1.1 0.1 0.907178 0.907178 0 1.1 1 0.589181 0.589181 0 1.1 5 0.168827 0.168827 5.55112e-17 2 0.15 0.452703 0.452703 0 2 1.5 0.196522 0.196522 0 2 7 0.0202593 0.0202593 6.93889e-18 6 2.5 0.000172118 0.000172118 8.13152e-20 6 12 3.28086e-07 3.28086e-07 3.17637e-22 11 16 5.2444e-14 5.2444e-14 9.46633e-29 26 25 2.01346e-37 2.01346e-37 4.17619e-53 41 45 1.23062e-68 1.23062e-68 1.77076e-82 R8_GAML_TEST: R8_GAML returns bounds for the argument of the gamma function. Lower limit XMIN = -170.433 Upper limit XMAX = 171.48 R8_GAMMA_TEST: R8_GAMMA evaluates the gamma function. X GAMMA(X) R8_GAMMA(X) Diff -0.5 -3.54491 -3.54491 0 -0.01 -100.587 -100.587 2.84217e-14 0.01 99.4326 99.4326 1.42109e-14 0.1 9.51351 9.51351 1.77636e-15 0.2 4.59084 4.59084 0 0.4 2.21816 2.21816 0 0.5 1.77245 1.77245 0 0.6 1.48919 1.48919 0 0.8 1.16423 1.16423 2.22045e-16 1 1 1 0 1.1 0.951351 0.951351 0 1.2 0.918169 0.918169 0 1.3 0.897471 0.897471 0 1.4 0.887264 0.887264 0 1.5 0.886227 0.886227 0 1.6 0.893515 0.893515 0 1.7 0.908639 0.908639 0 1.8 0.931384 0.931384 0 1.9 0.961766 0.961766 0 2 1 1 0 3 2 2 0 4 6 6 0 10 362880 362880 0 20 1.21645e+17 1.21645e+17 80 30 8.84176e+30 8.84176e+30 7.54353e+16 R8_GAMR_TEST: R8_GAMR evaluates 1.0/Gamma(x). X 1/GAMMA(X) R8_GAMR(X) Diff -0.5 -0.282095 -0.282095 0 -0.01 -0.00994162 -0.00994162 3.46945e-18 0.01 0.0100571 0.0100571 1.73472e-18 0.1 0.105114 0.105114 1.38778e-17 0.2 0.217825 0.217825 0 0.4 0.450824 0.450824 0 0.5 0.56419 0.56419 0 0.6 0.671505 0.671505 0 0.8 0.858937 0.858937 1.11022e-16 1 1 1 0 1.1 1.05114 1.05114 0 1.2 1.08912 1.08912 0 1.3 1.11424 1.11424 0 1.4 1.12706 1.12706 0 1.5 1.12838 1.12838 0 1.6 1.11917 1.11917 0 1.7 1.10055 1.10055 0 1.8 1.07367 1.07367 0 1.9 1.03975 1.03975 0 2 1 1 0 3 0.5 0.5 0 4 0.166667 0.166667 0 10 2.75573e-06 2.75573e-06 0 20 8.22064e-18 8.22064e-18 4.62223e-33 30 1.131e-31 1.131e-31 9.63393e-46 R8_INITS_TEST: R8_INITS determines the Chebyshev interpolant degree necessary to guarantee a desired accuracy level. Here, we use a 15 term Chebyshev expansion for the sine function. Accuracy Terms Needed 1 1 0.1 2 0.01 2 0.001 3 0.0001 3 1e-05 4 1e-06 4 1e-07 5 1e-08 5 1e-09 6 1e-10 6 1e-11 7 1e-12 7 1e-13 7 1e-14 8 1e-15 8 1e-16 9 1e-17 9 R8_INT_TEST: R8_INT rounds an R8 to an integer value. X INT(X) R8_INT(X) Diff -2.01 -2 -2 0 -1.99 -1 -1 0 -1.5 -1 -1 0 -1.1 -1 -1 0 -1.01 -1 -1 0 -1 -1 -1 0 -0.99 0 0 0 -0.9 0 0 0 -0.51 0 0 0 -0.5 0 0 0 -0.49 0 0 0 -0.01 0 0 0 0 0 0 0 0.01 0 0 0 0.49 0 0 0 0.5 0 0 0 0.51 0 0 0 0.9 0 0 0 0.99 0 0 0 1 1 1 0 1.01 1 1 0 1.1 1 1 0 1.5 1 1 0 1.99 1 1 0 2.01 2 2 0 R8_LBETA_TEST: R8_LBETA evaluates the logarithm of the Beta function. X LBETA(A,B) R8_LBETA(A,B) Diff 0.2 1 1.60944 1.60944 2.22045e-16 0.4 1 0.916291 0.916291 0 0.6 1 0.510826 0.510826 0 0.8 1 0.223144 0.223144 2.22045e-16 1 0.2 1.60944 1.60944 2.22045e-16 1 0.4 0.916291 0.916291 0 1 1 0 0 0 2 2 -1.79176 -1.79176 0 3 3 -3.4012 -3.4012 4.44089e-16 4 4 -4.94164 -4.94164 0 5 5 -6.44572 -6.44572 8.88178e-16 6 2 -3.73767 -3.73767 4.44089e-16 6 3 -5.12396 -5.12396 0 6 4 -6.22258 -6.22258 0 6 5 -7.13887 -7.13887 2.66454e-15 6 6 -7.92732 -7.92732 4.44089e-15 7 7 -9.39366 -9.39366 1.77636e-15 R8_LGAMS_TEST: R8_LGAMS evaluates the sign of Gamma(x) and the logarithm of the absolute value of Gamma(x). X LNGAM(X) Sign(Gamma(X)) ALNGAM Diff 0.2 1.52406 1 1.52406 4.44089e-16 0.4 0.796678 1 0.796678 1.11022e-16 0.6 0.398234 1 0.398234 1.11022e-16 0.8 0.15206 1 0.15206 0 1 0 1 0 0 1.1 -0.0498724 1 -0.0498724 1.38778e-17 1.2 -0.0853741 1 -0.0853741 2.77556e-17 1.3 -0.108175 1 -0.108175 6.93889e-17 1.4 -0.119613 1 -0.119613 1.38778e-16 1.5 -0.120782 1 -0.120782 2.77556e-17 1.6 -0.112592 1 -0.112592 6.93889e-17 1.7 -0.0958077 1 -0.0958077 1.38778e-17 1.8 -0.0710839 1 -0.0710839 1.38778e-17 1.9 -0.0389843 1 -0.0389843 1.38778e-17 2 0 1 0 0 3 0.693147 1 0.693147 0 4 1.79176 1 1.79176 0 10 12.8018 1 12.8018 0 20 39.3399 1 39.3399 0 30 71.257 1 71.257 1.42109e-14 R8_LGMC_TEST: R8_LGMC evaluates the correction log gamma factor. r8_lgmc(x) = log ( gamma ( x ) ) - log ( sqrt ( 2 * pi ) - ( x - 0.5 ) * log ( x ) + x X LGMC(X) R8_LGMC(X) Diff 10 0.00833056 0.00833056 1.62197e-15 20 0.00416632 0.00416632 6.67869e-17 30 0.00277767 0.00277767 7.41638e-15 R8_LI_TEST: R8_LI evaluates the logarithmic integral X LI(X) R8_LI(X) Diff 0 0 0 0 0.1 -0.0323898 -0.0323898 6.93889e-18 0.2 -0.0851265 -0.0851265 1.38778e-17 0.3 -0.157415 -0.157415 5.55112e-17 0.4 -0.252949 -0.252949 0 0.5 -0.378671 -0.378671 5.55112e-17 0.6 -0.546851 -0.546851 1.11022e-16 0.7 -0.780947 -0.780947 2.22045e-16 0.8 -1.13401 -1.13401 2.22045e-16 0.9 -1.7758 -1.7758 2.22045e-16 0.95 -2.44362 -2.44362 4.44089e-16 0.975 -3.12419 -3.12419 8.88178e-16 1.03125 -2.87294 -2.87294 4.44089e-16 1.0625 -2.16428 -2.16428 4.44089e-16 1.125 -1.44035 -1.44035 2.22045e-16 1.25 -0.686488 -0.686488 0 1.5 0.125065 0.125065 0 2 1.04516 1.04516 0 4 2.96759 2.96759 0 8 5.25372 5.25372 8.88178e-16 16 8.51972 8.51972 0 32 13.6051 13.6051 3.55271e-15 64 21.9347 21.9347 1.06581e-14 128 36.0425 36.0425 7.10543e-15 256 60.5131 60.5131 7.10543e-15 512 103.721 103.721 2.84217e-14 1024 181.078 181.078 2.84217e-14 2048 321.114 321.114 0 R8_LNGAM_TEST: Test R8_LNGAM X LNGAM(X) R8_LNGAM(X) Diff 0.2 1.52406 1.52406 4.44089e-16 0.4 0.796678 0.796678 1.11022e-16 0.6 0.398234 0.398234 1.11022e-16 0.8 0.15206 0.15206 0 1 0 0 0 1.1 -0.0498724 -0.0498724 1.38778e-17 1.2 -0.0853741 -0.0853741 2.77556e-17 1.3 -0.108175 -0.108175 6.93889e-17 1.4 -0.119613 -0.119613 1.38778e-16 1.5 -0.120782 -0.120782 2.77556e-17 1.6 -0.112592 -0.112592 6.93889e-17 1.7 -0.0958077 -0.0958077 1.38778e-17 1.8 -0.0710839 -0.0710839 1.38778e-17 1.9 -0.0389843 -0.0389843 1.38778e-17 2 0 0 0 3 0.693147 0.693147 0 4 1.79176 1.79176 0 10 12.8018 12.8018 0 20 39.3399 39.3399 0 30 71.257 71.257 1.42109e-14 R8_LNREL_TEST: R8_LNREL evaluates ln(1+x) X LOG(1+X) R8_LNREL(X) Diff -0.99999 -11.5129 -11.5129 4.55103e-12 -0.99 -4.60517 -4.60517 8.88178e-16 -0.9 -2.30259 -2.30259 0 -0.8 -1.60944 -1.60944 2.22045e-16 -0.7 -1.20397 -1.20397 0 -0.6 -0.916291 -0.916291 1.11022e-16 -0.5 -0.693147 -0.693147 0 -0.4 -0.510826 -0.510826 0 -0.3 -0.356675 -0.356675 1.11022e-16 -0.2 -0.223144 -0.223144 5.55112e-17 -0.1 -0.105361 -0.105361 2.77556e-17 0 0 0 0 1 0.693147 0.693147 0 2 1.09861 1.09861 0 2.14159 1.14473 1.14473 0 4 1.60944 1.60944 0 9 2.30259 2.30259 0 19 2.99573 2.99573 0 99 4.60517 4.60517 0 1.23457e+08 18.6314 18.6314 0 R8_LOG_TEST: R8_LOG evaluates ln(x) X LOG(X) R8_LOG(X) Diff 1e-05 -11.5129 -11.5129 0 0.01 -4.60517 -4.60517 0 0.1 -2.30259 -2.30259 4.44089e-16 0.2 -1.60944 -1.60944 0 0.3 -1.20397 -1.20397 2.22045e-16 0.4 -0.916291 -0.916291 1.11022e-16 0.5 -0.693147 -0.693147 0 0.6 -0.510826 -0.510826 0 0.7 -0.356675 -0.356675 5.55112e-17 0.8 -0.223144 -0.223144 1.11022e-16 0.9 -0.105361 -0.105361 6.93889e-17 1 0 0 0 2 0.693147 0.693147 0 3 1.09861 1.09861 2.22045e-16 3.14159 1.14473 1.14473 0 5 1.60944 1.60944 2.22045e-16 10 2.30259 2.30259 0 20 2.99573 2.99573 0 100 4.60517 4.60517 0 1.23457e+08 18.6314 18.6314 0 R8_LOG10_TEST: Test R8_LOG10 X LOG10(X) R8_LOG10(X) Diff 1e-05 -5 -5 0 0.01 -2 -2 2.22045e-16 0.1 -1 -1 1.11022e-16 0.2 -0.69897 -0.69897 1.11022e-16 0.3 -0.522879 -0.522879 1.11022e-16 0.4 -0.39794 -0.39794 5.55112e-17 0.5 -0.30103 -0.30103 0 0.6 -0.221849 -0.221849 0 0.7 -0.154902 -0.154902 2.77556e-17 0.8 -0.09691 -0.09691 2.77556e-17 0.9 -0.0457575 -0.0457575 6.93889e-18 1 0 0 0 2 0.30103 0.30103 0 3 0.477121 0.477121 5.55112e-17 3.14159 0.49715 0.49715 0 5 0.69897 0.69897 1.11022e-16 10 1 1 0 20 1.30103 1.30103 2.22045e-16 100 2 2 0 1.23457e+08 8.09151 8.09151 0 R8_MACH_TEST: R8_MACH evaluates double precision machine numbers. R8_MACH (1) = B^(EMIN-1), the smallest positive magnitude. R8_MACH (2) = B^EMAX*(1 - B^(-T)), the largest magnitude. R8_MACH (3) = B^(-T), the smallest relative spacing. R8_MACH (4) = B^(1-T), the largest relative spacing. R8_MACH (5) = LOG10(B) I R8_MACH(I) 1 4.45015e-308 2 8.98847e+307 3 1.11022e-16 4 2.22045e-16 5 0.30103 R8_PAK_TEST: R8_PAK converts a mantissa and base 2 exponent to an R8. Mantissa Exponent R8 0.5 7 64 0.5 8 128 -0.5 7 -64 0.75 7 96 0.9375 4 15 0.5 0 0.5 0.5 -1 0.25 0.625 0 0.625 0.504883 7 64.625 0.785398 2 3.14159 0 0 0 R8_POCH_TEST: R8_POCH evaluates the Pochhammer symbol.. X POCH(A,X) R8_POCH(A,X) Diff 5 4 1680 1680 0 5.25 4 1962.6 1962.6 0 5.5 4 2279.06 2279.06 0 5.75 4 2631.97 2631.97 0 6 4 3024 3024 0 7.5 0 1 1 0 7.5 1 7.5 7.5 0 7.5 2 63.75 63.75 0 7.5 3 605.625 605.625 0 7.5 4 6359.06 6359.06 0 7.5 5 73129.2 73129.2 0 7.5 6 914115 914115 0 7.5 7 1.23406e+07 1.23406e+07 0 7.5 8 1.78938e+08 1.78938e+08 5.96046e-08 7.5 9 2.77354e+09 2.77354e+09 0 R8_PSI_TEST: Test R8_PSI X PSI(X) R8_PSI(X) Diff 0.1 -10.4238 -10.4238 1.77636e-15 0.2 -5.28904 -5.28904 8.88178e-16 0.3 -3.50252 -3.50252 0 0.4 -2.56138 -2.56138 0 0.5 -1.96351 -1.96351 4.44089e-16 0.6 -1.54062 -1.54062 4.44089e-16 0.7 -1.22002 -1.22002 2.22045e-16 0.8 -0.965009 -0.965009 0 0.9 -0.754927 -0.754927 0 1 -0.577216 -0.577216 0 1.1 -0.423755 -0.423755 1.11022e-16 1.2 -0.28904 -0.28904 1.11022e-16 1.3 -0.169191 -0.169191 1.38778e-16 1.4 -0.0613845 -0.0613845 6.93889e-17 1.5 0.03649 0.03649 0 1.6 0.126047 0.126047 2.77556e-17 1.7 0.208548 0.208548 5.55112e-17 1.8 0.284991 0.284991 0 1.9 0.356184 0.356184 1.11022e-16 2 0.422784 0.422784 5.55112e-17 R8_RAND_TEST: R8_RAND returns a random R8 value. I R8_RAND Expected 1 0.000412703 0.000412703 2 0.675084 0.675084 3 0.161475 0.161475 4 0.90862 0.90862 10 0.552779 0.552779 100 0.360089 0.360089 1000 0.217699 0.217699 Average = 0.4999 0.5 Variance = 0.0832492 0.0833333 R8_RANDGS_TEST: R8_RANDGS generates a random normal R8. Mean = 3 Standard deviation = 2 I R8_RANDGS 1 4.40826 2 2.94895 3 2.13857 4 1.57902 5 8.76356 6 1.08396 7 3.87796 8 4.96727 9 0.673294 10 2.13816 Sequence mean = 2.98999 Sequence standard deviation = 2.03537 R8_RANDOM_TEST: R8_RANDOM returns a random R8 value. I R8_RANDOM 1 0.0913844 2 0.161475 3 0.658159 4 0.272802 10 0.174952 100 0.820123 1000 0.640682 Average = 0.499899 0.5 Variance = 0.0832494 0.0833333 R8_REN_TEST: R8_REN returns a random R8 value. I R8_REN Expected 1 0.470393 0.470393 2 0.799066 0.799066 3 0.883261 0.883261 4 0.407667 0.407667 10 0.955566 0.955566 100 0.173576 0.173576 1000 0.0121733 0.0121733 Average = 0.499767 0.5 Variance = 0.0833347 0.0833333 R8_SHI_TEST: Test R8_SHI. X SHI(X) R8_SHI(X) Diff 0.5 0.506997 0.506997 0 0.6 0.61213 0.61213 0 0.7 0.719338 0.719338 1.11022e-16 0.8 0.828997 0.828997 0 0.9 0.941498 0.941498 0 1 1.05725 1.05725 4.44089e-16 1.2 1.30025 1.30025 0 1.4 1.56171 1.56171 0 1.6 1.84581 1.84581 6.66134e-16 1.8 2.15729 2.15729 8.88178e-16 2 2.50157 2.50157 0 2.5 3.54934 3.54934 0 3 4.97344 4.97344 0 3.5 6.96616 6.96616 0 4 9.81733 9.81733 0 4.5 13.9679 13.9679 1.77636e-15 R8_SI_TEST: Test R8_SI. X SI(X) R8_SI(X) Diff 0.5 0.493107 0.493107 5.55112e-17 0.6 0.588129 0.588129 1.11022e-16 0.7 0.681222 0.681222 0 0.8 0.772096 0.772096 0 0.9 0.860471 0.860471 0 1 0.946083 0.946083 0 1.2 1.10805 1.10805 4.44089e-16 1.4 1.25623 1.25623 0 1.6 1.38918 1.38918 6.66134e-16 1.8 1.50582 1.50582 4.44089e-16 2 1.60541 1.60541 2.22045e-16 2.5 1.77852 1.77852 2.22045e-16 3 1.84865 1.84865 2.22045e-16 3.5 1.83313 1.83313 0 4 1.7582 1.7582 0 4.5 1.65414 1.65414 0 R8_SIN_TEST: Test R8_SIN. X SIN(X) R8_SIN(X) Diff 0 0 0 0 0.261799 0.258819 0.258819 5.55112e-17 0.5 0.479426 0.479426 0 0.523599 0.5 0.5 0 0.785398 0.707107 0.707107 1.11022e-16 1 0.841471 0.841471 0 1.0472 0.866025 0.866025 1.11022e-16 1.5708 1 1 1.11022e-16 2 0.909297 0.909297 0 3 0.14112 0.14112 0 3.14159 0 1.22515e-16 1.22515e-16 4 -0.756802 -0.756802 0 5 -0.958924 -0.958924 1.11022e-16 R8_SIN_DEG_TEST: Test R8_SIN_DEG. X SIN_DEG(X) R8_SIN_DEG(X) Diff -5 -0.0871557 -0.0871557 1.38778e-17 0 0 0 0 1 0.0174524 0.0174524 0 2 0.0348995 0.0348995 0 3 0.052336 0.052336 0 4 0.0697565 0.0697565 0 5 0.0871557 0.0871557 1.38778e-17 10 0.173648 0.173648 2.77556e-17 15 0.258819 0.258819 0 30 0.5 0.5 5.55112e-17 45 0.707107 0.707107 1.11022e-16 60 0.866025 0.866025 0 75 0.965926 0.965926 0 85 0.996195 0.996195 0 86 0.997564 0.997564 0 87 0.99863 0.99863 0 88 0.999391 0.999391 0 89 0.999848 0.999848 0 90 1 1 0 91 0.999848 0.999848 0 105 0.965926 0.965926 0 180 0 0 0 R8_SINH_TEST: Test R8_SINH X SINH(X) R8_SINH(X) Diff -5 -74.2032 -74.2032 0 -1 -1.1752 -1.1752 0 0 0 0 0 0.1 0.100167 0.100167 1.38778e-17 0.2 0.201336 0.201336 2.77556e-17 0.3 0.30452 0.30452 0 0.4 0.410752 0.410752 0 0.5 0.521095 0.521095 0 0.6 0.636654 0.636654 1.11022e-16 0.7 0.758584 0.758584 1.11022e-16 0.8 0.888106 0.888106 0 0.9 1.02652 1.02652 0 1 1.1752 1.1752 0 2 3.62686 3.62686 4.44089e-16 3 10.0179 10.0179 0 4 27.2899 27.2899 3.55271e-15 5 74.2032 74.2032 0 10 11013.2 11013.2 0 R8_SPENCE_TEST: Test R8_SPENCE X SPENCE(X) R8_SPENCE(X) Diff 0 0 0 0 0.05 0.0506393 0.0506393 0 0.1 0.102618 0.102618 2.77556e-17 0.15 0.156035 0.156035 0 0.2 0.211004 0.211004 2.77556e-17 0.25 0.267653 0.267653 0 0.3 0.32613 0.32613 0 0.35 0.386606 0.386606 0 0.4 0.449283 0.449283 5.55112e-17 0.45 0.514399 0.514399 1.11022e-16 0.5 0.582241 0.582241 0 0.55 0.653158 0.653158 0 0.6 0.727586 0.727586 0 0.65 0.806083 0.806083 0 0.7 0.889378 0.889378 2.22045e-16 0.75 0.978469 0.978469 1.11022e-16 0.8 1.07479 1.07479 4.44089e-16 0.85 1.18058 1.18058 0 0.9 1.29971 1.29971 2.22045e-16 0.95 1.44063 1.44063 4.44089e-16 1 1.64493 1.64493 4.44089e-16 R8_SQRT_TEST: Test R8_SQRT X SQRT(X) R8_SQRT(X) Diff 0 0 0 0 8.1e-09 9e-05 9e-05 1.35525e-20 0.09 0.3 0.3 5.55112e-17 0.1 0.316228 0.316228 1.11022e-16 0.4 0.632456 0.632456 1.11022e-16 1 1 1 2.22045e-16 2 1.41421 1.41421 2.22045e-16 3 1.73205 1.73205 2.22045e-16 3.14159 1.77245 1.77245 2.22045e-16 19 4.3589 4.3589 0 29 5.38516 5.38516 0 71 8.42615 8.42615 0 97 9.84886 9.84886 0 1.23457e+08 11111.1 11111.1 5.45697e-12 R8_TAN_TEST: Test R8_TAN. X TAN(X) R8_TAN(X) Diff 0 0 0 0 0.261799 0.267949 0.267949 5.55112e-17 0.5 0.546302 0.546302 0 0.523599 0.57735 0.57735 2.22045e-16 0.785398 1 1 0 1 1.55741 1.55741 2.22045e-16 1.0472 1.73205 1.73205 1.11022e-15 1.309 3.73205 3.73205 8.88178e-16 1.4399 7.59575 7.59575 7.99361e-15 1.50535 15.2571 15.2571 1.95399e-14 2 -2.18504 -2.18504 4.44089e-16 3 -0.142547 -0.142547 2.77556e-17 3.14159 0 -1.11022e-16 1.11022e-16 4 1.15782 1.15782 4.44089e-16 5 -3.38052 -3.38052 0 R8_TANH_TEST: Test R8_TANH X TANH(X) R8_TANH(X) Diff -5 -0.999909 -0.999909 1.11022e-16 -1 -0.761594 -0.761594 0 0 0 0 0 0.1 0.099668 0.099668 0 0.2 0.197375 0.197375 0 0.3 0.291313 0.291313 0 0.4 0.379949 0.379949 0 0.5 0.462117 0.462117 5.55112e-17 0.6 0.53705 0.53705 1.11022e-16 0.7 0.604368 0.604368 0 0.8 0.664037 0.664037 1.11022e-16 0.9 0.716298 0.716298 0 1 0.761594 0.761594 0 2 0.964028 0.964028 1.11022e-16 3 0.995055 0.995055 1.11022e-16 4 0.999329 0.999329 1.11022e-16 5 0.999909 0.999909 1.11022e-16 10 1 1 1.11022e-16 R8_UPAK_TEST: R8_UPAK converts an R8 to a mantissa and base 2 exponent. X Mantissa Exponent 64 0.5 7 128 0.5 8 -64 -0.5 7 96 0.75 7 15 0.9375 4 0.5 0.5 0 0.25 0.5 -1 0.625 0.625 0 64.625 0.504883 7 3.14159 0.785398 2 0 0 0 FN_TEST: Normal end of execution. 04 March 2022 02:13:07 PM