__complex128 catanhq (__complex128 x) { __complex128 res; int rcls = fpclassifyq (__real__ x); int icls = fpclassifyq (__imag__ x); if (rcls <= QUADFP_INFINITE || icls <= QUADFP_INFINITE) { if (icls == QUADFP_INFINITE) { __real__ res = copysignq (0.0, __real__ x); __imag__ res = copysignq (M_PI_2q, __imag__ x); } else if (rcls == QUADFP_INFINITE || rcls == QUADFP_ZERO) { __real__ res = copysignq (0.0, __real__ x); if (icls >= QUADFP_ZERO) __imag__ res = copysignq (M_PI_2q, __imag__ x); else __imag__ res = nanq (""); } else { __real__ res = nanq (""); __imag__ res = nanq (""); } } else if (rcls == QUADFP_ZERO && icls == QUADFP_ZERO) { res = x; } else { __float128 i2, num, den; i2 = __imag__ x * __imag__ x; num = 1.0 + __real__ x; num = i2 + num * num; den = 1.0 - __real__ x; den = i2 + den * den; __real__ res = 0.25 * (logq (num) - logq (den)); den = 1 - __real__ x * __real__ x - i2; __imag__ res = 0.5 * atan2q (2.0 * __imag__ x, den); } return res; }
__complex128 clogq (__complex128 z) { __complex128 v; COMPLEX_ASSIGN (v, logq (cabsq (z)), cargq (z)); return v; }
__float128 logxq (__float128 num, __float128 base) { return logq (num) / logq (base); }
inline void eval_log(float128_backend& result, const float128_backend& arg) { result.value() = logq(arg.value()); }
inline Quad Log( const Quad& alpha ) { return logq(alpha); }
tFloat p_u(tFloat s1, tFloat s2, tFloat s3, tFloat h1, tFloat h2, tFloat h3) { return logq((s2-s3)/(s1-s2))/logq(h3/h2); }