27 May 2021 03:02:01 PM TOMS743_TEST C version Test the TOMS743 library. Number of bits in mantissa - 1 = 52 TEST01 Compare WAPR(X) to stored values. Wp results for x near -exp(-1) Offset x W(x) (WAPR) W(x) (EXACT) Digits Correct 1.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 2.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 3.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 4.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 5.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 6.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 7.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 8.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 9.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 1.00000000e-39 -1.00000000e+00 -1.00000000e+00 16 1.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 2.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 3.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 4.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 5.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 6.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 7.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 8.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 9.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 1.00000000e-29 -1.00000000e+00 -1.00000000e+00 16 1.00000000e-20 -1.00000000e+00 -1.00000000e+00 16 2.00000000e-20 -1.00000000e+00 -1.00000000e+00 16 3.00000000e-20 -1.00000000e+00 -1.00000000e+00 16 4.00000000e-20 -1.00000000e+00 -1.00000000e+00 16 5.00000000e-20 -9.99999999e-01 -9.99999999e-01 16 6.00000000e-20 -9.99999999e-01 -9.99999999e-01 16 7.00000000e-20 -9.99999999e-01 -9.99999999e-01 16 8.00000000e-20 -9.99999999e-01 -9.99999999e-01 16 9.00000000e-20 -9.99999999e-01 -9.99999999e-01 16 1.00000000e-19 -9.99999999e-01 -9.99999999e-01 16 1.00000000e-10 -9.99976684e-01 -9.99976684e-01 16 2.00000000e-10 -9.99967026e-01 -9.99967026e-01 16 3.00000000e-10 -9.99959615e-01 -9.99959615e-01 16 4.00000000e-10 -9.99953368e-01 -9.99953368e-01 16 5.00000000e-10 -9.99947864e-01 -9.99947864e-01 16 6.00000000e-10 -9.99942888e-01 -9.99942888e-01 16 7.00000000e-10 -9.99938312e-01 -9.99938312e-01 16 8.00000000e-10 -9.99934053e-01 -9.99934053e-01 16 9.00000000e-10 -9.99930052e-01 -9.99930052e-01 16 1.00000000e-09 -9.99926269e-01 -9.99926269e-01 16 1.00000000e-05 -9.92644755e-01 -9.92644755e-01 16 2.00000000e-05 -9.89608643e-01 -9.89608643e-01 16 3.00000000e-05 -9.87283109e-01 -9.87283109e-01 16 4.00000000e-05 -9.85325390e-01 -9.85325390e-01 16 5.00000000e-05 -9.83602718e-01 -9.83602718e-01 16 6.00000000e-05 -9.82047003e-01 -9.82047003e-01 16 7.00000000e-05 -9.80617797e-01 -9.80617797e-01 16 8.00000000e-05 -9.79288746e-01 -9.79288746e-01 16 9.00000000e-05 -9.78041545e-01 -9.78041545e-01 16 1.00000000e-04 -9.76862866e-01 -9.76862866e-01 16 2.00000000e-04 -9.67382627e-01 -9.67382627e-01 15 3.00000000e-04 -9.60148542e-01 -9.60148542e-01 15 4.00000000e-04 -9.54076869e-01 -9.54076869e-01 15 5.00000000e-04 -9.48747869e-01 -9.48747869e-01 15 6.00000000e-04 -9.43946291e-01 -9.43946291e-01 15 7.00000000e-04 -9.39544278e-01 -9.39544278e-01 15 8.00000000e-04 -9.35458531e-01 -9.35458531e-01 15 9.00000000e-04 -9.31631195e-01 -9.31631195e-01 15 1.00000000e-03 -9.28020150e-01 -9.28020150e-01 15 2.00000000e-03 -8.99185766e-01 -8.99185766e-01 15 3.00000000e-03 -8.77428717e-01 -8.77428717e-01 15 4.00000000e-03 -8.59327504e-01 -8.59327504e-01 16 5.00000000e-03 -8.43558002e-01 -8.43558002e-01 15 6.00000000e-03 -8.29441686e-01 -8.29441686e-01 16 7.00000000e-03 -8.16575805e-01 -8.16575805e-01 16 8.00000000e-03 -8.04698156e-01 -8.04698156e-01 15 9.00000000e-03 -7.93626754e-01 -7.93626754e-01 16 1.00000000e-02 -7.83229199e-01 -7.83229199e-01 16 Wp results for x near 0 x W(x) (WAPR) W(x) (EXACT) Digits Correct 1.00000000e-09 9.99999999e-10 9.99999999e-10 16 2.00000000e-09 2.00000000e-09 2.00000000e-09 16 3.00000000e-09 2.99999999e-09 2.99999999e-09 16 4.00000000e-09 3.99999998e-09 3.99999998e-09 16 5.00000000e-09 4.99999998e-09 4.99999998e-09 16 6.00000000e-09 5.99999996e-09 5.99999996e-09 16 7.00000000e-09 6.99999995e-09 6.99999995e-09 16 8.00000000e-09 7.99999994e-09 7.99999994e-09 16 9.00000000e-09 8.99999992e-09 8.99999992e-09 16 1.00000000e-08 9.99999990e-09 9.99999990e-09 16 1.00000000e-02 9.90147384e-03 9.90147384e-03 16 2.00000000e-02 1.96115893e-02 1.96115893e-02 16 3.00000000e-02 2.91384592e-02 2.91384592e-02 16 4.00000000e-02 3.84896659e-02 3.84896659e-02 16 5.00000000e-02 4.76723086e-02 4.76723086e-02 16 6.00000000e-02 5.66930438e-02 5.66930438e-02 16 7.00000000e-02 6.55581227e-02 6.55581227e-02 16 8.00000000e-02 7.42734246e-02 7.42734246e-02 16 9.00000000e-02 8.28444857e-02 8.28444857e-02 16 1.00000000e-01 9.12765272e-02 9.12765272e-02 16 -1.00000000e-09 -1.00000000e-09 -1.00000000e-09 16 -2.00000000e-09 -2.00000000e-09 -2.00000000e-09 16 -3.00000000e-09 -3.00000001e-09 -3.00000001e-09 16 -4.00000000e-09 -4.00000002e-09 -4.00000002e-09 16 -5.00000000e-09 -5.00000002e-09 -5.00000003e-09 16 -6.00000000e-09 -6.00000004e-09 -6.00000004e-09 16 -7.00000000e-09 -7.00000005e-09 -7.00000005e-09 16 -8.00000000e-09 -8.00000006e-09 -8.00000006e-09 16 -9.00000000e-09 -9.00000008e-09 -9.00000008e-09 16 -1.00000000e-08 -1.00000001e-08 -1.00000001e-08 16 -1.00000000e-02 -1.01015272e-02 -1.01015272e-02 16 -2.00000000e-02 -2.04124441e-02 -2.04124441e-02 16 -3.00000000e-02 -3.09427950e-02 -3.09427950e-02 16 -4.00000000e-02 -4.17034084e-02 -4.17034084e-02 16 -5.00000000e-02 -5.27059836e-02 -5.27059836e-02 16 -6.00000000e-02 -6.39631894e-02 -6.39631894e-02 16 -7.00000000e-02 -7.54887789e-02 -7.54887789e-02 16 -8.00000000e-02 -8.72977209e-02 -8.72977209e-02 16 -9.00000000e-02 -9.94063528e-02 -9.94063528e-02 16 -1.00000000e-01 -1.11832559e-01 -1.11832559e-01 16 Other Wp results x W(x) (WAPR) W(x) (EXACT) Digits Correct 1.00000000e+01 1.74552800e+00 1.74552800e+00 16 1.00000000e+02 3.38563014e+00 3.38563014e+00 16 1.00000000e+03 5.24960285e+00 5.24960285e+00 16 1.00000000e+04 7.23184604e+00 7.23184604e+00 16 1.00000000e+05 9.28457143e+00 9.28457143e+00 16 1.00000000e+06 1.13833581e+01 1.13833581e+01 16 1.00000000e+07 1.35143440e+01 1.35143440e+01 16 1.00000000e+08 1.56689967e+01 1.56689967e+01 16 1.00000000e+09 1.78417260e+01 1.78417260e+01 16 1.00000000e+10 2.00286854e+01 2.00286854e+01 16 Wm results for x near 0 x W(x) (WAPR) W(x) (EXACT) Digits Correct 1.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 2.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 3.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 4.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 5.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 6.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 7.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 8.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 9.00000000e-40 -1.00000000e+00 -1.00000000e+00 16 1.00000000e-39 -1.00000000e+00 -1.00000000e+00 16 1.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 2.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 3.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 4.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 5.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 6.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 7.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 8.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 9.00000000e-30 -1.00000000e+00 -1.00000000e+00 16 1.00000000e-29 -1.00000000e+00 -1.00000000e+00 16 1.00000000e-20 -1.00000000e+00 -1.00000000e+00 16 2.00000000e-20 -1.00000000e+00 -1.00000000e+00 16 3.00000000e-20 -1.00000000e+00 -1.00000000e+00 16 4.00000000e-20 -1.00000000e+00 -1.00000000e+00 16 5.00000000e-20 -1.00000000e+00 -1.00000000e+00 16 6.00000000e-20 -1.00000000e+00 -1.00000000e+00 16 7.00000000e-20 -1.00000000e+00 -1.00000000e+00 16 8.00000000e-20 -1.00000000e+00 -1.00000000e+00 16 9.00000000e-20 -1.00000000e+00 -1.00000000e+00 16 1.00000000e-19 -1.00000000e+00 -1.00000000e+00 16 1.00000000e-10 -1.00002332e+00 -1.00002332e+00 16 2.00000000e-10 -1.00003297e+00 -1.00003297e+00 16 3.00000000e-10 -1.00004039e+00 -1.00004039e+00 16 4.00000000e-10 -1.00004663e+00 -1.00004663e+00 16 5.00000000e-10 -1.00005214e+00 -1.00005214e+00 16 6.00000000e-10 -1.00005711e+00 -1.00005711e+00 16 7.00000000e-10 -1.00006169e+00 -1.00006169e+00 16 8.00000000e-10 -1.00006595e+00 -1.00006595e+00 16 9.00000000e-10 -1.00006995e+00 -1.00006995e+00 16 1.00000000e-09 -1.00007373e+00 -1.00007373e+00 16 1.00000000e-05 -1.00739149e+00 -1.00739149e+00 16 2.00000000e-05 -1.01046385e+00 -1.01046385e+00 16 3.00000000e-05 -1.01282563e+00 -1.01282563e+00 16 4.00000000e-05 -1.01481959e+00 -1.01481959e+00 16 5.00000000e-05 -1.01657851e+00 -1.01657851e+00 16 6.00000000e-05 -1.01817048e+00 -1.01817048e+00 16 7.00000000e-05 -1.01963593e+00 -1.01963593e+00 16 8.00000000e-05 -1.02100123e+00 -1.02100123e+00 16 9.00000000e-05 -1.02228469e+00 -1.02228469e+00 16 1.00000000e-04 -1.02349962e+00 -1.02349962e+00 16 2.00000000e-04 -1.03334244e+00 -1.03334244e+00 15 3.00000000e-04 -1.04093919e+00 -1.04093919e+00 15 4.00000000e-04 -1.04737363e+00 -1.04737363e+00 15 5.00000000e-04 -1.05306550e+00 -1.05306550e+00 15 6.00000000e-04 -1.05823003e+00 -1.05823003e+00 15 7.00000000e-04 -1.06299509e+00 -1.06299509e+00 15 8.00000000e-04 -1.06744399e+00 -1.06744399e+00 16 9.00000000e-04 -1.07163456e+00 -1.07163456e+00 15 1.00000000e-03 -1.07560894e+00 -1.07560894e+00 15 2.00000000e-03 -1.10808188e+00 -1.10808188e+00 15 3.00000000e-03 -1.13348700e+00 -1.13348700e+00 15 4.00000000e-03 -1.15524585e+00 -1.15524585e+00 16 5.00000000e-03 -1.17468261e+00 -1.17468261e+00 15 6.00000000e-03 -1.19247585e+00 -1.19247585e+00 15 7.00000000e-03 -1.20902838e+00 -1.20902838e+00 16 8.00000000e-03 -1.22460245e+00 -1.22460245e+00 16 9.00000000e-03 -1.23938010e+00 -1.23938010e+00 15 1.00000000e-02 -1.25349379e+00 -1.25349379e+00 15 Wm results for x near -exp(-1) Offset x W(x) (WAPR) W(x) (EXACT) Digits Correct -1.00000000e-40 -9.66747560e+01 -9.66747560e+01 16 -2.00000000e-40 -9.59743374e+01 -9.59743374e+01 16 -3.00000000e-40 -9.55645938e+01 -9.55645938e+01 16 -4.00000000e-40 -9.52738649e+01 -9.52738649e+01 16 -5.00000000e-40 -9.50483515e+01 -9.50483515e+01 16 -6.00000000e-40 -9.48640895e+01 -9.48640895e+01 16 -7.00000000e-40 -9.47082952e+01 -9.47082952e+01 16 -8.00000000e-40 -9.45733378e+01 -9.45733378e+01 16 -9.00000000e-40 -9.44542952e+01 -9.44542952e+01 16 -1.00000000e-39 -9.43478067e+01 -9.43478067e+01 16 -1.00000000e-30 -7.33731103e+01 -7.33731103e+01 16 -2.00000000e-30 -7.26703389e+01 -7.26703389e+01 16 -3.00000000e-30 -7.22592002e+01 -7.22592002e+01 16 -4.00000000e-30 -7.19674727e+01 -7.19674727e+01 16 -5.00000000e-30 -7.17411798e+01 -7.17411798e+01 16 -6.00000000e-30 -7.15562776e+01 -7.15562776e+01 16 -7.00000000e-30 -7.13999397e+01 -7.13999397e+01 16 -8.00000000e-30 -7.12645097e+01 -7.12645097e+01 16 -9.00000000e-30 -7.11450489e+01 -7.11450489e+01 16 -1.00000000e-29 -7.10381852e+01 -7.10381852e+01 16 -1.00000000e-20 -4.99629843e+01 -4.99629843e+01 16 -2.00000000e-20 -4.92555773e+01 -4.92555773e+01 16 -3.00000000e-20 -4.88416735e+01 -4.88416735e+01 16 -4.00000000e-20 -4.85479597e+01 -4.85479597e+01 16 -5.00000000e-20 -4.83201118e+01 -4.83201118e+01 16 -6.00000000e-20 -4.81339297e+01 -4.81339297e+01 16 -7.00000000e-20 -4.79765031e+01 -4.79765031e+01 16 -8.00000000e-20 -4.78401250e+01 -4.78401250e+01 16 -9.00000000e-20 -4.77198242e+01 -4.77198242e+01 16 -1.00000000e-19 -4.76122059e+01 -4.76122059e+01 16 -1.00000000e-10 -2.62952388e+01 -2.62952388e+01 16 -2.00000000e-10 -2.55742914e+01 -2.55742914e+01 16 -3.00000000e-10 -2.51521833e+01 -2.51521833e+01 16 -4.00000000e-10 -2.48525155e+01 -2.48525155e+01 16 -5.00000000e-10 -2.46199710e+01 -2.46199710e+01 16 -6.00000000e-10 -2.44298992e+01 -2.44298992e+01 16 -7.00000000e-10 -2.42691466e+01 -2.42691466e+01 16 -8.00000000e-10 -2.41298594e+01 -2.41298594e+01 16 -9.00000000e-10 -2.40069706e+01 -2.40069706e+01 16 -1.00000000e-09 -2.38970196e+01 -2.38970196e+01 16 -1.00000000e-05 -1.41636008e+01 -1.41636008e+01 16 -2.00000000e-05 -1.34162445e+01 -1.34162445e+01 16 -3.00000000e-05 -1.29775328e+01 -1.29775328e+01 16 -4.00000000e-05 -1.26655140e+01 -1.26655140e+01 16 -5.00000000e-05 -1.24230404e+01 -1.24230404e+01 16 -6.00000000e-05 -1.22246178e+01 -1.22246178e+01 16 -7.00000000e-05 -1.20566300e+01 -1.20566300e+01 16 -8.00000000e-05 -1.19109413e+01 -1.19109413e+01 16 -9.00000000e-05 -1.17822992e+01 -1.17822992e+01 16 -1.00000000e-04 -1.16671145e+01 -1.16671145e+01 16 -2.00000000e-04 -1.09065574e+01 -1.09065574e+01 16 -3.00000000e-04 -1.04592111e+01 -1.04592111e+01 16 -4.00000000e-04 -1.01405924e+01 -1.01405924e+01 16 -5.00000000e-04 -9.89269952e+00 -9.89269952e+00 16 -6.00000000e-04 -9.68963797e+00 -9.68963797e+00 16 -7.00000000e-04 -9.51756976e+00 -9.51756976e+00 16 -8.00000000e-04 -9.36822217e+00 -9.36822217e+00 16 -9.00000000e-04 -9.23625197e+00 -9.23625197e+00 16 -1.00000000e-03 -9.11800647e+00 -9.11800647e+00 16 -2.00000000e-03 -8.33508138e+00 -8.33508138e+00 16 -3.00000000e-03 -7.87252138e+00 -7.87252138e+00 16 -4.00000000e-03 -7.54194042e+00 -7.54194042e+00 16 -5.00000000e-03 -7.28399714e+00 -7.28399714e+00 16 -6.00000000e-03 -7.07216205e+00 -7.07216205e+00 16 -7.00000000e-03 -6.89224149e+00 -6.89224149e+00 16 -8.00000000e-03 -6.73574166e+00 -6.73574166e+00 16 -9.00000000e-03 -6.59717173e+00 -6.59717173e+00 16 -1.00000000e-02 -6.47277512e+00 -6.47277512e+00 16 TEST02 Input X is the offset from -exp(-1). Both branches of the W function will be checked. Results for Wp(x): Offset x W(x) (WAPR) W(x) (BISECT) Digits Correct 1.00000000e-09 -9.99926269e-01 -9.99926269e-01 14 2.00000000e-09 -9.99895729e-01 -9.99895729e-01 12 3.00000000e-09 -9.99872296e-01 -9.99872296e-01 14 4.00000000e-09 -9.99852541e-01 -9.99852541e-01 13 5.00000000e-09 -9.99835137e-01 -9.99835137e-01 12 6.00000000e-09 -9.99819403e-01 -9.99819403e-01 12 7.00000000e-09 -9.99804933e-01 -9.99804933e-01 13 8.00000000e-09 -9.99791466e-01 -9.99791466e-01 12 9.00000000e-09 -9.99778817e-01 -9.99778817e-01 14 1.00000000e-08 -9.99766854e-01 -9.99766854e-01 13 1.10000000e-08 -9.99755475e-01 -9.99755475e-01 13 Results for Wm(x): Offset x W(x) (WAPR) W(x) (BISECT) Digits Correct 1.00000000e-09 -1.00007373e+00 -1.00007373e+00 12 2.00000000e-09 -1.00010428e+00 -1.00010428e+00 13 3.00000000e-09 -1.00012771e+00 -1.00012771e+00 12 4.00000000e-09 -1.00014747e+00 -1.00014747e+00 13 5.00000000e-09 -1.00016488e+00 -1.00016488e+00 12 6.00000000e-09 -1.00018062e+00 -1.00018062e+00 12 7.00000000e-09 -1.00019509e+00 -1.00019509e+00 12 8.00000000e-09 -1.00020856e+00 -1.00020856e+00 12 9.00000000e-09 -1.00022122e+00 -1.00022122e+00 12 1.00000000e-08 -1.00023318e+00 -1.00023318e+00 12 1.10000000e-08 -1.00024456e+00 -1.00024456e+00 13 TEST03 Input X is the argument. Wp has been selected (maximum x is > 0) Results for Wp(x): x W(x) (WAPR) W(x) (BISECT) Digits Correct 0.00000000e+00 0.00000000e+00 0.00000000e+00 16 5.00000000e+18 3.93826444e+01 3.93826444e+01 16 1.00000000e+19 4.00587692e+01 4.00587692e+01 16 1.50000000e+19 4.04544063e+01 4.04544063e+01 16 2.00000000e+19 4.07351721e+01 4.07351721e+01 16 2.50000000e+19 4.09529828e+01 4.09529828e+01 16 3.00000000e+19 4.11309677e+01 4.11309677e+01 16 3.50000000e+19 4.12814661e+01 4.12814661e+01 16 4.00000000e+19 4.14118442e+01 4.14118442e+01 16 4.50000000e+19 4.15268539e+01 4.15268539e+01 16 5.00000000e+19 4.16297399e+01 4.16297399e+01 16 5.50000000e+19 4.17228167e+01 4.17228167e+01 16 6.00000000e+19 4.18077935e+01 4.18077935e+01 16 6.50000000e+19 4.18859681e+01 4.18859681e+01 16 7.00000000e+19 4.19583495e+01 4.19583495e+01 16 7.50000000e+19 4.20257376e+01 4.20257376e+01 16 8.00000000e+19 4.20887772e+01 4.20887772e+01 16 8.50000000e+19 4.21479958e+01 4.21479958e+01 16 9.00000000e+19 4.22038304e+01 4.22038304e+01 16 9.50000000e+19 4.22566469e+01 4.22566469e+01 16 1.00000000e+20 4.23067551e+01 4.23067551e+01 16 TOMS743_TEST Normal end of execution. 27 May 2021 03:02:01 PM