log1p
double log1p(double);
extern double log1p(double);
extern double log1p(double);
return (log1p(t + sqrt(2.0 * t + t * t)));
w = log1p(fabs(x) + t / (one + sqrt(one + t)));
return (copysign(0.5, x) * log1p(t + t));
return ((float)log1p((double)x));
D_IM(ans) = log1p(xm1 + sqrt(xm1 * (x + one)));
D_IM(ans) = log1p((x - one) + t);
D_IM(ans) = ln2 + log(y) + half * log1p(t * t);
D_IM(ans) = half * log1p((y + y) * (y + A));
D_IM(ans) = log1p(Am1 + sqrt(Am1 * (A + one)));
D_IM(ans) = log1p(xm1 + sqrt(xm1 * (x + one)));
D_IM(ans) = log1p((x - one) + sqrt((x - one) *
D_IM(ans) = ln2 + log(y) + half * log1p(t * t);
D_IM(ans) = half * log1p((y + y) * (y + A));
D_IM(ans) = log1p(Am1 + sqrt(Am1 * (A + one)));
D_IM(ans) = half * log1p(two / (-t));
D_IM(ans) = half * log1p((ay + ay) / t);
D_IM(ans) = 0.25 * log1p(t * t);
D_IM(ans) = half * log1p(two / (ay - one));
D_IM(ans) = half * log1p((ay + ay) / (one - ay));
D_IM(ans) = 0.25 * log1p((4.0 * ay) / (ax * ax + t * t));
F_IM(ans) = (float)(0.25 * log1p(4.0 * dy /
D_RE(ans) = half * log1p(ay * ay);
D_RE(ans) = half * (log1p(ay * ay + (ax -
D_RE(ans) = log(ax) + half * log1p(t * t);
F_RE(ans) = (float)(0.5 * log1p(dy * dy));
F_RE(ans) = (float)(0.5 * log1p(dx * dx));
dt = dhalf * log1p((dx - 1.0) * (dx + 1.0) + dy * dy);
dt = dhalf * log1p((dy - 1.0) * (dy + 1.0) + dx * dx);