9 January 2016 11:14:35.346 AM TOMS597_PRB FORTRAN77 version Test the TOMS597 library. 1TEST OF I(X,ALPHA) VS SERIES EXPANSION 2000 RANDOM ARGUMENTS WERE TESTED FROM THE INTERVAL ( 0.0, 2.0) I(X,ALPHA) WAS LARGER 760 TIMES, AGREED 562 TIMES, AND WAS SMALLER 678 TIMES. THERE ARE 53 BASE 2 SIGNIFICANT DIGITS IN A FLOATING-POINT NUMBER THE MAXIMUM RELATIVE ERROR OF 0.2124E-14 = 2 ** -48.74 OCCURRED FOR X = 0.172002E+01 AND NU = 0.333546E-01 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 4.26 THE ROOT MEAN SQUARE RELATIVE ERROR WAS 0.5232E-15 = 2 ** -50.76 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 2.24 1TEST OF I(X,ALPHA) VS SERIES EXPANSION 2000 RANDOM ARGUMENTS WERE TESTED FROM THE INTERVAL ( 2.0, 4.0) I(X,ALPHA) WAS LARGER 854 TIMES, AGREED 360 TIMES, AND WAS SMALLER 786 TIMES. THERE ARE 53 BASE 2 SIGNIFICANT DIGITS IN A FLOATING-POINT NUMBER THE MAXIMUM RELATIVE ERROR OF 0.2292E-14 = 2 ** -48.63 OCCURRED FOR X = 0.367452E+01 AND NU = 0.125949E+00 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 4.37 THE ROOT MEAN SQUARE RELATIVE ERROR WAS 0.8137E-15 = 2 ** -50.13 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 2.87 1TEST OF I(X,ALPHA) VS SERIES EXPANSION 2000 RANDOM ARGUMENTS WERE TESTED FROM THE INTERVAL ( 4.0, 10.0) I(X,ALPHA) WAS LARGER 886 TIMES, AGREED 274 TIMES, AND WAS SMALLER 840 TIMES. THERE ARE 53 BASE 2 SIGNIFICANT DIGITS IN A FLOATING-POINT NUMBER THE MAXIMUM RELATIVE ERROR OF 0.4425E-14 = 2 ** -47.68 OCCURRED FOR X = 0.893431E+01 AND NU = 0.491335E-01 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 5.32 THE ROOT MEAN SQUARE RELATIVE ERROR WAS 0.1381E-14 = 2 ** -49.36 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 3.64 1TEST OF EXP(-X)*I(X,0.5) VS 1/SQRT(2*X*PI) 2000 RANDOM ARGUMENTS WERE TESTED FROM THE INTERVAL ( 20.0, 30.0) I(X,ALPHA) WAS LARGER 267 TIMES, AGREED 612 TIMES, AND WAS SMALLER 1121 TIMES. THERE ARE 53 BASE 2 SIGNIFICANT DIGITS IN A FLOATING-POINT NUMBER THE MAXIMUM RELATIVE ERROR OF 0.7387E-15 = 2 ** -50.27 OCCURRED FOR X = 0.281843E+02 AND NU = 0.500000E+00 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 2.73 THE ROOT MEAN SQUARE RELATIVE ERROR WAS 0.2156E-15 = 2 ** -52.04 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 0.96 1CHECK OF ERROR RETURNS THE FOLLOWING SUMMARIZES CALLS WITH INDICATED PARAMETERS NCALC DIFFERENT FROM NB INDICATES SOME FORM OF ERROR SEE DOCUMENTATION FOR RIBESL FOR DETAILS ARG ALPHA NB IZ RES NCALC 0.1000000E+01 0.5000000E+00 5 1 0.9376749E+00 5 -0.1000000E+01 0.5000000E+00 5 1 0.0000000E+00 -1 0.1000000E+01 0.1500000E+01 5 1 0.0000000E+00 -1 0.1000000E+01 0.5000000E+00 -5 1 0.0000000E+00 -6 0.1000000E+01 0.5000000E+00 5 5 0.0000000E+00 -1 0.4940656-323 0.0000000E+00 150 1 0.1000000E+01 1 0.3000000E+02 0.0000000E+00 150 1 0.7816723E+12 78 +Infinity 0.0000000E+00 150 1 0.0000000E+00 -1 1TEST OF GAMMA(X) VS DUPLICATION FORMULA 2000 RANDOM ARGUMENTS WERE TESTED FROM THE INTERVAL ( 0.0, 1.0) GAMMA(X) WAS LARGER 699 TIMES, AGREED 740 TIMES, AND WAS SMALLER 561 TIMES. THERE ARE 53 BASE 2 SIGNIFICANT DIGITS IN A FLOATING-POINT NUMBER THE MAXIMUM RELATIVE ERROR OF 0.1035E-14 = 2 ** -49.78 OCCURRED FOR X = 0.120583E-01 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 3.22 THE ROOT MEAN SQUARE RELATIVE ERROR WAS 0.2035E-15 = 2 ** -52.13 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 0.87 1TEST OF GAMMA(X) VS DUPLICATION FORMULA 2000 RANDOM ARGUMENTS WERE TESTED FROM THE INTERVAL ( 1.0, 2.0) GAMMA(X) WAS LARGER 706 TIMES, AGREED 782 TIMES, AND WAS SMALLER 512 TIMES. THERE ARE 53 BASE 2 SIGNIFICANT DIGITS IN A FLOATING-POINT NUMBER THE MAXIMUM RELATIVE ERROR OF 0.3759E-15 = 2 ** -51.24 OCCURRED FOR X = 0.142779E+01 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 1.76 THE ROOT MEAN SQUARE RELATIVE ERROR WAS 0.1204E-15 = 2 ** -52.88 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 0.12 1TEST OF GAMMA(X) VS DUPLICATION FORMULA 2000 RANDOM ARGUMENTS WERE TESTED FROM THE INTERVAL ( 2.0, 10.0) GAMMA(X) WAS LARGER 802 TIMES, AGREED 458 TIMES, AND WAS SMALLER 740 TIMES. THERE ARE 53 BASE 2 SIGNIFICANT DIGITS IN A FLOATING-POINT NUMBER THE MAXIMUM RELATIVE ERROR OF 0.3368E-14 = 2 ** -48.08 OCCURRED FOR X = 0.859977E+01 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 4.92 THE ROOT MEAN SQUARE RELATIVE ERROR WAS 0.7099E-15 = 2 ** -50.32 THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 2.68 1TEST OF GAMMA(X) VS DUPLICATION FORMULA 2000 RANDOM ARGUMENTS WERE TESTED FROM THE INTERVAL ( 10.0, NaN) GAMMA(X) WAS LARGER 0 TIMES, AGREED 2000 TIMES, AND WAS SMALLER 0 TIMES. THERE ARE 53 BASE 2 SIGNIFICANT DIGITS IN A FLOATING-POINT NUMBER THE MAXIMUM RELATIVE ERROR OF NaN = 2 ** NaN OCCURRED FOR X = NaN THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 0.00 THE ROOT MEAN SQUARE RELATIVE ERROR WAS NaN = 2 ** NaN THE ESTIMATED LOSS OF BASE 2 SIGNIFICANT DIGITS IS 0.00 1SPECIAL TESTS TEST OF SPECIAL ARGUMENTS GAMMA (-0.500000E+00) = -0.354491E+01 GAMMA ( 0.494066-323) = +Infinity GAMMA ( 0.100000E+01) = 0.100000E+01 GAMMA ( 0.200000E+01) = 0.100000E+01 GAMMA ( NaN) = NaN 1TEST OF ERROR RETURNS GAMMA WILL BE CALLED WITH THE ARGUMENT-0.100000E+01 THIS SHOULD TRIGGER AN ERROR MESSAGE GAMMA RETURNED THE VALUE NaN GAMMA WILL BE CALLED WITH THE ARGUMENT 0.000000E+00 THIS SHOULD TRIGGER AN ERROR MESSAGE GAMMA RETURNED THE VALUE +Infinity GAMMA WILL BE CALLED WITH THE ARGUMENT 0.494066-323 THIS SHOULD TRIGGER AN ERROR MESSAGE GAMMA RETURNED THE VALUE +Infinity GAMMA WILL BE CALLED WITH THE ARGUMENT NaN THIS SHOULD TRIGGER AN ERROR MESSAGE GAMMA RETURNED THE VALUE NaN THIS CONCLUDES THE TESTS