| Home | Libraries | People | FAQ | More

### Tables of Error Rates for all Functions

Table 23.3. Error rates for beta

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Beta Function: Small Values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = +INFε (Mean = +INFε) And other failures.) (Rmath 3.2.3: Max = 1.14ε (Mean = 0.574ε)) | Max = 2.86ε (Mean = 1.22ε) (<cmath>: Max = 364ε (Mean = 76.2ε)) | Max = 2.86ε (Mean = 1.22ε) | Max = 2.23ε (Mean = 1.14ε)
Beta Function: Medium Values | Max = 0.978ε (Mean = 0.0595ε) (GSL 2.1: Max = 1.18e+03ε (Mean = 238ε)) (Rmath 3.2.3: Max = 1.09e+03ε (Mean = 265ε)) | Max = 61.4ε (Mean = 19.4ε) (<cmath>: Max = 1.07e+03ε (Mean = 264ε)) | Max = 107ε (Mean = 24.5ε) | Max = 96.5ε (Mean = 22.4ε)
Beta Function: Divergent Values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 12.1ε (Mean = 1.99ε)) (Rmath 3.2.3: Max = 176ε (Mean = 28ε)) | Max = 8.99ε (Mean = 2.44ε) (<cmath>: Max = 128ε (Mean = 23.8ε)) | Max = 18.8ε (Mean = 2.71ε) | Max = 11.4ε (Mean = 2.19ε)

Table 23.4. Error rates for beta (incomplete)

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Incomplete Beta Function: Small Values | Max = 0ε (Mean = 0ε) | Max = 11.1ε (Mean = 2.32ε) | Max = 18.7ε (Mean = 3.19ε) | Max = 9.94ε (Mean = 2.17ε)
Incomplete Beta Function: Medium Values | Max = 0.568ε (Mean = 0.0254ε) | Max = 69.2ε (Mean = 13.4ε) | Max = 174ε (Mean = 25ε) | Max = 90ε (Mean = 12.7ε)
Incomplete Beta Function: Large and Diverse Values | Max = 0.999ε (Mean = 0.0325ε) | Max = 6.84e+04ε (Mean = 2.76e+03ε) | Max = 6.86e+04ε (Mean = 2.79e+03ε) | Max = 633ε (Mean = 29.7ε)
Incomplete Beta Function: Small Integer Values | Max = 0.786ε (Mean = 0.0323ε) | Max = 11.6ε (Mean = 3.6ε) | Max = 51.8ε (Mean = 11ε) | Max = 26ε (Mean = 6.28ε)

Table 23.5. Error rates for betac

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Incomplete Beta Function: Small Values | Max = 0.676ε (Mean = 0.0302ε) | Max = 9.92ε (Mean = 2.3ε) | Max = 11.2ε (Mean = 2.94ε) | Max = 8.94ε (Mean = 2.06ε)
Incomplete Beta Function: Medium Values | Max = 0.949ε (Mean = 0.098ε) | Max = 63.5ε (Mean = 13.5ε) | Max = 97.6ε (Mean = 24.3ε) | Max = 90.6ε (Mean = 14.8ε)
Incomplete Beta Function: Large and Diverse Values | Max = 1.12ε (Mean = 0.0458ε) | Max = 1.05e+05ε (Mean = 5.45e+03ε) | Max = 1.04e+05ε (Mean = 5.46e+03ε) | Max = 3.72e+03ε (Mean = 113ε)
Incomplete Beta Function: Small Integer Values | Max = 0.586ε (Mean = 0.0314ε) | Max = 11.1ε (Mean = 3.65ε) | Max = 103ε (Mean = 17.4ε) | Max = 26.2ε (Mean = 6.36ε)

Table 23.6. Error rates for binomial_coefficient

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Binomials: small arguments | Max = 1ε (Mean = 0.369ε) | Max = 1.5ε (Mean = 0.339ε) | Max = 1.5ε (Mean = 0.339ε) | Max = 1ε (Mean = 0.369ε)
Binomials: large arguments | Max = 0.939ε (Mean = 0.314ε) | Max = 26.6ε (Mean = 6.13ε) | Max = 53.2ε (Mean = 10.8ε) | Max = 37.2ε (Mean = 7.4ε)

Table 23.7. Error rates for boost::math::powm1

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
powm1 | (<math.h>: Max = 2.04ε (Mean = 0.493ε)) Max = 2.04ε (Mean = 0.493ε) (<cmath>: Max = 2.04ε (Mean = 0.493ε)) | (Rmath 3.2.3: Max = 1.06ε (Mean = 0.425ε)) (GSL 2.1: Max = 1.06ε (Mean = 0.425ε)) Max = 1.06ε (Mean = 0.425ε) | Max = 1.88ε (Mean = 0.49ε) (<math.h>: Max = 1.88ε (Mean = 0.49ε)) | (<math.h>: Max = 1.84ε (Mean = 0.486ε)) Max = 1.84ε (Mean = 0.486ε)

Table 23.8. Error rates for cbrt

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
cbrt Function | Max = 0ε (Mean = 0ε) | Max = 1.34ε (Mean = 0.471ε) (<cmath>: Max = 1.34ε (Mean = 0.471ε)) (<math.h>: Max = 1.34ε (Mean = 0.471ε)) | Max = 1.34ε (Mean = 0.471ε) (<math.h>: Max = 1.34ε (Mean = 0.471ε)) | Max = 1.7ε (Mean = 0.565ε) (<math.h>: Max = 1.7ε (Mean = 0.565ε))

Table 23.9. Error rates for cos_pi

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
sin_pi and cos_pi | Max = 0ε (Mean = 0ε) | Max = 0.991ε (Mean = 0.302ε) | Max = 0.991ε (Mean = 0.302ε) | Max = 0.996ε (Mean = 0.284ε)
sin_pi and cos_pi near integers and half integers | Max = 0ε (Mean = 0ε) | Max = 0.976ε (Mean = 0.28ε) | Max = 0.976ε (Mean = 0.28ε) | Max = 0.996ε (Mean = 0.298ε)

Table 23.10. Error rates for cyl_bessel_i

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel I0: Mathworld Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 270ε (Mean = 91.6ε) And other failures.) (Rmath 3.2.3: Max = 1.52ε (Mean = 0.622ε) And other failures.) | Max = 1.95ε (Mean = 0.738ε) (<cmath>: Max = 8.49ε (Mean = 3.46ε) And other failures.) | Max = 1.95ε (Mean = 0.661ε) | Max = 0.762ε (Mean = 0.329ε)
Bessel I1: Mathworld Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 128ε (Mean = 41ε) And other failures.) (Rmath 3.2.3: Max = 1.53ε (Mean = 0.483ε) And other failures.) | Max = 0.64ε (Mean = 0.202ε) (<cmath>: Max = 5ε (Mean = 2.15ε) And other failures.) | Max = 0.64ε (Mean = 0.202ε) | Max = 0.767ε (Mean = 0.398ε)
Bessel In: Mathworld Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 2.31ε (Mean = 0.838ε) And other failures.) (Rmath 3.2.3: Max = 1.73ε (Mean = 0.601ε) And other failures.) | Max = 1.8ε (Mean = 1.33ε) (<cmath>: Max = 430ε (Mean = 163ε) And other failures.) | Max = 463ε (Mean = 140ε) | Max = 3.46ε (Mean = 1.32ε)
Bessel Iv: Mathworld Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 5.95ε (Mean = 2.08ε) And other failures.) (Rmath 3.2.3: Max = 3.53ε (Mean = 1.39ε)) | Max = 4.12ε (Mean = 1.85ε) (<cmath>: Max = 616ε (Mean = 221ε) And other failures.) | Max = 4.12ε (Mean = 1.95ε) | Max = 2.97ε (Mean = 1.24ε)
Bessel In: Random Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 261ε (Mean = 53.2ε) And other failures.) (Rmath 3.2.3: Max = 7.37ε (Mean = 2.4ε)) | Max = 4.62ε (Mean = 1.06ε) (<cmath>: Max = 645ε (Mean = 132ε)) | Max = 176ε (Mean = 39.1ε) | Max = 9.67ε (Mean = 1.88ε)
Bessel Iv: Random Data | Max = 0.661ε (Mean = 0.0441ε) (GSL 2.1: Max = 6.18e+03ε (Mean = 1.55e+03ε) And other failures.) (Rmath 3.2.3: Max = 4.28e+08ε (Mean = 2.85e+07ε)) | Max = 8.35ε (Mean = 1.62ε) (<cmath>: Max = 1.05e+03ε (Mean = 224ε) And other failures.) | Max = 283ε (Mean = 88.4ε) | Max = 7.46ε (Mean = 1.71ε)
Bessel Iv: Mathworld Data (large values) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 37ε (Mean = 18ε) And other failures.) (Rmath 3.2.3: Max = 3.77e+168ε (Mean = 2.39e+168ε) And other failures.) | Max = 14.7ε (Mean = 6.66ε) (<cmath>: Max = 118ε (Mean = 57.2ε) And other failures.) | Max = 14.7ε (Mean = 6.59ε) | Max = 3.67ε (Mean = 1.64ε)

Table 23.11. Error rates for cyl_bessel_i (integer orders)

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel I0: Mathworld Data (Integer Version) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0.79ε (Mean = 0.482ε)) (Rmath 3.2.3: Max = 1.52ε (Mean = 0.622ε) And other failures.) | Max = 1.95ε (Mean = 0.738ε) (<cmath>: Max = 8.49ε (Mean = 3.46ε) And other failures.) | Max = 1.95ε (Mean = 0.661ε) | Max = 0.762ε (Mean = 0.329ε)
Bessel I1: Mathworld Data (Integer Version) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0.82ε (Mean = 0.456ε)) (Rmath 3.2.3: Max = 1.53ε (Mean = 0.483ε) And other failures.) | Max = 0.64ε (Mean = 0.202ε) (<cmath>: Max = 5ε (Mean = 2.15ε) And other failures.) | Max = 0.64ε (Mean = 0.202ε) | Max = 0.767ε (Mean = 0.398ε)
Bessel In: Mathworld Data (Integer Version) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 5.15ε (Mean = 2.13ε) And other failures.) (Rmath 3.2.3: Max = 1.73ε (Mean = 0.601ε) And other failures.) | Max = 1.8ε (Mean = 1.33ε) (<cmath>: Max = 430ε (Mean = 163ε) And other failures.) | Max = 463ε (Mean = 140ε) | Max = 3.46ε (Mean = 1.32ε)

Table 23.12. Error rates for cyl_bessel_i_prime

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel I'0: Mathworld Data | Max = 0ε (Mean = 0ε) | Max = 0.82ε (Mean = 0.259ε) | Max = 0.82ε (Mean = 0.259ε) | Max = 0.82ε (Mean = 0.354ε)
Bessel I'1: Mathworld Data | Max = 0ε (Mean = 0ε) | Max = 1.97ε (Mean = 0.757ε) | Max = 1.97ε (Mean = 0.757ε) | Max = 1.36ε (Mean = 0.782ε)
Bessel I'n: Mathworld Data | Max = 0ε (Mean = 0ε) | Max = 2.31ε (Mean = 1.41ε) | Max = 701ε (Mean = 212ε) | Max = 3.61ε (Mean = 1.22ε)
Bessel I'v: Mathworld Data | Max = 1.62ε (Mean = 0.512ε) | Max = 2.89e+03ε (Mean = 914ε) | Max = 2.89e+03ε (Mean = 914ε) | Max = 3.76e+03ε (Mean = 1.19e+03ε)
Bessel I'n: Random Data | Max = 0ε (Mean = 0ε) | Max = 3.95ε (Mean = 1.06ε) | Max = 195ε (Mean = 37.1ε) | Max = 9.85ε (Mean = 1.82ε)
Bessel I'v: Random Data | Max = 0ε (Mean = 0ε) | Max = 14.1ε (Mean = 2.93ε) | Max = 336ε (Mean = 68.5ε) | Max = 14ε (Mean = 2.5ε)
Bessel I'v: Mathworld Data (large values) | Max = 0ε (Mean = 0ε) | Max = 42.6ε (Mean = 20.2ε) | Max = 42.6ε (Mean = 20.2ε) | Max = 59.5ε (Mean = 26.6ε)

Table 23.13. Error rates for cyl_bessel_i_prime (integer orders)

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel I'0: Mathworld Data (Integer Version) | Max = 0ε (Mean = 0ε) | Max = 0.82ε (Mean = 0.259ε) | Max = 0.82ε (Mean = 0.259ε) | Max = 0.82ε (Mean = 0.354ε)
Bessel I'1: Mathworld Data (Integer Version) | Max = 0ε (Mean = 0ε) | Max = 1.97ε (Mean = 0.757ε) | Max = 1.97ε (Mean = 0.757ε) | Max = 1.36ε (Mean = 0.782ε)
Bessel I'n: Mathworld Data (Integer Version) | Max = 0ε (Mean = 0ε) | Max = 2.31ε (Mean = 1.41ε) | Max = 701ε (Mean = 212ε) | Max = 3.61ε (Mean = 1.22ε)

