qcomplex __cqpow(long double aqreal, long double aqimag, long double bqreal, long double bqimag) { long double logr, logi, x, y; qcomplex r, qc; logr = __qlog( __qhypot(aqreal, aqimag)); logi = __qatan2( aqimag, aqreal); x = __qexp( logr*bqreal - logi*bqimag); y = logr*bqimag + logi*bqreal; qc = __qcis(y); r.qreal = x*qc.qreal; r.qimag = x*qc.qimag; return r; }
long double __q_atn2(long double *y, long double *x) { return ( __qatan2(*y, *x) ); }