log1pl
long double log1pl(long double);
extern long double log1pl(long double);
return (primary(x)-log1pl(x))-logl(x);
return (log1pl(t * (t + sqrtl(x + one))));
return (copysignl(log1pl(w + w / (t + sqrtl(one + t * t))), x));
#define LOG1P log1pl
extern long double log1pl(long double);
return (primary(x)-log1pl(x))-logl(x);
return (log1pl(t * (t + sqrtl(x + one))));
return (copysignl(log1pl(w + w / (t + sqrtl(one + t * t))), x));
return (copysignl(half, x) * log1pl(t + t));
extern long double log1pl(long double);
LD_IM(ans) = log1pl(xm1 + sqrtl(xm1 * (x +
LD_IM(ans) = log1pl(t - (one - x));
LD_IM(ans) = ln2 + logl(y) + half * log1pl(t * t);
LD_IM(ans) = half * log1pl((y + y) * (y + A));
LD_IM(ans) = log1pl(Am1 + sqrtl(Am1 * (A + one)));
LD_IM(ans) = log1pl(xm1 + sqrtl(xm1 * (x +
LD_IM(ans) = log1pl((x - one) + sqrtl((x -
LD_IM(ans) = ln2 + logl(y) + half * log1pl(t * t);
LD_IM(ans) = half * log1pl((y + y) * (y + A));
LD_IM(ans) = log1pl(Am1 + sqrtl(Am1 * (A + one)));
LD_IM(ans) = half * log1pl(two / (-t));
LD_IM(ans) = half * log1pl((ay + ay) / t);
LD_IM(ans) = 0.25L * log1pl(t * t);
LD_IM(ans) = half * log1pl(two / (ay - one));
LD_IM(ans) = half * log1pl((ay + ay) / (one - ay));
LD_IM(ans) = 0.25L * log1pl((4.0L * ay) / (ax * ax + t * t));
LD_RE(ans) = logl(ax) + half * log1pl(t * t);
LD_RE(ans) = half * log1pl(ay * ay);
LD_RE(ans) = half * (log1pl(ay * ay + (ax -