Table 23.14. Error rates for cyl_bessel_j

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel J0: Mathworld Data | Max = 6.55ε (Mean = 2.86ε) (<cmath>: Max = 5.04ε (Mean = 1.78ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0.629ε (Mean = 0.223ε) And other failures.) (Rmath 3.2.3: Max = 0.629ε (Mean = 0.223ε) And other failures.) | Max = 6.55ε (Mean = 2.86ε) | Max = 2.52ε (Mean = 1.2ε)
Bessel J0: Mathworld Data (Tricky cases) | Max = 1.64e+08ε (Mean = 6.69e+07ε) (<cmath>: Max = 4.79e+08ε (Mean = 1.96e+08ε)) | Max = 8e+04ε (Mean = 3.27e+04ε) (GSL 2.1: Max = 6.5e+07ε (Mean = 2.66e+07ε)) (Rmath 3.2.3: Max = 1.04e+07ε (Mean = 4.29e+06ε)) | Max = 1.64e+08ε (Mean = 6.69e+07ε) | Max = 1e+07ε (Mean = 4.09e+06ε)
Bessel J1: Mathworld Data | Max = 3.59ε (Mean = 1.33ε) (<cmath>: Max = 6.1ε (Mean = 2.95ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 6.62ε (Mean = 2.35ε) And other failures.) (Rmath 3.2.3: Max = 0.946ε (Mean = 0.39ε) And other failures.) | Max = 1.44ε (Mean = 0.637ε) | Max = 1.73ε (Mean = 0.976ε)
Bessel J1: Mathworld Data (tricky cases) | Max = 2.18e+05ε (Mean = 9.76e+04ε) (<cmath>: Max = 2.15e+06ε (Mean = 1.58e+06ε)) | Max = 106ε (Mean = 47.5ε) (GSL 2.1: Max = 8.75e+05ε (Mean = 5.32e+05ε)) (Rmath 3.2.3: Max = 2.93e+06ε (Mean = 1.7e+06ε)) | Max = 2.18e+05ε (Mean = 9.76e+04ε) | Max = 3.23e+04ε (Mean = 1.45e+04ε)
Bessel JN: Mathworld Data | Max = 6.85ε (Mean = 3.35ε) (<cmath>: Max = 2.13e+19ε (Mean = 5.16e+18ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 6.9e+05ε (Mean = 2.15e+05ε) And other failures.) (Rmath 3.2.3: Max = +INFε (Mean = +INFε) And other failures.) | Max = 463ε (Mean = 112ε) | Max = 14.7ε (Mean = 5.4ε)
Bessel J: Mathworld Data | Max = 14.7ε (Mean = 4.11ε) (<cmath>: Max = 3.49e+05ε (Mean = 8.09e+04ε) And other failures.) | Max = 10ε (Mean = 2.24ε) (GSL 2.1: Max = 2.39e+05ε (Mean = 5.37e+04ε) And other failures.) (Rmath 3.2.3: Max = +INFε (Mean = +INFε) And other failures.) | Max = 14.7ε (Mean = 4.22ε) | Max = 14.9ε (Mean = 3.89ε)
Bessel J: Mathworld Data (large values) | Max = 607ε (Mean = 305ε) (<cmath>: Max = 34.9ε (Mean = 17.4ε) And other failures.) | Max = 0.536ε (Mean = 0.268ε) (GSL 2.1: Max = 4.91e+03ε (Mean = 2.46e+03ε) And other failures.) (Rmath 3.2.3: Max = 5.9ε (Mean = 3.76ε)) | Max = 607ε (Mean = 305ε) | Max = 9.31ε (Mean = 5.52ε)
Bessel JN: Random Data | Max = 50.8ε (Mean = 3.69ε) (<cmath>: Max = 1.12e+03ε (Mean = 88.7ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 75.7ε (Mean = 5.36ε)) (Rmath 3.2.3: Max = 3.93ε (Mean = 1.22ε)) | Max = 99.6ε (Mean = 22ε) | Max = 17.5ε (Mean = 1.46ε)
Bessel J: Random Data | Max = 11.4ε (Mean = 1.68ε) (<cmath>: Max = 501ε (Mean = 52.3ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 15.5ε (Mean = 3.33ε) And other failures.) (Rmath 3.2.3: Max = 6.74ε (Mean = 1.3ε)) | Max = 260ε (Mean = 34ε) | Max = 9.24ε (Mean = 1.17ε)
Bessel J: Random Data (Tricky large values) | Max = 785ε (Mean = 94.2ε) (<cmath>: Max = 5.01e+17ε (Mean = 6.23e+16ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 2.48e+05ε (Mean = 5.11e+04ε)) (Rmath 3.2.3: Max = 71.6ε (Mean = 11.7ε)) | Max = 785ε (Mean = 97.4ε) | Max = 59.2ε (Mean = 8.67ε)

Table 23.15. Error rates for cyl_bessel_j (integer orders)

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel J0: Mathworld Data (Integer Version) | Max = 6.55ε (Mean = 2.86ε) (<cmath>: Max = 5.04ε (Mean = 1.78ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.12ε (Mean = 0.488ε)) (Rmath 3.2.3: Max = 0.629ε (Mean = 0.223ε) And other failures.) | Max = 6.55ε (Mean = 2.86ε) | Max = 2.52ε (Mean = 1.2ε) (<math.h>: Max = 1.89ε (Mean = 0.988ε))
Bessel J0: Mathworld Data (Tricky cases) (Integer Version) | Max = 1.64e+08ε (Mean = 6.69e+07ε) (<cmath>: Max = 4.79e+08ε (Mean = 1.96e+08ε)) | Max = 8e+04ε (Mean = 3.27e+04ε) (GSL 2.1: Max = 1e+07ε (Mean = 4.11e+06ε)) (Rmath 3.2.3: Max = 1.04e+07ε (Mean = 4.29e+06ε)) | Max = 1.64e+08ε (Mean = 6.69e+07ε) | Max = 1e+07ε (Mean = 4.09e+06ε) (<math.h>: Max = 2.54e+08ε (Mean = 1.04e+08ε))
Bessel J1: Mathworld Data (Integer Version) | Max = 3.59ε (Mean = 1.33ε) (<cmath>: Max = 6.1ε (Mean = 2.95ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.89ε (Mean = 0.721ε)) (Rmath 3.2.3: Max = 0.946ε (Mean = 0.39ε) And other failures.) | Max = 1.44ε (Mean = 0.637ε) | Max = 1.73ε (Mean = 0.976ε) (<math.h>: Max = 11.4ε (Mean = 4.15ε))
Bessel J1: Mathworld Data (tricky cases) (Integer Version) | Max = 2.18e+05ε (Mean = 9.76e+04ε) (<cmath>: Max = 2.15e+06ε (Mean = 1.58e+06ε)) | Max = 106ε (Mean = 47.5ε) (GSL 2.1: Max = 1.26e+06ε (Mean = 6.28e+05ε)) (Rmath 3.2.3: Max = 2.93e+06ε (Mean = 1.7e+06ε)) | Max = 2.18e+05ε (Mean = 9.76e+04ε) | Max = 3.23e+04ε (Mean = 1.45e+04ε) (<math.h>: Max = 1.44e+07ε (Mean = 6.5e+06ε))
Bessel JN: Mathworld Data (Integer Version) | Max = 6.85ε (Mean = 3.35ε) (<cmath>: Max = 2.13e+19ε (Mean = 5.16e+18ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 6.9e+05ε (Mean = 2.53e+05ε) And other failures.) (Rmath 3.2.3: Max = +INFε (Mean = +INFε) And other failures.) | Max = 463ε (Mean = 112ε) | Max = 14.7ε (Mean = 5.4ε) (<math.h>: Max = +INFε (Mean = +INFε) And other failures.)

Table 23.16. Error rates for cyl_bessel_j_prime

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel J0': Mathworld Data | Max = 0ε (Mean = 0ε) | Max = 18.9ε (Mean = 6.82ε) | Max = 18.9ε (Mean = 6.72ε) | Max = 6.62ε (Mean = 2.55ε)
Bessel J0': Mathworld Data (Tricky cases) | Max = 0ε (Mean = 0ε) | Max = 7.44ε (Mean = 3.34ε) | Max = 7.44ε (Mean = 3.31ε) | Max = 3.67ε (Mean = 1.74ε)
Bessel J1': Mathworld Data | Max = 0ε (Mean = 0ε) | Max = 7.9ε (Mean = 3.37ε) | Max = 7.9ε (Mean = 3.37ε) | Max = 0.999ε (Mean = 0.627ε)
Bessel J1': Mathworld Data (tricky cases) | Max = 287ε (Mean = 129ε) | Max = 5.88e+05ε (Mean = 2.63e+05ε) | Max = 5.88e+05ε (Mean = 2.63e+05ε) | Max = 288ε (Mean = 129ε)
Bessel JN': Mathworld Data | Max = 0.527ε (Mean = 0.128ε) | Max = 1.29e+03ε (Mean = 312ε) | Max = 1.29e+03ε (Mean = 355ε) | Max = 14ε (Mean = 6.13ε)
Bessel J': Mathworld Data | Max = 21.5ε (Mean = 4.7ε) | Max = 42.5ε (Mean = 9.31ε) | Max = 42.5ε (Mean = 9.32ε) | Max = 23.7ε (Mean = 8ε)
Bessel J': Mathworld Data (large values) | Max = 0ε (Mean = 0ε) | Max = 989ε (Mean = 495ε) | Max = 989ε (Mean = 495ε) | Max = 2.9ε (Mean = 1.61ε)
Bessel JN': Random Data | Max = 0.593ε (Mean = 0.0396ε) | Max = 11.3ε (Mean = 1.85ε) | Max = 79.4ε (Mean = 16.2ε) | Max = 6.34ε (Mean = 0.999ε)
Bessel J': Random Data | Max = 0.885ε (Mean = 0.033ε) | Max = 139ε (Mean = 6.47ε) | Max = 279ε (Mean = 27.2ε) | Max = 176ε (Mean = 9.75ε)
Bessel J': Random Data (Tricky large values) | Max = 0ε (Mean = 0ε) | Max = 474ε (Mean = 62.2ε) | Max = 474ε (Mean = 64.5ε) | Max = 379ε (Mean = 45.4ε)

Table 23.17. Error rates for cyl_bessel_j_prime (integer orders)

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel J0': Mathworld Data (Integer Version) | Max = 0ε (Mean = 0ε) | Max = 18.9ε (Mean = 6.82ε) | Max = 18.9ε (Mean = 6.72ε) | Max = 6.62ε (Mean = 2.55ε)
Bessel J0': Mathworld Data (Tricky cases) (Integer Version) | Max = 0ε (Mean = 0ε) | Max = 7.44ε (Mean = 3.34ε) | Max = 7.44ε (Mean = 3.31ε) | Max = 3.67ε (Mean = 1.74ε)
Bessel J1': Mathworld Data (Integer Version) | Max = 0ε (Mean = 0ε) | Max = 7.9ε (Mean = 3.37ε) | Max = 7.9ε (Mean = 3.37ε) | Max = 0.999ε (Mean = 0.627ε)
Bessel J1': Mathworld Data (tricky cases) (Integer Version) | Max = 287ε (Mean = 129ε) | Max = 5.88e+05ε (Mean = 2.63e+05ε) | Max = 5.88e+05ε (Mean = 2.63e+05ε) | Max = 288ε (Mean = 129ε)
Bessel JN': Mathworld Data (Integer Version) | Max = 0.527ε (Mean = 0.128ε) | Max = 1.29e+03ε (Mean = 312ε) | Max = 1.29e+03ε (Mean = 355ε) | Max = 14ε (Mean = 6.13ε)

Table 23.18. Error rates for cyl_bessel_k

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel K0: Mathworld Data | Max = 0.833ε (Mean = 0.436ε) (<cmath>: Max = 9.33ε (Mean = 3.25ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 6.04ε (Mean = 2.16ε)) (Rmath 3.2.3: Max = 0.833ε (Mean = 0.601ε)) | Max = 0.833ε (Mean = 0.436ε) | Max = 0.833ε (Mean = 0.552ε)
Bessel K1: Mathworld Data | Max = 0.786ε (Mean = 0.329ε) (<cmath>: Max = 8.94ε (Mean = 3.19ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 6.26ε (Mean = 2.21ε)) (Rmath 3.2.3: Max = 0.894ε (Mean = 0.516ε)) | Max = 0.786ε (Mean = 0.329ε) | Max = 0.786ε (Mean = 0.39ε)
Bessel Kn: Mathworld Data | Max = 2.6ε (Mean = 1.21ε) (<cmath>: Max = 12.9ε (Mean = 4.91ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 3.36ε (Mean = 1.43ε) And other failures.) (Rmath 3.2.3: Max = 8.48ε (Mean = 2.98ε)) | Max = 2.6ε (Mean = 1.21ε) | Max = 3.63ε (Mean = 1.46ε)
Bessel Kv: Mathworld Data | Max = 3.58ε (Mean = 2.39ε) (<cmath>: Max = 13ε (Mean = 4.81ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 5.47ε (Mean = 2.04ε) And other failures.) (Rmath 3.2.3: Max = 3.15ε (Mean = 1.35ε)) | Max = 5.21ε (Mean = 2.53ε) | Max = 4.78ε (Mean = 2.19ε)
Bessel Kv: Mathworld Data (large values) | Max = 42.3ε (Mean = 21ε) (<cmath>: Max = 42.3ε (Mean = 19.8ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 308ε (Mean = 142ε) And other failures.) (Rmath 3.2.3: Max = 84.6ε (Mean = 37.8ε)) | Max = 42.3ε (Mean = 21ε) | Max = 59.8ε (Mean = 26.9ε)
Bessel Kn: Random Data | Max = 4.55ε (Mean = 1.12ε) (<cmath>: Max = 13.9ε (Mean = 2.91ε)) | Max = 0.764ε (Mean = 0.0348ε) (GSL 2.1: Max = 8.71ε (Mean = 1.76ε) And other failures.) (Rmath 3.2.3: Max = 7.47ε (Mean = 1.34ε)) | Max = 4.55ε (Mean = 1.12ε) | Max = 9.34ε (Mean = 1.7ε)
Bessel Kv: Random Data | Max = 7.88ε (Mean = 1.48ε) (<cmath>: Max = 13.6ε (Mean = 2.68ε) And other failures.) | Max = 0.507ε (Mean = 0.0313ε) (GSL 2.1: Max = 9.71ε (Mean = 1.47ε) And other failures.) (Rmath 3.2.3: Max = 7.37ε (Mean = 1.49ε)) | Max = 7.88ε (Mean = 1.47ε) | Max = 8.33ε (Mean = 1.62ε)

Table 23.19. Error rates for cyl_bessel_k (integer orders)

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel K0: Mathworld Data (Integer Version) | Max = 0.833ε (Mean = 0.436ε) (<cmath>: Max = 9.33ε (Mean = 3.25ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.2ε (Mean = 0.733ε)) (Rmath 3.2.3: Max = 0.833ε (Mean = 0.601ε)) | Max = 0.833ε (Mean = 0.436ε) | Max = 0.833ε (Mean = 0.552ε)
Bessel K1: Mathworld Data (Integer Version) | Max = 0.786ε (Mean = 0.329ε) (<cmath>: Max = 8.94ε (Mean = 3.19ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0.626ε (Mean = 0.333ε)) (Rmath 3.2.3: Max = 0.894ε (Mean = 0.516ε)) | Max = 0.786ε (Mean = 0.329ε) | Max = 0.786ε (Mean = 0.39ε)
Bessel Kn: Mathworld Data (Integer Version) | Max = 2.6ε (Mean = 1.21ε) (<cmath>: Max = 12.9ε (Mean = 4.91ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 168ε (Mean = 59.5ε)) (Rmath 3.2.3: Max = 8.48ε (Mean = 2.98ε)) | Max = 2.6ε (Mean = 1.21ε) | Max = 3.63ε (Mean = 1.46ε)

Table 23.20. Error rates for cyl_bessel_k_prime

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel K'0: Mathworld Data | Max = 0ε (Mean = 0ε) | Max = 0.786ε (Mean = 0.329ε) | Max = 0.786ε (Mean = 0.329ε) | Max = 0.786ε (Mean = 0.39ε)
Bessel K'1: Mathworld Data | Max = 0ε (Mean = 0ε) | Max = 0.736ε (Mean = 0.389ε) | Max = 0.736ε (Mean = 0.389ε) | Max = 0.761ε (Mean = 0.444ε)
Bessel K'n: Mathworld Data | Max = 0ε (Mean = 0ε) | Max = 2.16ε (Mean = 1.08ε) | Max = 2.16ε (Mean = 1.08ε) | Max = 4.17ε (Mean = 1.75ε)
Bessel K'v: Mathworld Data | Max = 0ε (Mean = 0ε) | Max = 3.94ε (Mean = 2.44ε) | Max = 3.94ε (Mean = 2.34ε) | Max = 3.94ε (Mean = 1.47ε)
Bessel K'v: Mathworld Data (large values) | Max = 0ε (Mean = 0ε) | Max = 59.2ε (Mean = 42.9ε) | Max = 58.7ε (Mean = 42.6ε) | Max = 18.6ε (Mean = 11.8ε)
Bessel K'n: Random Data | Max = 0ε (Mean = 0ε) | Max = 4.45ε (Mean = 1.19ε) | Max = 4.45ε (Mean = 1.19ε) | Max = 9.67ε (Mean = 1.73ε)
Bessel K'v: Random Data | Max = 0ε (Mean = 0ε) | Max = 7.95ε (Mean = 1.53ε) | Max = 7.95ε (Mean = 1.52ε) | Max = 8.32ε (Mean = 1.65ε)

Table 23.21. Error rates for cyl_bessel_k_prime (integer orders)

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel K'0: Mathworld Data (Integer Version) | Max = 0ε (Mean = 0ε) | Max = 0.786ε (Mean = 0.329ε) | Max = 0.786ε (Mean = 0.329ε) | Max = 0.786ε (Mean = 0.39ε)
Bessel K'1: Mathworld Data (Integer Version) | Max = 0ε (Mean = 0ε) | Max = 0.736ε (Mean = 0.389ε) | Max = 0.736ε (Mean = 0.389ε) | Max = 0.761ε (Mean = 0.444ε)
Bessel K'n: Mathworld Data (Integer Version) | Max = 0ε (Mean = 0ε) | Max = 2.16ε (Mean = 1.08ε) | Max = 2.16ε (Mean = 1.08ε) | Max = 4.17ε (Mean = 1.75ε)

Table 23.22. Error rates for cyl_neumann

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Y0: Mathworld Data | Max = 5.53ε (Mean = 2.4ε) (<cmath>: Max = 2.05e+05ε (Mean = 6.87e+04ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 60.9ε (Mean = 20.4ε)) (Rmath 3.2.3: Max = 167ε (Mean = 56.5ε)) | Max = 5.53ε (Mean = 2.4ε) | Max = 4.61ε (Mean = 2.29ε)
Y1: Mathworld Data | Max = 6.33ε (Mean = 2.25ε) (<cmath>: Max = 9.71e+03ε (Mean = 4.08e+03ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 23.4ε (Mean = 8.1ε)) (Rmath 3.2.3: Max = 193ε (Mean = 64.4ε)) | Max = 6.33ε (Mean = 2.29ε) | Max = 4.75ε (Mean = 1.72ε)
Yn: Mathworld Data | Max = 55.2ε (Mean = 17.8ε) (<cmath>: Max = 2.2e+20ε (Mean = 6.97e+19ε) And other failures.) | Max = 0.993ε (Mean = 0.314ε) (GSL 2.1: Max = 2.41e+05ε (Mean = 7.62e+04ε) And other failures.) (Rmath 3.2.3: Max = 1.24e+04ε (Mean = 4e+03ε)) | Max = 55.2ε (Mean = 17.8ε) | Max = 35ε (Mean = 11.9ε)
Yv: Mathworld Data | Max = 10.7ε (Mean = 4.93ε) (<cmath>: Max = 3.49e+15ε (Mean = 1.05e+15ε) And other failures.) | Max = 10ε (Mean = 3.02ε) (GSL 2.1: Max = 1.07e+05ε (Mean = 3.22e+04ε) And other failures.) (Rmath 3.2.3: Max = 243ε (Mean = 73.9ε)) | Max = 10.7ε (Mean = 5.1ε) | Max = 7.89ε (Mean = 3.27ε)
Yv: Mathworld Data (large values) | Max = 1.7ε (Mean = 1.33ε) (<cmath>: Max = 43.2ε (Mean = 16.3ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 60.8ε (Mean = 23ε) And other failures.) (Rmath 3.2.3: Max = 0.682ε (Mean = 0.335ε)) | Max = 1.7ε (Mean = 1.33ε) | Max = 0.682ε (Mean = 0.423ε)
Y0 and Y1: Random Data | Max = 10.8ε (Mean = 3.04ε) (<cmath>: Max = 2.59e+03ε (Mean = 500ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 34.4ε (Mean = 8.9ε)) (Rmath 3.2.3: Max = 83ε (Mean = 14.2ε)) | Max = 10.8ε (Mean = 3.04ε) | Max = 4.17ε (Mean = 1.24ε)
Yn: Random Data | Max = 338ε (Mean = 27.5ε) (<cmath>: Max = 4.01e+03ε (Mean = 348ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 500ε (Mean = 47.8ε)) (Rmath 3.2.3: Max = 691ε (Mean = 67.9ε)) | Max = 338ε (Mean = 27.5ε) | Max = 117ε (Mean = 10.2ε)
Yv: Random Data | Max = 2.08e+03ε (Mean = 149ε) (<cmath>: Max = +INFε (Mean = +INFε) And other failures.) | Max = 1.53ε (Mean = 0.102ε) (GSL 2.1: Max = 1.41e+06ε (Mean = 7.67e+04ε)) (Rmath 3.2.3: Max = 1.79e+05ε (Mean = 9.64e+03ε)) | Max = 2.08e+03ε (Mean = 149ε) | Max = 1.23e+03ε (Mean = 69.9ε)

Table 23.23. Error rates for cyl_neumann (integer orders)

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Y0: Mathworld Data (Integer Version) | Max = 5.53ε (Mean = 2.4ε) (<cmath>: Max = 2.05e+05ε (Mean = 6.87e+04ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 6.46ε (Mean = 2.38ε)) (Rmath 3.2.3: Max = 167ε (Mean = 56.5ε)) | Max = 5.53ε (Mean = 2.4ε) | Max = 4.61ε (Mean = 2.29ε) (<math.h>: Max = 5.37e+03ε (Mean = 1.81e+03ε))
Y1: Mathworld Data (Integer Version) | Max = 6.33ε (Mean = 2.25ε) (<cmath>: Max = 9.71e+03ε (Mean = 4.08e+03ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.51ε (Mean = 0.839ε)) (Rmath 3.2.3: Max = 193ε (Mean = 64.4ε)) | Max = 6.33ε (Mean = 2.29ε) | Max = 4.75ε (Mean = 1.72ε) (<math.h>: Max = 1.86e+04ε (Mean = 6.2e+03ε))
Yn: Mathworld Data (Integer Version) | Max = 55.2ε (Mean = 17.8ε) (<cmath>: Max = 2.2e+20ε (Mean = 6.97e+19ε) And other failures.) | Max = 0.993ε (Mean = 0.314ε) (GSL 2.1: Max = 2.41e+05ε (Mean = 7.62e+04ε)) (Rmath 3.2.3: Max = 1.24e+04ε (Mean = 4e+03ε)) | Max = 55.2ε (Mean = 17.8ε) | Max = 35ε (Mean = 11.9ε) (<math.h>: Max = 2.49e+05ε (Mean = 8.14e+04ε))

Table 23.24. Error rates for cyl_neumann_prime

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Y'0: Mathworld Data | Max = 0ε (Mean = 0ε) | Max = 6.33ε (Mean = 3.12ε) | Max = 6.33ε (Mean = 3.14ε) | Max = 4.75ε (Mean = 1.75ε)
Y'1: Mathworld Data | Max = 0.58ε (Mean = 0.193ε) | Max = 37.1ε (Mean = 12.8ε) | Max = 34ε (Mean = 11.8ε) | Max = 3.08ε (Mean = 1.2ε)
Y'n: Mathworld Data | Max = 2.05ε (Mean = 0.677ε) | Max = 56ε (Mean = 18.2ε) | Max = 56ε (Mean = 21.3ε) | Max = 563ε (Mean = 178ε)
Y'v: Mathworld Data | Max = 21.5ε (Mean = 6.49ε) | Max = 42.5ε (Mean = 13.4ε) | Max = 42.5ε (Mean = 13.6ε) | Max = 23.7ε (Mean = 10.1ε)
Y'v: Mathworld Data (large values) | Max = 0ε (Mean = 0ε) | Max = 1.57ε (Mean = 1.24ε) | Max = 1.57ε (Mean = 1.24ε) | Max = 0.627ε (Mean = 0.237ε)
Y'0 and Y'1: Random Data | Max = 0ε (Mean = 0ε) | Max = 23.8ε (Mean = 3.69ε) | Max = 23.8ε (Mean = 3.69ε) | Max = 5.95ε (Mean = 1.36ε)
Y'n: Random Data | Max = 1.53ε (Mean = 0.0885ε) | Max = 2.35e+03ε (Mean = 136ε) | Max = 2.35e+03ε (Mean = 136ε) | Max = 621ε (Mean = 36ε)
Y'v: Random Data | Max = 56.8ε (Mean = 2.59ε) | Max = 1.16e+05ε (Mean = 5.28e+03ε) | Max = 1.16e+05ε (Mean = 5.28e+03ε) | Max = 3.23e+04ε (Mean = 1.13e+03ε)

Table 23.25. Error rates for cyl_neumann_prime (integer orders)

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Y'0: Mathworld Data (Integer Version) | Max = 0ε (Mean = 0ε) | Max = 6.33ε (Mean = 3.12ε) | Max = 6.33ε (Mean = 3.14ε) | Max = 4.75ε (Mean = 1.75ε)
Y'1: Mathworld Data (Integer Version) | Max = 0.58ε (Mean = 0.193ε) | Max = 37.1ε (Mean = 12.8ε) | Max = 34ε (Mean = 11.8ε) | Max = 3.08ε (Mean = 1.2ε)
Y'n: Mathworld Data (Integer Version) | Max = 2.05ε (Mean = 0.677ε) | Max = 56ε (Mean = 18.2ε) | Max = 56ε (Mean = 21.3ε) | Max = 563ε (Mean = 178ε)

Table 23.26. Error rates for digamma

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Digamma Function: Large Values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.84ε (Mean = 0.71ε)) (Rmath 3.2.3: Max = 1.18ε (Mean = 0.331ε)) | Max = 1.39ε (Mean = 0.413ε) | Max = 1.39ε (Mean = 0.413ε) | Max = 0.98ε (Mean = 0.369ε)
Digamma Function: Near the Positive Root | Max = 0.891ε (Mean = 0.0995ε) (GSL 2.1: Max = 135ε (Mean = 11.9ε)) (Rmath 3.2.3: Max = 2.02e+03ε (Mean = 256ε)) | Max = 1.37ε (Mean = 0.477ε) | Max = 1.31ε (Mean = 0.471ε) | Max = 0.997ε (Mean = 0.527ε)
Digamma Function: Near Zero | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0.953ε (Mean = 0.348ε)) (Rmath 3.2.3: Max = 1.17ε (Mean = 0.564ε)) | Max = 0.984ε (Mean = 0.361ε) | Max = 0.984ε (Mean = 0.361ε) | Max = 0.953ε (Mean = 0.337ε)
Digamma Function: Negative Values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 4.56e+04ε (Mean = 3.91e+03ε)) (Rmath 3.2.3: Max = 4.6e+04ε (Mean = 3.94e+03ε)) | Max = 180ε (Mean = 13ε) | Max = 180ε (Mean = 13ε) | Max = 214ε (Mean = 16.1ε)
Digamma Function: Values near 0 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0.866ε (Mean = 0.387ε)) (Rmath 3.2.3: Max = 3.58e+05ε (Mean = 1.6e+05ε)) | Max = 1ε (Mean = 0.592ε) | Max = 1ε (Mean = 0.592ε) | Max = 0ε (Mean = 0ε)
Digamma Function: Integer arguments | Max = 0.992ε (Mean = 0.215ε) (GSL 2.1: Max = 1.18ε (Mean = 0.607ε)) (Rmath 3.2.3: Max = 4.33ε (Mean = 0.982ε)) | Max = 0.888ε (Mean = 0.403ε) | Max = 0.888ε (Mean = 0.403ε) | Max = 0.992ε (Mean = 0.452ε)
Digamma Function: Half integer arguments | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.09ε (Mean = 0.531ε)) (Rmath 3.2.3: Max = 46.2ε (Mean = 7.24ε)) | Max = 0.906ε (Mean = 0.409ε) | Max = 0.906ε (Mean = 0.409ε) | Max = 0.78ε (Mean = 0.314ε)

Table 23.27. Error rates for ellint_1

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Elliptic Integral F: Mathworld Data | Max = 0.94ε (Mean = 0.509ε) (<cmath>: Max = +INFε (Mean = +INFε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0.919ε (Mean = 0.544ε)) | Max = 0.94ε (Mean = 0.509ε) | Max = 0.919ε (Mean = 0.542ε)
Elliptic Integral F: Random Data | Max = 1.57ε (Mean = 0.56ε) (<cmath>: Max = 2.56ε (Mean = 0.816ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 2.99ε (Mean = 0.797ε)) | Max = 1.57ε (Mean = 0.561ε) | Max = 2.26ε (Mean = 0.631ε)

Table 23.28. Error rates for ellint_1 (complete)

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Elliptic Integral K: Mathworld Data | Max = 0.887ε (Mean = 0.296ε) (<cmath>: Max = 1.19ε (Mean = 0.765ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0.623ε (Mean = 0.393ε)) | Max = 0.887ε (Mean = 0.296ε) | Max = 0.915ε (Mean = 0.547ε)
Elliptic Integral K: Random Data | Max = 1.27ε (Mean = 0.473ε) (<cmath>: Max = 2.19ε (Mean = 0.694ε)) | Max = 0.851ε (Mean = 0.0851ε) (GSL 2.1: Max = 2.32ε (Mean = 0.688ε)) | Max = 1.27ε (Mean = 0.473ε) | Max = 0.958ε (Mean = 0.408ε)

Table 23.29. Error rates for ellint_2

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Elliptic Integral E: Mathworld Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0.63ε (Mean = 0.325ε)) | Max = 0.656ε (Mean = 0.317ε) (<cmath>: Max = +INFε (Mean = +INFε) And other failures.) | Max = 0.656ε (Mean = 0.317ε) | Max = 1.31ε (Mean = 0.727ε)
Elliptic Integral E: Random Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 4.4ε (Mean = 1.16ε)) | Max = 2.05ε (Mean = 0.632ε) (<cmath>: Max = 3.08e+04ε (Mean = 3.84e+03ε)) | Max = 2.05ε (Mean = 0.632ε) | Max = 2.23ε (Mean = 0.639ε)
Elliptic Integral E: Small Angles | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0.5ε (Mean = 0.118ε)) | Max = 1ε (Mean = 0.283ε) (<cmath>: Max = 2ε (Mean = 0.333ε)) | Max = 1ε (Mean = 0.283ε) | Max = 1ε (Mean = 0.421ε)

Table 23.30. Error rates for ellint_2 (complete)

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Elliptic Integral E: Mathworld Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 3.09ε (Mean = 1.04ε) And other failures.) | Max = 0.836ε (Mean = 0.469ε) (<cmath>: Max = 170ε (Mean = 55.1ε)) | Max = 0.836ε (Mean = 0.469ε) | Max = 1.3ε (Mean = 0.615ε)
Elliptic Integral E: Random Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 4.34ε (Mean = 1.18ε)) | Max = 1.97ε (Mean = 0.629ε) (<cmath>: Max = 2.49e+04ε (Mean = 3.39e+03ε)) | Max = 1.97ε (Mean = 0.629ε) | Max = 1.71ε (Mean = 0.553ε)

Table 23.31. Error rates for ellint_3

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Elliptic Integral PI: Mathworld Data | Max = 475ε (Mean = 86.3ε) (<cmath>: Max = +INFε (Mean = +INFε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.48e+05ε (Mean = 2.54e+04ε) And other failures.) | Max = 475ε (Mean = 86.3ε) | Max = 565ε (Mean = 102ε)
Elliptic Integral PI: Random Data | Max = 4.54ε (Mean = 0.895ε) (<cmath>: Max = 3.37e+20ε (Mean = 3.47e+19ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 633ε (Mean = 50.1ε)) | Max = 4.49ε (Mean = 0.885ε) | Max = 8.33ε (Mean = 0.971ε)
Elliptic Integral PI: Large Random Data | Max = 3.7ε (Mean = 0.893ε) (<cmath>: Max = 2.52e+18ε (Mean = 4.83e+17ε) And other failures.) | Max = 0.557ε (Mean = 0.0389ε) (GSL 2.1: Max = 40.1ε (Mean = 7.77ε)) | Max = 3.7ε (Mean = 0.892ε) | Max = 2.86ε (Mean = 0.944ε)

Table 23.32. Error rates for ellint_3 (complete)

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Complete Elliptic Integral PI: Mathworld Data | Max = 1.4ε (Mean = 0.575ε) (<cmath>: Max = 6.31e+20ε (Mean = 1.53e+20ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 6.33e+04ε (Mean = 1.54e+04ε) And other failures.) | Max = 1.4ε (Mean = 0.575ε) | Max = 0.971ε (Mean = 0.464ε)
Complete Elliptic Integral PI: Random Data | Max = 2.45ε (Mean = 0.696ε) (<cmath>: Max = 8.78e+20ε (Mean = 1.02e+20ε) And other failures.) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 24ε (Mean = 2.99ε)) | Max = 2.4ε (Mean = 0.677ε) | Max = 2.46ε (Mean = 0.657ε)

Table 23.33. Error rates for ellint_d

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Elliptic Integral E: Mathworld Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0.862ε (Mean = 0.568ε)) | Max = 1.3ε (Mean = 0.813ε) | Max = 1.3ε (Mean = 0.813ε) | Max = 0.862ε (Mean = 0.457ε)
Elliptic Integral D: Random Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 3.01ε (Mean = 0.928ε)) | Max = 2.51ε (Mean = 0.883ε) | Max = 2.51ε (Mean = 0.883ε) | Max = 2.87ε (Mean = 0.805ε)

Table 23.34. Error rates for ellint_d (complete)

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Elliptic Integral E: Mathworld Data | Max = 0.637ε (Mean = 0.368ε) | Max = 1.27ε (Mean = 0.735ε) | Max = 1.27ε (Mean = 0.735ε) | Max = 0.637ε (Mean = 0.368ε)
Elliptic Integral D: Random Data | Max = 0ε (Mean = 0ε) | Max = 1.27ε (Mean = 0.334ε) | Max = 1.27ε (Mean = 0.334ε) | Max = 1.27ε (Mean = 0.355ε)

Table 23.35. Error rates for ellint_rc

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | ---
RC: Random data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 2.4ε (Mean = 0.624ε)) | Max = 0.995ε (Mean = 0.433ε) | Max = 0.962ε (Mean = 0.407ε)

Table 23.36. Error rates for ellint_rd

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | ---
RD: Random data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 2.59ε (Mean = 0.878ε)) | Max = 2.73ε (Mean = 0.831ε) | Max = 2.16ε (Mean = 0.803ε)
RD: y = z | Max = 0.896ε (Mean = 0.022ε) (GSL 2.1: Max = 2.88ε (Mean = 0.839ε)) | Max = 2.65ε (Mean = 0.82ε) | Max = 16.5ε (Mean = 0.843ε)
RD: x = y | Max = 0.824ε (Mean = 0.0272ε) (GSL 2.1: Max = 3.74ε (Mean = 0.84ε)) | Max = 2.85ε (Mean = 0.865ε) | Max = 3.51ε (Mean = 0.816ε)
RD: x = 0, y = z | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 2ε (Mean = 0.656ε)) | Max = 1.19ε (Mean = 0.522ε) | Max = 1.16ε (Mean = 0.497ε)
RD: x = y = z | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.03ε (Mean = 0.418ε)) | Max = 0.998ε (Mean = 0.387ε) | Max = 1.03ε (Mean = 0.418ε)
RD: x = 0 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 2.85ε (Mean = 0.781ε)) | Max = 2.79ε (Mean = 0.883ε) | Max = 2.64ε (Mean = 0.894ε)

Table 23.37. Error rates for ellint_rf

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | ---
RF: Random data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 2.73ε (Mean = 0.804ε)) | Max = 2.54ε (Mean = 0.674ε) | Max = 2.02ε (Mean = 0.677ε)
RF: x = y = z | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0.999ε (Mean = 0.34ε)) | Max = 0.991ε (Mean = 0.345ε) | Max = 0.999ε (Mean = 0.34ε)
RF: x = y or y = z or x = z | Max = 0.536ε (Mean = 0.00658ε) (GSL 2.1: Max = 2.89ε (Mean = 0.749ε)) | Max = 1.95ε (Mean = 0.418ε) | Max = 1.21ε (Mean = 0.394ε)
RF: x = 0, y = z | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.29ε (Mean = 0.527ε)) | Max = 0.894ε (Mean = 0.338ε) | Max = 0.999ε (Mean = 0.407ε)
RF: z = 0 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 2.54ε (Mean = 0.781ε)) | Max = 1.7ε (Mean = 0.539ε) | Max = 1.89ε (Mean = 0.587ε)

Table 23.38. Error rates for ellint_rg

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | ---
RG: Random Data | Max = 0.983ε (Mean = 0.0172ε) (GSL 2.1: Max = 0.983ε (Mean = 0.0172ε)) | Max = 3.95ε (Mean = 0.951ε) | Max = 3.65ε (Mean = 0.929ε)
RG: two values 0 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0ε (Mean = 0ε)) | Max = 0ε (Mean = 0ε) | Max = 0ε (Mean = 0ε)
RG: All values the same or zero | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0ε (Mean = 0ε)) | Max = 0.992ε (Mean = 0.288ε) | Max = 1.06ε (Mean = 0.348ε)
RG: two values the same | Max = 0.594ε (Mean = 0.0103ε) (GSL 2.1: Max = 0.594ε (Mean = 0.0103ε)) | Max = 1.51ε (Mean = 0.404ε) | Max = 1.96ε (Mean = 0.374ε)
RG: one value zero | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0ε (Mean = 0ε)) | Max = 2.14ε (Mean = 0.722ε) | Max = 1.96ε (Mean = 0.674ε)

Table 23.39. Error rates for ellint_rj

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | ---
RJ: Random data | Max = 0.52ε (Mean = 0.0184ε) (GSL 2.1: Max = 3.57ε (Mean = 0.704ε) And other failures.) | Max = 186ε (Mean = 6.67ε) | Max = 215ε (Mean = 7.66ε)
RJ: 4 Equal Values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.03ε (Mean = 0.418ε)) | Max = 0.998ε (Mean = 0.387ε) | Max = 1.03ε (Mean = 0.418ε)
RJ: 3 Equal Values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 3.96ε (Mean = 1.06ε)) | Max = 20.8ε (Mean = 0.986ε) | Max = 39.9ε (Mean = 1.17ε)
RJ: 2 Equal Values | Max = 0.6ε (Mean = 0.0228ε) (GSL 2.1: Max = 2.57ε (Mean = 0.754ε)) | Max = 220ε (Mean = 6.64ε) | Max = 214ε (Mean = 5.28ε)
RJ: Equal z and p | Max = 0.742ε (Mean = 0.0166ε) (GSL 2.1: Max = 2.62ε (Mean = 0.699ε)) | Max = 17.2ε (Mean = 1.16ε) | Max = 16.1ε (Mean = 1.14ε)

Table 23.40. Error rates for erf

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Erf Function: Small Values | Max = 0.925ε (Mean = 0.193ε) (<cmath>: Max = 0.944ε (Mean = 0.191ε)) (<math.h>: Max = 0.944ε (Mean = 0.191ε)) | Max = 0.841ε (Mean = 0.0687ε) (GSL 2.1: Max = 2.06ε (Mean = 0.319ε)) | Max = 0.925ε (Mean = 0.193ε) (<math.h>: Max = 0.944ε (Mean = 0.194ε)) | Max = 0.996ε (Mean = 0.182ε) (<math.h>: Max = 1.57ε (Mean = 0.317ε))
Erf Function: Medium Values | Max = 1.5ε (Mean = 0.193ε) (<cmath>: Max = 0.921ε (Mean = 0.0723ε)) (<math.h>: Max = 0.921ε (Mean = 0.0723ε)) | Max = 1ε (Mean = 0.119ε) (GSL 2.1: Max = 2.31ε (Mean = 0.368ε)) | Max = 1.5ε (Mean = 0.197ε) (<math.h>: Max = 0.921ε (Mean = 0.071ε)) | Max = 1ε (Mean = 0.171ε) (<math.h>: Max = 1.19ε (Mean = 0.244ε))
Erf Function: Large Values | Max = 0ε (Mean = 0ε) (<cmath>: Max = 0ε (Mean = 0ε)) (<math.h>: Max = 0ε (Mean = 0ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0ε (Mean = 0ε)) | Max = 0ε (Mean = 0ε) (<math.h>: Max = 0ε (Mean = 0ε)) | Max = 0ε (Mean = 0ε) (<math.h>: Max = 0ε (Mean = 0ε))

Table 23.41. Error rates for erf_inv

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Inverse Erf Function | Max = 0ε (Mean = 0ε) | Max = 0.996ε (Mean = 0.389ε) | Max = 1.08ε (Mean = 0.395ε) | Max = 1.09ε (Mean = 0.502ε)

Table 23.42. Error rates for erfc

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Erf Function: Small Values | Max = 0ε (Mean = 0ε) (<cmath>: Max = 0ε (Mean = 0ε)) (<math.h>: Max = 0ε (Mean = 0ε)) | Max = 0.658ε (Mean = 0.0537ε) (GSL 2.1: Max = 1.01ε (Mean = 0.485ε)) | Max = 0ε (Mean = 0ε) (<math.h>: Max = 0ε (Mean = 0ε)) | Max = 0ε (Mean = 0ε) (<math.h>: Max = 0ε (Mean = 0ε))
Erf Function: Medium Values | Max = 1.76ε (Mean = 0.365ε) (<cmath>: Max = 1.35ε (Mean = 0.307ε)) (<math.h>: Max = 1.35ε (Mean = 0.307ε)) | Max = 0.983ε (Mean = 0.213ε) (GSL 2.1: Max = 2.64ε (Mean = 0.662ε)) | Max = 1.76ε (Mean = 0.38ε) (<math.h>: Max = 2.81ε (Mean = 0.739ε)) | Max = 1.65ε (Mean = 0.373ε) (<math.h>: Max = 2.36ε (Mean = 0.539ε))
Erf Function: Large Values | Max = 1.57ε (Mean = 0.542ε) (<cmath>: Max = 1.26ε (Mean = 0.441ε)) (<math.h>: Max = 1.26ε (Mean = 0.441ε)) | Max = 0.868ε (Mean = 0.147ε) (GSL 2.1: Max = 3.9ε (Mean = 0.472ε)) | Max = 1.57ε (Mean = 0.564ε) (<math.h>: Max = 4.91ε (Mean = 1.54ε)) | Max = 1.14ε (Mean = 0.248ε) (<math.h>: Max = 1.84ε (Mean = 0.331ε))

Table 23.43. Error rates for erfc_inv

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Inverse Erfc Function | Max = 0ε (Mean = 0ε) | Max = 0.996ε (Mean = 0.397ε) | Max = 1.08ε (Mean = 0.403ε) | Max = 1ε (Mean = 0.491ε)
Inverse Erfc Function: extreme values |  | Max = 1.62ε (Mean = 0.383ε) | Max = 1.62ε (Mean = 0.383ε) |

Table 23.44. Error rates for expint (Ei)

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Exponential Integral Ei | Max = 5.05ε (Mean = 0.821ε) (<cmath>: Max = 14.1ε (Mean = 2.43ε) And other failures.) | Max = 0.994ε (Mean = 0.142ε) (GSL 2.1: Max = 8.96ε (Mean = 0.703ε)) | Max = 5.05ε (Mean = 0.835ε) | Max = 1.43ε (Mean = 0.54ε)
Exponential Integral Ei: double exponent range | Max = 1.72ε (Mean = 0.593ε) (<cmath>: Max = 3.11ε (Mean = 1.13ε)) | Max = 0.998ε (Mean = 0.156ε) (GSL 2.1: Max = 1.5ε (Mean = 0.612ε)) | Max = 1.72ε (Mean = 0.607ε) | Max = 1.7ε (Mean = 0.66ε)
Exponential Integral Ei: long exponent range | Max = 1.98ε (Mean = 0.595ε) (<cmath>: Max = 1.93ε (Mean = 0.855ε)) |  | Max = 1.98ε (Mean = 0.575ε) |

Table 23.45. Error rates for expint (En)

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Exponential Integral En | Max = 0.589ε (Mean = 0.0331ε) (GSL 2.1: Max = 58.5ε (Mean = 17.1ε)) | Max = 9.97ε (Mean = 2.13ε) | Max = 9.97ε (Mean = 2.13ε) | Max = 7.16ε (Mean = 1.85ε)
Exponential Integral En: small z values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 115ε (Mean = 23.6ε)) | Max = 1.99ε (Mean = 0.559ε) | Max = 1.99ε (Mean = 0.559ε) | Max = 2.62ε (Mean = 0.531ε)
Exponential Integral E1 | Max = 0.556ε (Mean = 0.0625ε) (GSL 2.1: Max = 0.988ε (Mean = 0.469ε)) | Max = 0.965ε (Mean = 0.414ε) | Max = 0.965ε (Mean = 0.408ε) | Max = 0.988ε (Mean = 0.486ε)

Table 23.46. Error rates for expm1

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Random test data | Max = 0.992ε (Mean = 0.402ε) (<cmath>: Max = 0.992ε (Mean = 0.402ε)) (<math.h>: Max = 0.992ε (Mean = 0.402ε)) | Max = 0.793ε (Mean = 0.126ε) (Rmath 3.2.3: Max = 0.793ε (Mean = 0.126ε)) | Max = 1.31ε (Mean = 0.428ε) (<math.h>: Max = 0.996ε (Mean = 0.426ε)) | Max = 1.31ε (Mean = 0.496ε) (<math.h>: Max = 1.31ε (Mean = 0.496ε))

Table 23.47. Error rates for gamma_p

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
tgamma(a, z) medium values | Max = 0.955ε (Mean = 0.05ε) (GSL 2.1: Max = 342ε (Mean = 45.8ε)) (Rmath 3.2.3: Max = 389ε (Mean = 44ε)) | Max = 41.6ε (Mean = 8.09ε) | Max = 239ε (Mean = 30.2ε) | Max = 35.1ε (Mean = 6.98ε)
tgamma(a, z) small values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 4.82ε (Mean = 0.758ε)) (Rmath 3.2.3: Max = 1.01ε (Mean = 0.306ε)) | Max = 2ε (Mean = 0.464ε) | Max = 2ε (Mean = 0.461ε) | Max = 1.54ε (Mean = 0.439ε)
tgamma(a, z) large values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.02e+03ε (Mean = 105ε)) (Rmath 3.2.3: Max = 1.11e+03ε (Mean = 67.5ε)) | Max = 3.08e+04ε (Mean = 1.86e+03ε) | Max = 3.02e+04ε (Mean = 1.91e+03ε) | Max = 243ε (Mean = 20.2ε)
tgamma(a, z) integer and half integer values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 128ε (Mean = 22.6ε)) (Rmath 3.2.3: Max = 66.2ε (Mean = 12.2ε)) | Max = 11.8ε (Mean = 2.66ε) | Max = 71.6ε (Mean = 9.47ε) | Max = 13ε (Mean = 2.97ε)

