The casin functions compute the complex arc sine of z, with
    branch cuts outside the interval [-1, +1] along the real axis.