The casinh functions compute the complex arc hyperbolic sine
    of z, with branch cuts outside the interval [-i, +i] along the
    imaginary axis.