Table 23.48. Error rates for gamma_p_inv

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
incomplete gamma inverse(a, z) medium values | Max = 0.993ε (Mean = 0.15ε) (Rmath 3.2.3: Max = 4.88ε (Mean = 0.868ε)) | Max = 1.8ε (Mean = 0.406ε) | Max = 1.89ε (Mean = 0.466ε) | Max = 1.71ε (Mean = 0.34ε)
incomplete gamma inverse(a, z) large values | Max = 0ε (Mean = 0ε) (Rmath 3.2.3: Max = 0.816ε (Mean = 0.0874ε)) | Max = 0.509ε (Mean = 0.0447ε) | Max = 0.509ε (Mean = 0.0447ε) | Max = 0.924ε (Mean = 0.108ε)
incomplete gamma inverse(a, z) small values | Max = 441ε (Mean = 53.9ε) (Rmath 3.2.3: Max = 547ε (Mean = 61.6ε)) | Max = 9.17e+03ε (Mean = 1.45e+03ε) | Max = 1.09e+04ε (Mean = 1.3e+03ε) | Max = 1.1e+03ε (Mean = 131ε)

Table 23.49. Error rates for gamma_p_inva

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Incomplete gamma inverses. | Max = 0ε (Mean = 0ε) | Max = 7.87ε (Mean = 1.15ε) | Max = 4.08ε (Mean = 1.12ε) | Max = 4.92ε (Mean = 1.03ε)

