sqrtl
long double sqrtl(long double);
extern long double sqrtl(long double);
t = sqrtl(x - one);
return (log1pl(t * (t + sqrtl(x + one))));
return (copysignl(log1pl(w + w / (t + sqrtl(one + t * t))), x));
x = sqrtl(t1 * y1 - (w * (-w) - (t2 * y1 + y2 * x)));
x = sqrtl(t1 * t1 - (y * (-y) - t2 * (x + t1)));
return (invsqrtpi*cc)/sqrtl(x);
return (invsqrtpi*(u*cc-v*ss)/sqrtl(x));
return (invsqrtpi*ss)/sqrtl(x);
return (invsqrtpi*(pzero(x)*ss+qzero(x)*cc)/sqrtl(x));
return (invsqrtpi*cc)/sqrtl(x);
d = invsqrtpi*(pone(x)*cc-qone(x)*ss)/sqrtl(x);
return (invsqrtpi*ss)/sqrtl(x);
return (invsqrtpi*(pone(x)*ss+qone(x)*cc)/sqrtl(x));
b = invsqrtpi*temp/sqrtl(x);
b = invsqrtpi*temp/sqrtl(x);
extern long double sqrtl(long double);
t = sqrtl(x - one);
return (log1pl(t * (t + sqrtl(x + one))));
x = atanl(sqrtl((one - x) / (one + x)));
return (copysignl(log1pl(w + w / (t + sqrtl(one + t * t))), x));
return (atanl(x / sqrtl(one - x * x)));
return (atanl(x / sqrtl(w - t * t)));
x = sqrtl(t1 * t1 - (y * (-y) - t2 * (x + t1)));
x = sqrtl(t1 * y1 - (w * (-w) - (t2 * y1 + y2 * x)));
if (x>1.0e120L) return (invsqrtpi*cc)/sqrtl(x);
return invsqrtpi*(u*cc-v*ss)/sqrtl(x);
if (x>1.0e120L) return (invsqrtpi*ss)/sqrtl(x);
return invsqrtpi*(pzero(x)*ss+qzero(x)*cc)/sqrtl(x);
if (x>1.0e120L) return (invsqrtpi*cc)/sqrtl(x);
d = invsqrtpi*(pone(x)*cc-qone(x)*ss)/sqrtl(x);
if (x>1.0e91L) return (invsqrtpi*ss)/sqrtl(x);
return invsqrtpi*(pone(x)*ss+qone(x)*cc)/sqrtl(x);
b = invsqrtpi*temp/sqrtl(x);
b = invsqrtpi*temp/sqrtl(x);
extern long double sqrtl(long double);
return (sqrtl(x));
LD_IM(ans) = logl(x + sqrtl((x - one) * (x +
LD_IM(ans) = log1pl(xm1 + sqrtl(xm1 * (x +
LD_IM(ans) = y / sqrtl((one + x) * (one - x));
t = sqrtl((x - one) * (x + one));
t = sqrtl(y);
A = sqrtl(one + y * y);
R = sqrtl(xp1 * xp1 + y2);
S = sqrtl(xm1 * xm1 + y2);
LD_RE(ans) = atan2l(sqrtl(half * Apx * (y2 /
LD_RE(ans) = atan2l((y * sqrtl(half * (Apx /
LD_IM(ans) = log1pl(Am1 + sqrtl(Am1 * (A + one)));
LD_IM(ans) = logl(A + sqrtl(A * A - one));
LD_IM(ans) = logl(x + sqrtl((x - one) * (x +
LD_IM(ans) = log1pl(xm1 + sqrtl(xm1 * (x +
LD_IM(ans) = y / sqrtl((one + x) * (one - x));
LD_IM(ans) = logl(x + sqrtl((x - one) * (x +
LD_IM(ans) = log1pl((x - one) + sqrtl((x -
t = sqrtl(y);
A = sqrtl(one + y * y);
R = sqrtl(xp1 * xp1 + y2);
S = sqrtl(xm1 * xm1 + y2);
LD_RE(ans) = atanl(x / sqrtl(half * Apx * (y2 /
LD_RE(ans) = atanl(x / (y * sqrtl(half * (Apx /
LD_IM(ans) = log1pl(Am1 + sqrtl(Am1 * (A + one)));
LD_IM(ans) = logl(A + sqrtl(A * A - one));
t = sqrtl(half * (ax + sqrtl(ax * ax + y * y)));
t = sqrtl(half * ay);
t = sqrtl(half * ay + half * ax);
t = half * (sqrtl(two * (ax + ay)));
t = sqrtl(half * (ax + ay));
t = two4500 * sqrtl(ax + sqrtl(ax * ax + y * y));
t = twom4500 * sqrtl(ax + sqrtl(ax * ax + y * y));
t = sqrtl(half * (ax + sqrtl(ax * ax + y * y)));
LD_RE(ans) = sqrtl(ax);
LD_IM(ans) = sqrtl(ax);
t = sqrtl(ax);
t = two4500 * sqrtl(ax + sqrtl(ax * ax + y * y));
t = twom4500 * sqrtl(ax + sqrtl(ax * ax + y * y));