s01bac computes values of
, retaining full relative precision even when
is small. The function is based on the Chebyshev expansion
Setting
, and choosing
,
the expansion is valid in the domain
.
Background information to multithreading can be found in the
Multithreading documentation.