Table 23.50. Error rates for gamma_q

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
tgamma(a, z) medium values | Max = 0.927ε (Mean = 0.035ε) (GSL 2.1: Max = 201ε (Mean = 13.5ε)) (Rmath 3.2.3: Max = 131ε (Mean = 12.7ε)) | Max = 32.3ε (Mean = 6.61ε) | Max = 199ε (Mean = 26.6ε) | Max = 23.7ε (Mean = 4ε)
tgamma(a, z) small values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.38e+10ε (Mean = 1.05e+09ε)) (Rmath 3.2.3: Max = 65.6ε (Mean = 11ε)) | Max = 2.45ε (Mean = 0.885ε) | Max = 2.45ε (Mean = 0.819ε) | Max = 2.26ε (Mean = 0.74ε)
tgamma(a, z) large values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 2.71e+04ε (Mean = 2.16e+03ε)) (Rmath 3.2.3: Max = 1.02e+03ε (Mean = 62.7ε)) | Max = 6.82e+03ε (Mean = 414ε) | Max = 1.15e+04ε (Mean = 733ε) | Max = 469ε (Mean = 31.5ε)
tgamma(a, z) integer and half integer values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 118ε (Mean = 12.5ε)) (Rmath 3.2.3: Max = 138ε (Mean = 16.9ε)) | Max = 11.1ε (Mean = 2.07ε) | Max = 54.7ε (Mean = 6.16ε) | Max = 8.72ε (Mean = 1.48ε)

