  We present an algorithm to approximate the real trilogarithm for a real
argument with IEEE 754-1985 double precision accuracy. The approximation is
structured such that it can make use of instruction-level parallelism when
executed on appropriate CPUs.
