bin/170206: complex arcsinh, log, etc.

[Bruce Evans]

On Sat, 28 Jul 2012, Stephen Montgomery-Smith wrote:

> OK.  This clog really seems to work.
>
> x*x + y*y - 1 is computed with a ULP less than 0.8.  The rest of the
> errors seem to be due to the implementation of log1p.  The ULP of the
> final answer seems to be never bigger than a little over 2.

I really don't like this version.  If we want bignums and slowness, we
can use lib[g]mp or even exec a mostly-interactive calculator (bc for
portability, closed-source bigware for unportability).

This version would have been useful in development to verify that the
doubled or tripled floating point was as exact as intended.  Now writing
it in your favourite interactive calculator is easiest, except for the
problem of comparing the results.

Bruce