Table 23.51. Error rates for gamma_q_inv

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
incomplete gamma inverse(a, z) medium values | Max = 0.912ε (Mean = 0.154ε) (Rmath 3.2.3: Max = 4.66ε (Mean = 0.792ε)) | Max = 6.2ε (Mean = 0.627ε) | Max = 6.2ε (Mean = 0.683ε) | Max = 2.88ε (Mean = 0.469ε)
incomplete gamma inverse(a, z) large values | Max = 0.894ε (Mean = 0.0915ε) (Rmath 3.2.3: Max = 0.894ε (Mean = 0.106ε)) | Max = 0ε (Mean = 0ε) | Max = 0ε (Mean = 0ε) | Max = 0.814ε (Mean = 0.0856ε)
incomplete gamma inverse(a, z) small values | Max = 292ε (Mean = 36.4ε) (Rmath 3.2.3: Max = 415ε (Mean = 48.7ε)) | Max = 8.28e+03ε (Mean = 1.09e+03ε) | Max = 8.98e+03ε (Mean = 877ε) | Max = 451ε (Mean = 64.7ε)

Table 23.52. Error rates for gamma_q_inva

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Incomplete gamma inverses. | Max = 0ε (Mean = 0ε) | Max = 8.42ε (Mean = 1.3ε) | Max = 7.86ε (Mean = 1.24ε) | Max = 5.05ε (Mean = 1.08ε)

