sqrtf
*rx = sqrtf(x);
*rx = log1pf(Am1 + sqrtf(Am1 * (A + 1)));
*rx = x / sqrtf((1 - y) * (1 + y));
*rx = log1pf((y - 1) + sqrtf((y - 1) * (y + 1)));
*rx = logf(A + sqrtf(A * A - 1));
*sqrt_A2my2 = sqrtf(x) * sqrtf((A + y) / 2);
*sqrt_A2my2 = sqrtf(Amy * (A + y));
sqrtf((y + 1) * (y - 1));
*sqrt_A2my2 = sqrtf((1 - y) * (1 + y));
s = sqrtf(z);
s = sqrtf(z);
return logf((float)2.0*x-one/(x+sqrtf(t-one)));
return log1pf(t+sqrtf((float)2.0*t+t*t));
w = sqrtf(t1*t1-(b*(-b)-t2*(a+t1)));
w = sqrtf(t1*y1-(w*(-w)-(t1*y2+t2*b)));
if(ix>0x58000000) z = (invsqrtpi*ss)/sqrtf(x); /* |x|>2**49 */
z = invsqrtpi*(u*ss+v*cc)/sqrtf(x);
if(ix>0x58000000) z = (invsqrtpi*cc)/sqrtf(x); /* |x|>2**49 */
z = invsqrtpi*(u*cc-v*ss)/sqrtf(x);
if(ix>0x58000000) z = (invsqrtpi*ss)/sqrtf(x); /* |x|>2**49 */
z = invsqrtpi*(u*ss+v*cc)/sqrtf(x);
if(ix>0x58000000) z = (invsqrtpi*cc)/sqrtf(y); /* |x|>2**49 */
z = invsqrtpi*(u*cc-v*ss)/sqrtf(y);
return sqrtf(x);
float sqrtf(float);
w = logf((float)2.0*t+one/(sqrtf(x*x+one)+t));
w =log1pf(fabsf(x)+t/(one+sqrtf(one+t)));
rho = sqrtf(1 + s * s);