F_COMPLEX cpow (F_TYPE a, F_TYPE b, F_TYPE c, F_TYPE d)
cpow(a,b,c,d) is defined as e**((c + id) ln(a + ib)).
Exceptional Argument Routine Behavior sqrt (a**2 + b**2) > max_float Overflow c/2 * ln(a**2 + b**2) > max_float Overflow c/2 * ln(a**2 + b**2) - (d * atan2(b,c)) Overflow > max_float