Table 23.53. Error rates for hermite

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Hermite Polynomials | Max = 0ε (Mean = 0ε) | Max = 6.24ε (Mean = 2.07ε) | Max = 6.24ε (Mean = 2.07ε) | Max = 4.46ε (Mean = 1.41ε)

Table 23.54. Error rates for heuman_lambda

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Elliptic Integral Jacobi Zeta: Mathworld Data | Max = 0ε (Mean = 0ε) | Max = 1.89ε (Mean = 0.887ε) | Max = 1.89ε (Mean = 0.887ε) | Max = 1.08ε (Mean = 0.734ε)
Elliptic Integral Heuman Lambda: Random Data | Max = 0ε (Mean = 0ε) | Max = 3.82ε (Mean = 0.609ε) | Max = 3.82ε (Mean = 0.608ε) | Max = 2.12ε (Mean = 0.588ε)

Table 23.55. Error rates for ibeta

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Incomplete Beta Function: Small Values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 682ε (Mean = 32.6ε)) (Rmath 3.2.3: Max = 22.9ε (Mean = 3.35ε)) | Max = 8.97ε (Mean = 2.09ε) | Max = 21.3ε (Mean = 2.75ε) | Max = 8.4ε (Mean = 1.93ε)
Incomplete Beta Function: Medium Values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 690ε (Mean = 151ε)) (Rmath 3.2.3: Max = 232ε (Mean = 27.9ε)) | Max = 50ε (Mean = 12.1ε) | Max = 124ε (Mean = 18.4ε) | Max = 106ε (Mean = 16.3ε)
Incomplete Beta Function: Large and Diverse Values | Max = 1.26ε (Mean = 0.063ε) (GSL 2.1: Max = 3.9e+05ε (Mean = 1.82e+04ε) And other failures.) (Rmath 3.2.3: Max = 574ε (Mean = 49.4ε)) | Max = 1.96e+04ε (Mean = 997ε) | Max = 4.98e+04ε (Mean = 2.07e+03ε) | Max = 1.32e+03ε (Mean = 68.5ε)
Incomplete Beta Function: Small Integer Values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 254ε (Mean = 50.9ε)) (Rmath 3.2.3: Max = 62.2ε (Mean = 8.95ε)) | Max = 4.45ε (Mean = 0.814ε) | Max = 44.5ε (Mean = 10.1ε) | Max = 3.85ε (Mean = 0.791ε)

Table 23.56. Error rates for ibeta_inv

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Inverse incomplete beta | Max = 11ε (Mean = 0.345ε) (Rmath 3.2.3: Max = 1.14e+121ε (Mean = 3.28e+119ε) And other failures.) | Max = 3.8e+04ε (Mean = 2.66e+03ε) | Max = 4.07e+04ε (Mean = 2.86e+03ε) | Max = 8.59e+03ε (Mean = 277ε)

Table 23.57. Error rates for ibeta_inva

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Inverse incomplete beta | Max = 0.602ε (Mean = 0.0239ε) | Max = 377ε (Mean = 24.4ε) | Max = 438ε (Mean = 31.3ε) | Max = 242ε (Mean = 22.9ε)

Table 23.58. Error rates for ibeta_invb

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Inverse incomplete beta | Max = 0.765ε (Mean = 0.0422ε) | Max = 407ε (Mean = 27.2ε) | Max = 407ε (Mean = 24.4ε) | Max = 409ε (Mean = 19.3ε)

Table 23.59. Error rates for ibetac

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Incomplete Beta Function: Small Values | Max = 0ε (Mean = 0ε) (Rmath 3.2.3: Max = 22.4ε (Mean = 3.67ε)) | Max = 10.6ε (Mean = 2.22ε) | Max = 13.8ε (Mean = 2.68ε) | Max = 6.94ε (Mean = 1.71ε)
Incomplete Beta Function: Medium Values | Max = 0ε (Mean = 0ε) (Rmath 3.2.3: Max = 204ε (Mean = 25.8ε)) | Max = 73.9ε (Mean = 11.2ε) | Max = 132ε (Mean = 19.8ε) | Max = 56.7ε (Mean = 14.3ε)
Incomplete Beta Function: Large and Diverse Values | Max = 0.981ε (Mean = 0.0573ε) (Rmath 3.2.3: Max = 889ε (Mean = 68.4ε)) | Max = 3.45e+04ε (Mean = 1.32e+03ε) | Max = 6.31e+04ε (Mean = 2.04e+03ε) | Max = 1.88e+03ε (Mean = 82.7ε)
Incomplete Beta Function: Small Integer Values | Max = 0ε (Mean = 0ε) (Rmath 3.2.3: Max = 84.6ε (Mean = 18ε)) | Max = 5.34ε (Mean = 1.11ε) | Max = 107ε (Mean = 17.1ε) | Max = 6.37ε (Mean = 1.03ε)

Table 23.60. Error rates for ibetac_inv

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Inverse incomplete beta | Max = 0.977ε (Mean = 0.0976ε) (Rmath 3.2.3: Max = 3.01e+132ε (Mean = 8.65e+130ε) And other failures.) | Max = 4.88e+04ε (Mean = 3.16e+03ε) | Max = 5.05e+04ε (Mean = 3.33e+03ε) | Max = 2.93e+03ε (Mean = 198ε)

Table 23.61. Error rates for ibetac_inva

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Inverse incomplete beta | Max = 0.683ε (Mean = 0.0314ε) | Max = 382ε (Mean = 22.2ε) | Max = 315ε (Mean = 23.7ε) | Max = 408ε (Mean = 26.7ε)

Table 23.62. Error rates for ibetac_invb

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Inverse incomplete beta | Max = 0.724ε (Mean = 0.0303ε) | Max = 317ε (Mean = 19.8ε) | Max = 369ε (Mean = 22.6ε) | Max = 271ε (Mean = 16.4ε)

Table 23.63. Error rates for jacobi_cn

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Jacobi Elliptic: Mathworld Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 17.3ε (Mean = 4.29ε) And other failures.) | Max = 71.6ε (Mean = 19.3ε) | Max = 71.6ε (Mean = 19.4ε) | Max = 45.8ε (Mean = 11.4ε)
Jacobi Elliptic: Random Data | Max = 0.816ε (Mean = 0.0563ε) (GSL 2.1: Max = 2.43ε (Mean = 0.803ε)) | Max = 1.68ε (Mean = 0.443ε) | Max = 1.68ε (Mean = 0.454ε) | Max = 1.83ε (Mean = 0.455ε)
Jacobi Elliptic: Random Small Values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 55.2ε (Mean = 1.64ε) And other failures.) | Max = 10.4ε (Mean = 0.594ε) | Max = 10.4ε (Mean = 0.602ε) | Max = 26.2ε (Mean = 1.17ε)
Jacobi Elliptic: Modulus near 1 | Max = 0.919ε (Mean = 0.127ε) (GSL 2.1: Max = 0ε (Mean = 0ε) And other failures.) | Max = 675ε (Mean = 87.1ε) | Max = 675ε (Mean = 86.8ε) | Max = 513ε (Mean = 126ε)
Jacobi Elliptic: Large Phi | Max = 14.2ε (Mean = 0.927ε) (GSL 2.1: Max = 5.92e+03ε (Mean = 477ε)) | Max = 2.97e+04ε (Mean = 1.9e+03ε) | Max = 2.97e+04ε (Mean = 1.9e+03ε) | Max = 3.27e+04ε (Mean = 1.93e+03ε)

Table 23.64. Error rates for jacobi_dn

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Jacobi Elliptic: Mathworld Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 2.82ε (Mean = 1.18ε) And other failures.) | Max = 49ε (Mean = 14ε) | Max = 49ε (Mean = 14ε) | Max = 34.3ε (Mean = 8.71ε)
Jacobi Elliptic: Random Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 3ε (Mean = 0.61ε)) | Max = 1.53ε (Mean = 0.473ε) | Max = 1.53ε (Mean = 0.481ε) | Max = 1.52ε (Mean = 0.466ε)
Jacobi Elliptic: Random Small Values | Max = 0.5ε (Mean = 0.0122ε) (GSL 2.1: Max = 1.5ε (Mean = 0.391ε) And other failures.) | Max = 22.4ε (Mean = 0.777ε) | Max = 22.4ε (Mean = 0.763ε) | Max = 16.1ε (Mean = 0.685ε)
Jacobi Elliptic: Modulus near 1 | Max = 2.28ε (Mean = 0.194ε) (GSL 2.1: Max = 0ε (Mean = 0ε) And other failures.) | Max = 3.75e+03ε (Mean = 293ε) | Max = 3.75e+03ε (Mean = 293ε) | Max = 6.24e+03ε (Mean = 482ε)
Jacobi Elliptic: Large Phi | Max = 14.1ε (Mean = 0.897ε) (GSL 2.1: Max = 121ε (Mean = 22ε)) | Max = 2.82e+04ε (Mean = 1.79e+03ε) | Max = 2.82e+04ε (Mean = 1.79e+03ε) | Max = 1.67e+04ε (Mean = 1e+03ε)

Table 23.65. Error rates for jacobi_sn

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Jacobi Elliptic: Mathworld Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 588ε (Mean = 146ε) And other failures.) | Max = 341ε (Mean = 80.7ε) | Max = 341ε (Mean = 80.7ε) | Max = 481ε (Mean = 113ε)
Jacobi Elliptic: Random Data | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 4.02ε (Mean = 1.07ε)) | Max = 2.01ε (Mean = 0.584ε) | Max = 2.01ε (Mean = 0.593ε) | Max = 1.92ε (Mean = 0.567ε)
Jacobi Elliptic: Random Small Values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 11.7ε (Mean = 1.65ε) And other failures.) | Max = 1.99ε (Mean = 0.347ε) | Max = 1.99ε (Mean = 0.347ε) | Max = 2.11ε (Mean = 0.385ε)
Jacobi Elliptic: Modulus near 1 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0ε (Mean = 0ε) And other failures.) | Max = 109ε (Mean = 7.35ε) | Max = 109ε (Mean = 7.38ε) | Max = 23.2ε (Mean = 1.85ε)
Jacobi Elliptic: Large Phi | Max = 12ε (Mean = 0.771ε) (GSL 2.1: Max = 4.54e+04ε (Mean = 2.63e+03ε)) | Max = 2.45e+04ε (Mean = 1.51e+03ε) | Max = 2.45e+04ε (Mean = 1.51e+03ε) | Max = 4.36e+04ε (Mean = 2.54e+03ε)

Table 23.66. Error rates for jacobi_zeta

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Elliptic Integral Jacobi Zeta: Mathworld Data | Max = 0ε (Mean = 0ε) | Max = 1.66ε (Mean = 0.48ε) | Max = 1.66ε (Mean = 0.48ε) | Max = 1.52ε (Mean = 0.357ε)
Elliptic Integral Jacobi Zeta: Random Data | Max = 0ε (Mean = 0ε) | Max = 2.99ε (Mean = 0.824ε) | Max = 3.96ε (Mean = 1.06ε) | Max = 3.89ε (Mean = 0.824ε)
Elliptic Integral Jacobi Zeta: Large Phi Values | Max = 0ε (Mean = 0ε) | Max = 2.92ε (Mean = 0.951ε) | Max = 3.05ε (Mean = 1.13ε) | Max = 2.52ε (Mean = 0.977ε)

Table 23.67. Error rates for laguerre(n, m, x)

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Associated Laguerre Polynomials | Max = 0.84ε (Mean = 0.0358ε) (GSL 2.1: Max = 434ε (Mean = 10.7ε)) | Max = 167ε (Mean = 6.38ε) (<cmath>: Max = 206ε (Mean = 6.86ε)) | Max = 167ε (Mean = 6.38ε) | Max = 434ε (Mean = 11.1ε)

Table 23.68. Error rates for laguerre(n, x)

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Laguerre Polynomials | Max = 6.82ε (Mean = 0.408ε) (GSL 2.1: Max = 3.1e+03ε (Mean = 185ε)) | Max = 1.39e+04ε (Mean = 828ε) (<cmath>: Max = 4.2e+03ε (Mean = 251ε)) | Max = 1.39e+04ε (Mean = 828ε) | Max = 3.1e+03ε (Mean = 185ε)

Table 23.69. Error rates for legendre_p

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Legendre Polynomials: Small Values | Max = 0.732ε (Mean = 0.0619ε) (GSL 2.1: Max = 211ε (Mean = 20.4ε)) | Max = 69.2ε (Mean = 9.58ε) (<cmath>: Max = 124ε (Mean = 13.2ε)) | Max = 69.2ε (Mean = 9.58ε) | Max = 211ε (Mean = 20.4ε)
Legendre Polynomials: Large Values | Max = 0.632ε (Mean = 0.0693ε) (GSL 2.1: Max = 300ε (Mean = 33.2ε)) | Max = 699ε (Mean = 59.6ε) (<cmath>: Max = 343ε (Mean = 32.1ε)) | Max = 699ε (Mean = 59.6ε) | Max = 300ε (Mean = 33.2ε)

Table 23.70. Error rates for legendre_p (associated)

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Associated Legendre Polynomials: Small Values | Max = 0.999ε (Mean = 0.05ε) (GSL 2.1: Max = 121ε (Mean = 6.75ε) And other failures.) | Max = 175ε (Mean = 9.88ε) (<cmath>: Max = 175ε (Mean = 9.36ε) And other failures.) | Max = 77.7ε (Mean = 5.59ε) | Max = 121ε (Mean = 7.14ε)

Table 23.71. Error rates for legendre_q

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Legendre Polynomials: Small Values | Max = 0.612ε (Mean = 0.0517ε) (GSL 2.1: Max = 46.4ε (Mean = 7.46ε)) | Max = 50.9ε (Mean = 9ε) | Max = 50.9ε (Mean = 8.98ε) | Max = 46.4ε (Mean = 7.32ε)
Legendre Polynomials: Large Values | Max = 2.49ε (Mean = 0.202ε) (GSL 2.1: Max = 4.6e+03ε (Mean = 366ε)) | Max = 5.98e+03ε (Mean = 478ε) | Max = 5.98e+03ε (Mean = 478ε) | Max = 4.6e+03ε (Mean = 366ε)

Table 23.72. Error rates for lgamma

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
factorials | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 33.6ε (Mean = 2.78ε)) (Rmath 3.2.3: Max = 1.55ε (Mean = 0.592ε)) | Max = 0.991ε (Mean = 0.308ε) (<cmath>: Max = 1.67ε (Mean = 0.487ε)) (<math.h>: Max = 1.67ε (Mean = 0.487ε)) | Max = 0.991ε (Mean = 0.383ε) (<math.h>: Max = 1.36ε (Mean = 0.476ε)) | Max = 0.914ε (Mean = 0.175ε) (<math.h>: Max = 0.958ε (Mean = 0.38ε))
near 0 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 5.21ε (Mean = 1.57ε)) (Rmath 3.2.3: Max = 0ε (Mean = 0ε)) | Max = 1.42ε (Mean = 0.566ε) (<cmath>: Max = 0.964ε (Mean = 0.543ε)) (<math.h>: Max = 0.964ε (Mean = 0.543ε)) | Max = 1.42ε (Mean = 0.566ε) (<math.h>: Max = 0.964ε (Mean = 0.543ε)) | Max = 0.964ε (Mean = 0.462ε) (<math.h>: Max = 0.962ε (Mean = 0.372ε))
near 1 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 442ε (Mean = 88.8ε)) (Rmath 3.2.3: Max = 7.99e+04ε (Mean = 1.68e+04ε)) | Max = 0.948ε (Mean = 0.36ε) (<cmath>: Max = 0.615ε (Mean = 0.096ε)) (<math.h>: Max = 0.615ε (Mean = 0.096ε)) | Max = 0.948ε (Mean = 0.36ε) (<math.h>: Max = 1.71ε (Mean = 0.581ε)) | Max = 0.867ε (Mean = 0.468ε) (<math.h>: Max = 0.906ε (Mean = 0.565ε))
near 2 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.17e+03ε (Mean = 274ε)) (Rmath 3.2.3: Max = 2.63e+05ε (Mean = 5.84e+04ε)) | Max = 0.878ε (Mean = 0.242ε) (<cmath>: Max = 0.741ε (Mean = 0.263ε)) (<math.h>: Max = 0.741ε (Mean = 0.263ε)) | Max = 0.878ε (Mean = 0.242ε) (<math.h>: Max = 0.598ε (Mean = 0.235ε)) | Max = 0.591ε (Mean = 0.159ε) (<math.h>: Max = 0.741ε (Mean = 0.473ε))
near -10 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 24.9ε (Mean = 4.6ε)) (Rmath 3.2.3: Max = 4.22ε (Mean = 1.26ε)) | Max = 3.81ε (Mean = 1.01ε) (<cmath>: Max = 0.997ε (Mean = 0.412ε)) (<math.h>: Max = 0.997ε (Mean = 0.412ε)) | Max = 3.81ε (Mean = 1.01ε) (<math.h>: Max = 3.04ε (Mean = 1.01ε)) | Max = 4.22ε (Mean = 1.33ε) (<math.h>: Max = 0.997ε (Mean = 0.444ε))
near -55 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 7.02ε (Mean = 1.47ε)) (Rmath 3.2.3: Max = 250ε (Mean = 60.9ε)) | Max = 0.821ε (Mean = 0.513ε) (<cmath>: Max = 1.58ε (Mean = 0.672ε)) (<math.h>: Max = 1.58ε (Mean = 0.672ε)) | Max = 1.59ε (Mean = 0.587ε) (<math.h>: Max = 0.821ε (Mean = 0.674ε)) | Max = 0.821ε (Mean = 0.419ε) (<math.h>: Max = 249ε (Mean = 43.1ε))

Table 23.73. Error rates for log1p

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Random test data | Max = 0.818ε (Mean = 0.227ε) (<cmath>: Max = 0.818ε (Mean = 0.227ε)) (<math.h>: Max = 0.818ε (Mean = 0.227ε)) | Max = 0.846ε (Mean = 0.153ε) (Rmath 3.2.3: Max = 0.846ε (Mean = 0.153ε)) | Max = 2.3ε (Mean = 0.66ε) (<math.h>: Max = 0.818ε (Mean = 0.249ε)) | Max = 0.509ε (Mean = 0.057ε) (<math.h>: Max = 0.509ε (Mean = 0.057ε))

Table 23.74. Error rates for non central beta CDF

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Non Central Beta, medium parameters | Max = 0.998ε (Mean = 0.0649ε) (Rmath 3.2.3: Max = 1.46e+26ε (Mean = 3.5e+24ε) And other failures.) | Max = 824ε (Mean = 27.4ε) | Max = 832ε (Mean = 38.1ε) | Max = 242ε (Mean = 31ε)
Non Central Beta, large parameters | Max = 1.18ε (Mean = 0.175ε) (Rmath 3.2.3: Max = 1.01e+36ε (Mean = 1.19e+35ε) And other failures.) | Max = 2.5e+04ε (Mean = 3.78e+03ε) | Max = 2.57e+04ε (Mean = 4.45e+03ε) | Max = 3.66e+03ε (Mean = 500ε)

Table 23.75. Error rates for non central beta CDF complement

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Non Central Beta, medium parameters | Max = 0.998ε (Mean = 0.0936ε) (Rmath 3.2.3: Max = 7.5e+97ε (Mean = 1.37e+96ε) And other failures.) | Max = 396ε (Mean = 50.7ε) | Max = 554ε (Mean = 57.2ε) | Max = 624ε (Mean = 62.7ε)
Non Central Beta, large parameters | Max = 0.986ε (Mean = 0.188ε) (Rmath 3.2.3: Max = +INFε (Mean = +INFε) And other failures.) | Max = 6.83e+03ε (Mean = 993ε) | Max = 3.56e+03ε (Mean = 707ε) | Max = 1.25e+04ε (Mean = 1.49e+03ε)

Table 23.76. Error rates for non central chi squared CDF

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Non Central Chi Squared, medium parameters | Max = 0.99ε (Mean = 0.0544ε) (Rmath 3.2.3: Max = 727ε (Mean = 121ε)) | Max = 46.5ε (Mean = 10.3ε) | Max = 115ε (Mean = 13.9ε) | Max = 48.9ε (Mean = 10ε)
Non Central Chi Squared, large parameters | Max = 1.07ε (Mean = 0.102ε) (Rmath 3.2.3: Max = 3.27e+08ε (Mean = 2.23e+07ε)) | Max = 3.07e+03ε (Mean = 336ε) | Max = 6.17e+03ε (Mean = 677ε) | Max = 9.79e+03ε (Mean = 723ε)

Table 23.77. Error rates for non central chi squared CDF complement

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Non Central Chi Squared, medium parameters | Max = 0.96ε (Mean = 0.0635ε) (Rmath 3.2.3: Max = +INFε (Mean = +INFε) And other failures.) | Max = 107ε (Mean = 17.2ε) | Max = 171ε (Mean = 22.8ε) | Max = 98.6ε (Mean = 15.8ε)
Non Central Chi Squared, large parameters | Max = 2.11ε (Mean = 0.278ε) (Rmath 3.2.3: Max = +INFε (Mean = +INFε) And other failures.) | Max = 5.02e+03ε (Mean = 630ε) | Max = 5.1e+03ε (Mean = 577ε) | Max = 5.43e+03ε (Mean = 705ε)

Table 23.78. Error rates for non central t CDF

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Non Central T | Max = 0.796ε (Mean = 0.0691ε) (Rmath 3.2.3: Max = 5.28e+15ε (Mean = 8.49e+14ε) And other failures.) | Max = 139ε (Mean = 31ε) | Max = 145ε (Mean = 30.9ε) | Max = 135ε (Mean = 32.1ε)
Non Central T (small non-centrality) | Max = 0ε (Mean = 0ε) (Rmath 3.2.3: Max = 2.09e+03ε (Mean = 244ε)) | Max = 3.86ε (Mean = 1.4ε) | Max = 9.15ε (Mean = 2.13ε) | Max = 6.17ε (Mean = 1.45ε)
Non Central T (large parameters) | Max = 257ε (Mean = 72.1ε) (Rmath 3.2.3: Max = 2.46ε (Mean = 0.657ε)) | Max = 5.26e+05ε (Mean = 1.48e+05ε) | Max = 5.24e+05ε (Mean = 1.47e+05ε) | Max = 286ε (Mean = 62.8ε)

Table 23.79. Error rates for non central t CDF complement

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Non Central T | Max = 0.707ε (Mean = 0.0497ε) (Rmath 3.2.3: Max = 6.19e+15ε (Mean = 6.72e+14ε) And other failures.) | Max = 201ε (Mean = 31.7ε) | Max = 340ε (Mean = 43.2ε) | Max = 154ε (Mean = 32.1ε)
Non Central T (small non-centrality) | Max = 0ε (Mean = 0ε) (Rmath 3.2.3: Max = 1.87e+03ε (Mean = 263ε)) | Max = 10.5ε (Mean = 2.13ε) | Max = 10.5ε (Mean = 2.39ε) | Max = 4.6ε (Mean = 1.63ε)
Non Central T (large parameters) | Max = 478ε (Mean = 96.3ε) (Rmath 3.2.3: Max = 2.24ε (Mean = 0.945ε)) | Max = 9.79e+05ε (Mean = 1.97e+05ε) | Max = 9.79e+05ε (Mean = 1.97e+05ε) | Max = 227ε (Mean = 50.4ε)

Table 23.80. Error rates for owens_t

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Owens T (medium small values) | Max = 0ε (Mean = 0ε) | Max = 3.34ε (Mean = 0.944ε) | Max = 3.34ε (Mean = 0.911ε) | Max = 4.37ε (Mean = 0.98ε)
Owens T (large and diverse values) | Max = 0ε (Mean = 0ε) | Max = 49ε (Mean = 2.16ε) | Max = 24.5ε (Mean = 1.39ε) | Max = 3.78ε (Mean = 0.621ε)

Table 23.81. Error rates for polygamma

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Mathematica Data | Max = 0.824ε (Mean = 0.0574ε) (GSL 2.1: Max = 62.9ε (Mean = 12.8ε)) (Rmath 3.2.3: Max = 108ε (Mean = 15.2ε)) | Max = 7.38ε (Mean = 1.84ε) | Max = 34.3ε (Mean = 7.65ε) | Max = 9.32ε (Mean = 1.95ε)
Mathematica Data - large arguments | Max = 0.998ε (Mean = 0.0592ε) (GSL 2.1: Max = 244ε (Mean = 32.8ε) And other failures.) (Rmath 3.2.3: Max = 1.71e+56ε (Mean = 1.01e+55ε) And other failures.) | Max = 2.23ε (Mean = 0.323ε) | Max = 11.1ε (Mean = 0.848ε) | Max = 150ε (Mean = 13.9ε)
Mathematica Data - negative arguments | Max = 0.516ε (Mean = 0.022ε) (GSL 2.1: Max = 36.6ε (Mean = 3.04ε) And other failures.) (Rmath 3.2.3: Max = 0ε (Mean = 0ε) And other failures.) | Max = 269ε (Mean = 87.7ε) | Max = 269ε (Mean = 88.4ε) | Max = 497ε (Mean = 129ε)
Mathematica Data - large negative arguments | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.79ε (Mean = 0.197ε) And other failures.) (Rmath 3.2.3: Max = 0ε (Mean = 0ε) And other failures.) | Max = 155ε (Mean = 96.4ε) | Max = 155ε (Mean = 96.4ε) | Max = 162ε (Mean = 101ε)
Mathematica Data - small arguments | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 15.2ε (Mean = 5.03ε)) (Rmath 3.2.3: Max = 106ε (Mean = 20ε)) | Max = 3.33ε (Mean = 0.75ε) | Max = 3.33ε (Mean = 0.75ε) | Max = 3ε (Mean = 0.496ε)
Mathematica Data - Large orders and other bug cases | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 151ε (Mean = 39.3ε) And other failures.) (Rmath 3.2.3: Max = +INFε (Mean = +INFε) And other failures.) | Max = 54.5ε (Mean = 13.3ε) | Max = 145ε (Mean = 55.9ε) | Max = 200ε (Mean = 57.2ε)

Table 23.82. Error rates for powm1

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
powm1 | Max = 1.06ε (Mean = 0.425ε) | Max = 2.04ε (Mean = 0.493ε) | Max = 1.88ε (Mean = 0.49ε) | Max = 1.84ε (Mean = 0.486ε)

Table 23.83. Error rates for sin_pi

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
sin_pi and cos_pi | Max = 0ε (Mean = 0ε) | Max = 0.996ε (Mean = 0.335ε) | Max = 0.996ε (Mean = 0.336ε) | Max = 0.99ε (Mean = 0.328ε)
sin_pi and cos_pi near integers and half integers | Max = 0ε (Mean = 0ε) | Max = 0.976ε (Mean = 0.293ε) | Max = 0.976ε (Mean = 0.293ε) | Max = 0.996ε (Mean = 0.343ε)

Table 23.84. Error rates for sph_bessel

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel j: Random Data | Max = 243ε (Mean = 13.3ε) (<cmath>: Max = 1.91e+06ε (Mean = 1.09e+05ε)) | Max = 0.978ε (Mean = 0.0445ε) (GSL 2.1: Max = 1.79e+03ε (Mean = 107ε)) | Max = 243ε (Mean = 33.7ε) | Max = 245ε (Mean = 16.3ε)

Table 23.85. Error rates for sph_bessel_prime

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Bessel j': Random Data | Max = 0.753ε (Mean = 0.0343ε) | Max = 167ε (Mean = 12ε) | Max = 167ε (Mean = 33.2ε) | Max = 307ε (Mean = 25.2ε)

Table 23.86. Error rates for sph_neumann

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
y: Random Data | Max = 234ε (Mean = 19.5ε) (<cmath>: Max = 1.6e+06ε (Mean = 1.4e+05ε)) | Max = 0.995ε (Mean = 0.0665ε) (GSL 2.1: Max = 8.5e+04ε (Mean = 5.33e+03ε)) | Max = 234ε (Mean = 19.8ε) | Max = 281ε (Mean = 31.1ε)

Table 23.87. Error rates for sph_neumann_prime

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
y': Random Data | Max = 0.988ε (Mean = 0.0869ε) | Max = 158ε (Mean = 18.8ε) | Max = 158ε (Mean = 20.2ε) | Max = 296ε (Mean = 25.6ε)

Table 23.88. Error rates for spherical_harmonic_i

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Spherical Harmonics | Max = 1.36ε (Mean = 0.0765ε) | Max = 2.89e+03ε (Mean = 108ε) | Max = 1.03e+04ε (Mean = 327ε) | Max = 2.27e+04ε (Mean = 725ε)

Table 23.89. Error rates for spherical_harmonic_r

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Spherical Harmonics | Max = 1.58ε (Mean = 0.0707ε) | Max = 2.89e+03ε (Mean = 108ε) | Max = 1.03e+04ε (Mean = 327ε) | Max = 2.27e+04ε (Mean = 725ε)

Table 23.90. Error rates for sqrt1pm1

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
sqrt1pm1 | Max = 1.3ε (Mean = 0.404ε) | Max = 1.33ε (Mean = 0.404ε) | Max = 1.54ε (Mean = 0.563ε) | Max = 1.35ε (Mean = 0.497ε)

Table 23.91. Error rates for tgamma

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
factorials | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 3.95ε (Mean = 0.783ε)) (Rmath 3.2.3: Max = 314ε (Mean = 93.4ε)) | Max = 2.67ε (Mean = 0.617ε) (<cmath>: Max = 1.66ε (Mean = 0.584ε)) (<math.h>: Max = 1.66ε (Mean = 0.584ε)) | Max = 172ε (Mean = 41ε) (<math.h>: Max = 0ε (Mean = 0ε)) | Max = 1.85ε (Mean = 0.566ε) (<math.h>: Max = 3.17ε (Mean = 0.928ε))
near 0 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 4.51ε (Mean = 1.92ε)) (Rmath 3.2.3: Max = 1ε (Mean = 0.335ε)) | Max = 2ε (Mean = 0.608ε) (<cmath>: Max = 1ε (Mean = 0.376ε)) (<math.h>: Max = 1ε (Mean = 0.376ε)) | Max = 2ε (Mean = 0.647ε) (<math.h>: Max = 0.5ε (Mean = 0.0791ε)) | Max = 1.5ε (Mean = 0.635ε) (<math.h>: Max = 1ε (Mean = 0.405ε))
near 1 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 4.41ε (Mean = 1.81ε)) (Rmath 3.2.3: Max = 1ε (Mean = 0.32ε)) | Max = 2.51ε (Mean = 1.02ε) (<cmath>: Max = 0.918ε (Mean = 0.203ε)) (<math.h>: Max = 0.918ε (Mean = 0.203ε)) | Max = 3.01ε (Mean = 1.06ε) (<math.h>: Max = 1ε (Mean = 0.175ε)) | Max = 1.1ε (Mean = 0.59ε) (<math.h>: Max = 1ε (Mean = 0.4ε))
near 2 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 7.95ε (Mean = 3.12ε)) (Rmath 3.2.3: Max = 1ε (Mean = 0.191ε)) | Max = 4.1ε (Mean = 1.55ε) (<cmath>: Max = 0.558ε (Mean = 0.298ε)) (<math.h>: Max = 0.558ε (Mean = 0.298ε)) | Max = 5.01ε (Mean = 1.89ε) (<math.h>: Max = 0ε (Mean = 0ε)) | Max = 2ε (Mean = 0.733ε) (<math.h>: Max = 0ε (Mean = 0ε))
near -10 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 2.6ε (Mean = 1.05ε)) (Rmath 3.2.3: Max = 34.9ε (Mean = 9.2ε)) | Max = 1.75ε (Mean = 0.895ε) (<cmath>: Max = 2.26ε (Mean = 1.08ε)) (<math.h>: Max = 2.26ε (Mean = 1.08ε)) | Max = 1.75ε (Mean = 0.819ε) (<math.h>: Max = 0ε (Mean = 0ε)) | Max = 1.86ε (Mean = 0.881ε) (<math.h>: Max = 0.866ε (Mean = 0.445ε))
near -55 | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 1.8ε (Mean = 0.782ε)) (Rmath 3.2.3: Max = 3.89e+04ε (Mean = 9.52e+03ε)) | Max = 2.69ε (Mean = 1.09ε) (<cmath>: Max = 1.79ε (Mean = 0.75ε)) (<math.h>: Max = 1.79ε (Mean = 0.75ε)) | Max = 98.5ε (Mean = 53.4ε) (<math.h>: Max = 0ε (Mean = 0ε)) | Max = 2.7ε (Mean = 1.35ε) (<math.h>: Max = 3.87e+04ε (Mean = 6.71e+03ε))

Table 23.92. Error rates for tgamma1pm1

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
tgamma1pm1(dz) | Max = 0ε (Mean = 0ε) | Max = 1.12ε (Mean = 0.49ε) | Max = 6.61ε (Mean = 0.84ε) | Max = 3.31ε (Mean = 0.517ε)

Table 23.93. Error rates for tgamma_delta_ratio

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
tgamma + small delta ratios | Max = 0ε (Mean = 0ε) | Max = 5.83ε (Mean = 1.3ε) | Max = 15.4ε (Mean = 2.09ε) | Max = 7.56ε (Mean = 1.31ε)
tgamma + small delta ratios (negative delta) | Max = 0ε (Mean = 0ε) | Max = 7.94ε (Mean = 1.4ε) | Max = 18.3ε (Mean = 2.03ε) | Max = 7.43ε (Mean = 1.42ε)
tgamma + small integer ratios | Max = 0ε (Mean = 0ε) | Max = 1.96ε (Mean = 0.677ε) | Max = 1.96ε (Mean = 0.677ε) | Max = 2.74ε (Mean = 0.736ε)
tgamma + small integer ratios (negative delta) | Max = 0ε (Mean = 0ε) | Max = 1.62ε (Mean = 0.451ε) | Max = 1.62ε (Mean = 0.451ε) | Max = 2.15ε (Mean = 0.685ε)
integer tgamma ratios | Max = 0ε (Mean = 0ε) | Max = 0.997ε (Mean = 0.4ε) | Max = 0.997ε (Mean = 0.4ε) | Max = 0.968ε (Mean = 0.386ε)
integer tgamma ratios (negative delta) | Max = 0ε (Mean = 0ε) | Max = 0.853ε (Mean = 0.176ε) | Max = 0.853ε (Mean = 0.176ε) | Max = 0.974ε (Mean = 0.175ε)

Table 23.94. Error rates for tgamma (incomplete)

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
tgamma(a, z) medium values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 200ε (Mean = 13.3ε)) | Max = 8.47ε (Mean = 1.9ε) | Max = 412ε (Mean = 95.5ε) | Max = 8.14ε (Mean = 1.76ε)
tgamma(a, z) small values | Max = 0.753ε (Mean = 0.0474ε) (GSL 2.1: Max = 1.38e+10ε (Mean = 1.05e+09ε)) | Max = 2.31ε (Mean = 0.775ε) | Max = 2.13ε (Mean = 0.717ε) | Max = 2.53ε (Mean = 0.66ε)
tgamma(a, z) integer and half integer values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 117ε (Mean = 12.5ε)) | Max = 5.52ε (Mean = 1.48ε) | Max = 79.6ε (Mean = 20.9ε) | Max = 5.16ε (Mean = 1.33ε)

Table 23.95. Error rates for tgamma_lower

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
tgamma(a, z) medium values | Max = 0.833ε (Mean = 0.0315ε) (GSL 2.1: Max = 0.833ε (Mean = 0.0315ε)) | Max = 6.79ε (Mean = 1.46ε) | Max = 363ε (Mean = 63.8ε) | Max = 5.62ε (Mean = 1.49ε)
tgamma(a, z) small values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0ε (Mean = 0ε)) | Max = 1.97ε (Mean = 0.555ε) | Max = 1.97ε (Mean = 0.558ε) | Max = 1.57ε (Mean = 0.525ε)
tgamma(a, z) integer and half integer values | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0ε (Mean = 0ε)) | Max = 4.83ε (Mean = 1.15ε) | Max = 84.7ε (Mean = 17.5ε) | Max = 2.69ε (Mean = 0.849ε)

Table 23.96. Error rates for tgamma_ratio

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
tgamma ratios | Max = 0.694ε (Mean = 0.0347ε) | Max = 2.99ε (Mean = 1.15ε) | Max = 174ε (Mean = 61.2ε) | Max = 3.28ε (Mean = 1.12ε)

Table 23.97. Error rates for trigamma

| GNU C++ version 7.1.0 linux double | GNU C++ version 7.1.0 linux long double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Mathematica Data | Max = 0.998ε (Mean = 0.105ε) (GSL 2.1: Max = 1.34e+04ε (Mean = 1.49e+03ε)) (Rmath 3.2.3: Max = 1.34e+04ε (Mean = 1.51e+03ε)) | Max = 1.28ε (Mean = 0.449ε) | Max = 1.28ε (Mean = 0.449ε) | Max = 1ε (Mean = 0.382ε)

Table 23.98. Error rates for zeta

| GNU C++ version 7.1.0 linux long double | GNU C++ version 7.1.0 linux double | Sun compiler version 0x5150 Sun Solaris long double | Microsoft Visual C++ version 14.1 Win32 double
--- | --- | --- | --- | ---
Zeta: Random values greater than 1 | Max = 0.846ε (Mean = 0.0833ε) (<cmath>: Max = 5.45ε (Mean = 1ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 8.69ε (Mean = 1.03ε)) | Max = 0.846ε (Mean = 0.0833ε) | Max = 0.836ε (Mean = 0.093ε)
Zeta: Random values less than 1 | Max = 7.03ε (Mean = 2.93ε) (<cmath>: Max = 538ε (Mean = 59.3ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 137ε (Mean = 13.8ε)) | Max = 70.1ε (Mean = 17.1ε) | Max = 6.84ε (Mean = 3.12ε)
Zeta: Values close to and greater than 1 | Max = 0.995ε (Mean = 0.5ε) (<cmath>: Max = 1.9e+06ε (Mean = 5.11e+05ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 7.73ε (Mean = 4.07ε)) | Max = 0.995ε (Mean = 0.5ε) | Max = 0.994ε (Mean = 0.421ε)
Zeta: Values close to and less than 1 | Max = 0.998ε (Mean = 0.508ε) (<cmath>: Max = 8.53e+06ε (Mean = 1.87e+06ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 0.991ε (Mean = 0.28ε)) | Max = 0.998ε (Mean = 0.508ε) | Max = 0.991ε (Mean = 0.375ε)
Zeta: Integer arguments | Max = 9ε (Mean = 3.06ε) (<cmath>: Max = 70.3ε (Mean = 17.4ε)) | Max = 0ε (Mean = 0ε) (GSL 2.1: Max = 3.75ε (Mean = 1.1ε)) | Max = 28ε (Mean = 5.62ε) | Max = 9ε (Mean = 3ε)