Lines Matching refs:fp0
233 floatx80 fp0, fp1, fp2, fp3, f, logof2, klog2, saveu; in floatx80_lognp1() local
276 fp0 = a; /* Z */ in floatx80_lognp1()
279 fp0 = floatx80_add(fp0, float32_to_floatx80(make_float32(0x3F800000), in floatx80_lognp1()
282 aExp = extractFloatx80Exp(fp0); in floatx80_lognp1()
283 aSig = extractFloatx80Frac(fp0); in floatx80_lognp1()
296 fp0 = packFloatx80(0, 0x3FFF, aSig); /* Y */ in floatx80_lognp1()
298 fp0 = floatx80_sub(fp0, f, status); /* Y-F */ in floatx80_lognp1()
302 fp0 = floatx80_mul(fp0, log_tbl[j], status); /* FP0 IS U = (Y-F)/F */ in floatx80_lognp1()
305 fp2 = floatx80_mul(fp0, fp0, status); /* FP2 IS V=U*U */ in floatx80_lognp1()
332 fp1 = floatx80_mul(fp1, fp0, status); /* U*V*(A2+V*(A4+V*A6)) */ in floatx80_lognp1()
333 fp0 = floatx80_add(fp0, fp2, status); /* U+V*(A1+V*(A3+V*A5)) */ in floatx80_lognp1()
337 fp0 = floatx80_add(fp0, fp1, status); /* FP0 IS LOG(F) + LOG(1+U) */ in floatx80_lognp1()
342 a = floatx80_add(fp0, klog2, status); in floatx80_lognp1()
356 fp0 = floatx80_sub(float32_to_floatx80(make_float32(0x3F800000), in floatx80_lognp1()
358 fp0 = floatx80_add(fp0, fp1, status); /* FP0 IS Y-F = (1-F)+Z */ in floatx80_lognp1()
362 fp0 = floatx80_sub(float32_to_floatx80(make_float32(0x40000000), in floatx80_lognp1()
365 fp0 = floatx80_add(fp0, fp1, status); /* FP0 IS Y-F = (2-F)+2Z */ in floatx80_lognp1()
372 fp0 = floatx80_add(fp0, float32_to_floatx80(make_float32(0x3F800000), in floatx80_lognp1()
376 fp1 = floatx80_div(fp1, fp0, status); /* U */ in floatx80_lognp1()
378 fp0 = floatx80_mul(fp1, fp1, status); /* FP0 IS V = U*U */ in floatx80_lognp1()
379 fp1 = floatx80_mul(fp0, fp0, status); /* FP1 IS W = V*V */ in floatx80_lognp1()
394 fp2 = floatx80_mul(fp2, fp0, status); /* V*(B2+W*B4) */ in floatx80_lognp1()
399 fp0 = floatx80_mul(fp0, saveu, status); /* FP0 IS U*V */ in floatx80_lognp1()
402 fp0 = floatx80_mul(fp0, fp1, in floatx80_lognp1()
408 a = floatx80_add(fp0, saveu, status); in floatx80_lognp1()
432 floatx80 fp0, fp1, fp2, fp3, f, logof2, klog2, saveu; in floatx80_logn() local
484 fp0 = packFloatx80(0, 0x3FFF, aSig); /* Y */ in floatx80_logn()
486 fp0 = floatx80_sub(fp0, f, status); /* Y-F */ in floatx80_logn()
489 fp0 = floatx80_mul(fp0, log_tbl[j], status); /* FP0 IS U = (Y-F)/F */ in floatx80_logn()
492 fp2 = floatx80_mul(fp0, fp0, status); /* FP2 IS V=U*U */ in floatx80_logn()
519 fp1 = floatx80_mul(fp1, fp0, status); /* U*V*(A2+V*(A4+V*A6)) */ in floatx80_logn()
520 fp0 = floatx80_add(fp0, fp2, status); /* U+V*(A1+V*(A3+V*A5)) */ in floatx80_logn()
524 fp0 = floatx80_add(fp0, fp1, status); /* FP0 IS LOG(F) + LOG(1+U) */ in floatx80_logn()
529 a = floatx80_add(fp0, klog2, status); in floatx80_logn()
535 fp0 = a; in floatx80_logn()
539 fp0 = floatx80_add(fp0, float32_to_floatx80(make_float32(0x3F800000), in floatx80_logn()
544 fp1 = floatx80_div(fp1, fp0, status); /* U */ in floatx80_logn()
546 fp0 = floatx80_mul(fp1, fp1, status); /* FP0 IS V = U*U */ in floatx80_logn()
547 fp1 = floatx80_mul(fp0, fp0, status); /* FP1 IS W = V*V */ in floatx80_logn()
562 fp2 = floatx80_mul(fp2, fp0, status); /* V*(B2+W*B4) */ in floatx80_logn()
567 fp0 = floatx80_mul(fp0, saveu, status); /* FP0 IS U*V */ in floatx80_logn()
569 fp0 = floatx80_mul(fp0, fp1, in floatx80_logn()
575 a = floatx80_add(fp0, saveu, status); in floatx80_logn()
598 floatx80 fp0, fp1; in floatx80_log10() local
628 fp0 = floatx80_logn(a, status); in floatx80_log10()
634 a = floatx80_mul(fp0, fp1, status); /* LOGN(X)*INV_L10 */ in floatx80_log10()
654 floatx80 fp0, fp1; in floatx80_log2() local
693 fp0 = floatx80_logn(a, status); in floatx80_log2()
699 a = floatx80_mul(fp0, fp1, status); /* LOGN(X)*INV_L2 */ in floatx80_log2()
721 floatx80 fp0, fp1, fp2, fp3, l2, scale, adjscale; in floatx80_etox() local
753 fp0 = a; in floatx80_etox()
755 fp0 = floatx80_mul(fp0, float32_to_floatx80( in floatx80_etox()
759 n = floatx80_to_int32(fp0, status); /* int(64/log2*X) */ in floatx80_etox()
760 fp0 = int32_to_floatx80(n, status); in floatx80_etox()
774 fp2 = fp0; /* N */ in floatx80_etox()
775 fp0 = floatx80_mul(fp0, float32_to_floatx80( in floatx80_etox()
780 fp0 = floatx80_add(fp0, fp1, status); /* X + N*L1 */ in floatx80_etox()
781 fp0 = floatx80_add(fp0, fp2, status); /* R */ in floatx80_etox()
783 fp1 = floatx80_mul(fp0, fp0, status); /* S = R*R */ in floatx80_etox()
801 fp3 = floatx80_mul(fp3, fp0, status); /* fp3 IS R*S*(A2+S*A4) */ in floatx80_etox()
804 fp0 = floatx80_add(fp0, fp3, status); /* fp0 IS R+R*S*(A2+S*A4) */ in floatx80_etox()
805 fp0 = floatx80_add(fp0, fp2, status); /* fp0 IS EXP(R) - 1 */ in floatx80_etox()
808 fp0 = floatx80_mul(fp0, fp1, status); /* 2^(J/64)*(Exp(R)-1) */ in floatx80_etox()
809 fp0 = floatx80_add(fp0, float32_to_floatx80(exp_tbl2[j], status), in floatx80_etox()
811 fp0 = floatx80_add(fp0, fp1, in floatx80_etox()
817 fp0 = floatx80_mul(fp0, adjscale, status); in floatx80_etox()
823 a = floatx80_mul(fp0, scale, status); in floatx80_etox()
845 fp0 = a; in floatx80_etox()
847 fp0 = floatx80_mul(fp0, float32_to_floatx80( in floatx80_etox()
851 n = floatx80_to_int32(fp0, status); /* int(64/log2*X) */ in floatx80_etox()
852 fp0 = int32_to_floatx80(n, status); in floatx80_etox()
905 floatx80 fp0, fp1, fp2, fp3, adjfact, fact1, fact2; in floatx80_twotox() local
930 fp0 = a; in floatx80_twotox()
951 a = floatx80_add(fp0, float32_to_floatx80( in floatx80_twotox()
960 fp1 = fp0; /* X */ in floatx80_twotox()
997 fp0 = floatx80_sub(fp0, fp1, status); /* X - (1/64)*INT(64 X) */ in floatx80_twotox()
999 fp0 = floatx80_mul(fp0, fp2, status); /* R */ in floatx80_twotox()
1002 fp1 = floatx80_mul(fp0, fp0, status); /* S = R*R */ in floatx80_twotox()
1020 fp3 = floatx80_mul(fp3, fp0, status); /* R*S*(A2+S*A4) */ in floatx80_twotox()
1023 fp0 = floatx80_add(fp0, fp3, status); /* R+R*S*(A2+S*A4) */ in floatx80_twotox()
1024 fp0 = floatx80_add(fp0, fp2, status); /* EXP(R) - 1 */ in floatx80_twotox()
1026 fp0 = floatx80_mul(fp0, fact1, status); in floatx80_twotox()
1027 fp0 = floatx80_add(fp0, fact2, status); in floatx80_twotox()
1028 fp0 = floatx80_add(fp0, fact1, status); in floatx80_twotox()
1033 a = floatx80_mul(fp0, adjfact, status); in floatx80_twotox()
1055 floatx80 fp0, fp1, fp2, fp3, adjfact, fact1, fact2; in floatx80_tentox() local
1080 fp0 = a; in floatx80_tentox()
1101 a = floatx80_add(fp0, float32_to_floatx80( in floatx80_tentox()
1110 fp1 = fp0; /* X */ in floatx80_tentox()
1151 fp0 = floatx80_sub(fp0, fp1, status); /* X - N L_LEAD */ in floatx80_tentox()
1152 fp0 = floatx80_sub(fp0, fp2, status); /* X - N L_TRAIL */ in floatx80_tentox()
1154 fp0 = floatx80_mul(fp0, fp2, status); /* R */ in floatx80_tentox()
1157 fp1 = floatx80_mul(fp0, fp0, status); /* S = R*R */ in floatx80_tentox()
1175 fp3 = floatx80_mul(fp3, fp0, status); /* R*S*(A2+S*A4) */ in floatx80_tentox()
1178 fp0 = floatx80_add(fp0, fp3, status); /* R+R*S*(A2+S*A4) */ in floatx80_tentox()
1179 fp0 = floatx80_add(fp0, fp2, status); /* EXP(R) - 1 */ in floatx80_tentox()
1181 fp0 = floatx80_mul(fp0, fact1, status); in floatx80_tentox()
1182 fp0 = floatx80_add(fp0, fact2, status); in floatx80_tentox()
1183 fp0 = floatx80_add(fp0, fact1, status); in floatx80_tentox()
1188 a = floatx80_mul(fp0, adjfact, status); in floatx80_tentox()
1210 floatx80 fp0, fp1, fp2, fp3, fp4, fp5, invtwopi, twopi1, twopi2; in floatx80_tan() local
1237 fp0 = a; in floatx80_tan()
1249 fp0 = floatx80_add(fp0, twopi1, status); in floatx80_tan()
1250 fp1 = fp0; in floatx80_tan()
1251 fp0 = floatx80_add(fp0, twopi2, status); in floatx80_tan()
1252 fp1 = floatx80_sub(fp1, fp0, status); in floatx80_tan()
1256 xSign = extractFloatx80Sign(fp0); in floatx80_tan()
1257 xExp = extractFloatx80Exp(fp0); in floatx80_tan()
1274 fp2 = floatx80_mul(fp0, invtwopi, status); in floatx80_tan()
1283 fp0 = floatx80_sub(fp0, fp3, status); /* FP0 is A := R - P */ in floatx80_tan()
1285 fp3 = fp0; /* FP3 is A */ in floatx80_tan()
1287 fp0 = floatx80_add(fp0, fp1, status); /* FP0 is R := A+a */ in floatx80_tan()
1293 fp3 = floatx80_sub(fp3, fp0, status); /* A-R */ in floatx80_tan()
1307 fp1 = floatx80_mul(fp0, float64_to_floatx80( in floatx80_tan()
1314 fp0 = floatx80_sub(fp0, pi_tbl[j], status); /* X-Y1 */ in floatx80_tan()
1315 fp0 = floatx80_sub(fp0, float32_to_floatx80(pi_tbl2[j], status), in floatx80_tan()
1321 fp1 = fp0; /* R */ in floatx80_tan()
1322 fp0 = floatx80_mul(fp0, fp0, status); /* S = R*R */ in floatx80_tan()
1327 fp3 = floatx80_mul(fp3, fp0, status); /* SQ4 */ in floatx80_tan()
1328 fp2 = floatx80_mul(fp2, fp0, status); /* SP3 */ in floatx80_tan()
1334 fp3 = floatx80_mul(fp3, fp0, status); /* S(Q3+SQ4) */ in floatx80_tan()
1335 fp2 = floatx80_mul(fp2, fp0, status); /* S(P2+SP3) */ in floatx80_tan()
1340 fp3 = floatx80_mul(fp3, fp0, status); /* S(Q2+S(Q3+SQ4)) */ in floatx80_tan()
1341 fp2 = floatx80_mul(fp2, fp0, status); /* S(P1+S(P2+SP3)) */ in floatx80_tan()
1345 fp0 = floatx80_mul(fp0, fp3, status); /* S(Q1+S(Q2+S(Q3+SQ4))) */ in floatx80_tan()
1347 fp0 = floatx80_add(fp0, float32_to_floatx80( in floatx80_tan()
1360 a = floatx80_div(fp0, fp1, status); in floatx80_tan()
1366 fp1 = floatx80_mul(fp0, fp0, status); /* S = R*R */ in floatx80_tan()
1388 fp2 = floatx80_mul(fp2, fp0, status); /* RS(P1+S(P2+SP3)) */ in floatx80_tan()
1390 fp0 = floatx80_add(fp0, fp2, status); /* R+RS(P1+S(P2+SP3)) */ in floatx80_tan()
1398 a = floatx80_div(fp0, fp1, status); in floatx80_tan()
1421 floatx80 fp0, fp1, fp2, fp3, fp4, fp5, x, invtwopi, twopi1, twopi2; in floatx80_sin() local
1448 fp0 = a; in floatx80_sin()
1460 fp0 = floatx80_add(fp0, twopi1, status); in floatx80_sin()
1461 fp1 = fp0; in floatx80_sin()
1462 fp0 = floatx80_add(fp0, twopi2, status); in floatx80_sin()
1463 fp1 = floatx80_sub(fp1, fp0, status); in floatx80_sin()
1467 xSign = extractFloatx80Sign(fp0); in floatx80_sin()
1468 xExp = extractFloatx80Exp(fp0); in floatx80_sin()
1485 fp2 = floatx80_mul(fp0, invtwopi, status); in floatx80_sin()
1494 fp0 = floatx80_sub(fp0, fp3, status); /* FP0 is A := R - P */ in floatx80_sin()
1496 fp3 = fp0; /* FP3 is A */ in floatx80_sin()
1498 fp0 = floatx80_add(fp0, fp1, status); /* FP0 is R := A+a */ in floatx80_sin()
1504 fp3 = floatx80_sub(fp3, fp0, status); /* A-R */ in floatx80_sin()
1509 fp0 = float32_to_floatx80(make_float32(0x3F800000), in floatx80_sin()
1522 fp1 = floatx80_mul(fp0, float64_to_floatx80( in floatx80_sin()
1529 fp0 = floatx80_sub(fp0, pi_tbl[j], status); /* X-Y1 */ in floatx80_sin()
1530 fp0 = floatx80_sub(fp0, float32_to_floatx80(pi_tbl2[j], status), in floatx80_sin()
1536 fp0 = floatx80_mul(fp0, fp0, status); /* FP0 IS S */ in floatx80_sin()
1537 fp1 = floatx80_mul(fp0, fp0, status); /* FP1 IS T */ in floatx80_sin()
1543 xSign = extractFloatx80Sign(fp0); /* X IS S */ in floatx80_sin()
1544 xExp = extractFloatx80Exp(fp0); in floatx80_sin()
1545 xSig = extractFloatx80Frac(fp0); in floatx80_sin()
1577 fp0 = floatx80_mul(fp0, fp2, status); /* S(B2+T(B4+T(B6+TB8))) */ in floatx80_sin()
1578 fp0 = floatx80_add(fp0, fp1, status); /* [B1+T(B3+T(B5+TB7))]+ in floatx80_sin()
1583 fp0 = floatx80_mul(fp0, x, status); in floatx80_sin()
1588 a = floatx80_add(fp0, float32_to_floatx80(posneg1, status), status); in floatx80_sin()
1595 xSign = extractFloatx80Sign(fp0); /* X IS R */ in floatx80_sin()
1596 xExp = extractFloatx80Exp(fp0); in floatx80_sin()
1597 xSig = extractFloatx80Frac(fp0); in floatx80_sin()
1601 fp0 = floatx80_mul(fp0, fp0, status); /* FP0 IS S */ in floatx80_sin()
1602 fp1 = floatx80_mul(fp0, fp0, status); /* FP1 IS T */ in floatx80_sin()
1623 fp2 = floatx80_mul(fp2, fp0, status); /* S(A2+T(A4+TA6)) */ in floatx80_sin()
1632 fp0 = floatx80_mul(fp0, x, status); /* R'*S */ in floatx80_sin()
1633 fp0 = floatx80_mul(fp0, fp1, status); /* SIN(R')-R' */ in floatx80_sin()
1638 a = floatx80_add(fp0, x, status); in floatx80_sin()
1661 floatx80 fp0, fp1, fp2, fp3, fp4, fp5, x, invtwopi, twopi1, twopi2; in floatx80_cos() local
1688 fp0 = a; in floatx80_cos()
1700 fp0 = floatx80_add(fp0, twopi1, status); in floatx80_cos()
1701 fp1 = fp0; in floatx80_cos()
1702 fp0 = floatx80_add(fp0, twopi2, status); in floatx80_cos()
1703 fp1 = floatx80_sub(fp1, fp0, status); in floatx80_cos()
1707 xSign = extractFloatx80Sign(fp0); in floatx80_cos()
1708 xExp = extractFloatx80Exp(fp0); in floatx80_cos()
1725 fp2 = floatx80_mul(fp0, invtwopi, status); in floatx80_cos()
1734 fp0 = floatx80_sub(fp0, fp3, status); /* FP0 is A := R - P */ in floatx80_cos()
1736 fp3 = fp0; /* FP3 is A */ in floatx80_cos()
1738 fp0 = floatx80_add(fp0, fp1, status); /* FP0 is R := A+a */ in floatx80_cos()
1744 fp3 = floatx80_sub(fp3, fp0, status); /* A-R */ in floatx80_cos()
1749 fp0 = float32_to_floatx80(make_float32(0x3F800000), status); /* 1 */ in floatx80_cos()
1755 a = floatx80_sub(fp0, float32_to_floatx80( in floatx80_cos()
1763 fp1 = floatx80_mul(fp0, float64_to_floatx80( in floatx80_cos()
1770 fp0 = floatx80_sub(fp0, pi_tbl[j], status); /* X-Y1 */ in floatx80_cos()
1771 fp0 = floatx80_sub(fp0, float32_to_floatx80(pi_tbl2[j], status), in floatx80_cos()
1777 fp0 = floatx80_mul(fp0, fp0, status); /* FP0 IS S */ in floatx80_cos()
1778 fp1 = floatx80_mul(fp0, fp0, status); /* FP1 IS T */ in floatx80_cos()
1784 xSign = extractFloatx80Sign(fp0); /* X IS S */ in floatx80_cos()
1785 xExp = extractFloatx80Exp(fp0); in floatx80_cos()
1786 xSig = extractFloatx80Frac(fp0); in floatx80_cos()
1818 fp0 = floatx80_mul(fp0, fp2, status); /* S(B2+T(B4+T(B6+TB8))) */ in floatx80_cos()
1819 fp0 = floatx80_add(fp0, fp1, status); in floatx80_cos()
1823 fp0 = floatx80_mul(fp0, x, status); in floatx80_cos()
1828 a = floatx80_add(fp0, float32_to_floatx80(posneg1, status), status); in floatx80_cos()
1835 xSign = extractFloatx80Sign(fp0); /* X IS R */ in floatx80_cos()
1836 xExp = extractFloatx80Exp(fp0); in floatx80_cos()
1837 xSig = extractFloatx80Frac(fp0); in floatx80_cos()
1841 fp0 = floatx80_mul(fp0, fp0, status); /* FP0 IS S */ in floatx80_cos()
1842 fp1 = floatx80_mul(fp0, fp0, status); /* FP1 IS T */ in floatx80_cos()
1863 fp2 = floatx80_mul(fp2, fp0, status); /* S(A2+T(A4+TA6)) */ in floatx80_cos()
1870 fp0 = floatx80_mul(fp0, x, status); /* R'*S */ in floatx80_cos()
1871 fp0 = floatx80_mul(fp0, fp1, status); /* SIN(R')-R' */ in floatx80_cos()
1876 a = floatx80_add(fp0, x, status); in floatx80_cos()
1899 floatx80 fp0, fp1, fp2, fp3, xsave; in floatx80_atan() local
1929 fp0 = packFloatx80(aSign, piby2_exp, pi_sig); in floatx80_atan()
1935 a = floatx80_sub(fp0, fp1, status); in floatx80_atan()
1941 fp0 = a; in floatx80_atan()
1943 fp1 = floatx80_div(fp1, fp0, status); /* X' = -1/X */ in floatx80_atan()
1945 fp0 = floatx80_mul(fp1, fp1, status); /* Y = X'*X' */ in floatx80_atan()
1946 fp1 = floatx80_mul(fp0, fp0, status); /* Z = Y*Y */ in floatx80_atan()
1960 fp2 = floatx80_mul(fp2, fp0, status); /* Y*(C2+Z*C4) */ in floatx80_atan()
1964 fp0 = floatx80_mul(fp0, xsave, status); /* X'*Y */ in floatx80_atan()
1968 fp0 = floatx80_mul(fp0, fp1, status); in floatx80_atan()
1969 fp0 = floatx80_add(fp0, xsave, status); in floatx80_atan()
1975 a = floatx80_add(fp0, fp1, status); in floatx80_atan()
1992 fp0 = a; in floatx80_atan()
1994 fp0 = floatx80_mul(fp0, fp0, status); /* Y = X*X */ in floatx80_atan()
1995 fp1 = floatx80_mul(fp0, fp0, status); /* Z = Y*Y */ in floatx80_atan()
2016 fp2 = floatx80_mul(fp2, fp0, status); /* Y*(B2+Z*(B4+Z*B6)) */ in floatx80_atan()
2017 fp0 = floatx80_mul(fp0, xsave, status); /* X*Y */ in floatx80_atan()
2021 fp0 = floatx80_mul(fp0, fp1, status); in floatx80_atan()
2026 a = floatx80_add(fp0, xsave, status); in floatx80_atan()
2037 fp0 = a; in floatx80_atan()
2041 fp0 = floatx80_sub(fp0, xsave, status); /* X-F */ in floatx80_atan()
2043 fp0 = floatx80_div(fp0, fp1, status); /* U = (X-F)/(1+X*F) */ in floatx80_atan()
2057 fp1 = floatx80_mul(fp0, fp0, status); /* V = U*U */ in floatx80_atan()
2062 fp1 = floatx80_mul(fp1, fp0, status); /* U*V */ in floatx80_atan()
2070 fp0 = floatx80_add(fp0, fp1, status); /* ATAN(U) */ in floatx80_atan()
2075 a = floatx80_add(fp0, fp3, status); /* ATAN(X) */ in floatx80_atan()
2097 floatx80 fp0, fp1, fp2, one; in floatx80_asin() local
2131 fp0 = a; in floatx80_asin()
2133 fp1 = floatx80_sub(one, fp0, status); /* 1 - X */ in floatx80_asin()
2134 fp2 = floatx80_add(one, fp0, status); /* 1 + X */ in floatx80_asin()
2137 fp0 = floatx80_div(fp0, fp1, status); /* X/SQRT((1+X)*(1-X)) */ in floatx80_asin()
2142 a = floatx80_atan(fp0, status); /* ATAN(X/SQRT((1+X)*(1-X))) */ in floatx80_asin()
2163 floatx80 fp0, fp1, one; in floatx80_acos() local
2201 fp0 = a; in floatx80_acos()
2203 fp1 = floatx80_add(one, fp0, status); /* 1 + X */ in floatx80_acos()
2204 fp0 = floatx80_sub(one, fp0, status); /* 1 - X */ in floatx80_acos()
2205 fp0 = floatx80_div(fp0, fp1, status); /* (1-X)/(1+X) */ in floatx80_acos()
2206 fp0 = floatx80_sqrt(fp0, status); /* SQRT((1-X)/(1+X)) */ in floatx80_acos()
2207 fp0 = floatx80_atan(fp0, status); /* ATAN(SQRT((1-X)/(1+X))) */ in floatx80_acos()
2212 a = floatx80_add(fp0, fp0, status); /* 2 * ATAN(SQRT((1-X)/(1+X))) */ in floatx80_acos()
2233 floatx80 fp0, fp1, fp2, one; in floatx80_atanh() local
2266 fp0 = packFloatx80(0, aExp, aSig); /* Y = |X| */ in floatx80_atanh()
2268 fp0 = floatx80_add(fp0, fp0, status); /* 2Y */ in floatx80_atanh()
2270 fp0 = floatx80_div(fp0, fp1, status); /* Z = 2Y/(1-Y) */ in floatx80_atanh()
2271 fp0 = floatx80_lognp1(fp0, status); /* LOG1P(Z) */ in floatx80_atanh()
2276 a = floatx80_mul(fp0, fp2, in floatx80_atanh()
2298 floatx80 fp0, fp1, fp2, fp3, l2, sc, onebysc; in floatx80_etoxm1() local
2327 fp0 = a; in floatx80_etoxm1()
2329 fp0 = floatx80_mul(fp0, float32_to_floatx80( in floatx80_etoxm1()
2332 n = floatx80_to_int32(fp0, status); /* int(64/log2*X) */ in floatx80_etoxm1()
2333 fp0 = int32_to_floatx80(n, status); in floatx80_etoxm1()
2348 fp2 = fp0; /* N */ in floatx80_etoxm1()
2349 fp0 = floatx80_mul(fp0, float32_to_floatx80( in floatx80_etoxm1()
2354 fp0 = floatx80_add(fp0, fp1, status); /* X + N*L1 */ in floatx80_etoxm1()
2355 fp0 = floatx80_add(fp0, fp2, status); /* R */ in floatx80_etoxm1()
2357 fp1 = floatx80_mul(fp0, fp0, status); /* S = R*R */ in floatx80_etoxm1()
2381 fp2 = floatx80_mul(fp2, fp0, in floatx80_etoxm1()
2383 fp0 = floatx80_add(fp0, fp1, in floatx80_etoxm1()
2385 fp0 = floatx80_add(fp0, fp2, status); /* fp0 IS EXP(R) - 1 */ in floatx80_etoxm1()
2387 fp0 = floatx80_mul(fp0, exp_tbl[j], in floatx80_etoxm1()
2394 fp0 = floatx80_add(fp0, fp1, status); in floatx80_etoxm1()
2395 fp0 = floatx80_add(fp0, exp_tbl[j], status); in floatx80_etoxm1()
2397 fp0 = floatx80_add(fp0, float32_to_floatx80(exp_tbl2[j], in floatx80_etoxm1()
2399 fp0 = floatx80_add(fp0, exp_tbl[j], status); in floatx80_etoxm1()
2401 fp0 = floatx80_add(fp0, onebysc, status); in floatx80_etoxm1()
2404 fp0 = floatx80_add(fp0, float32_to_floatx80(exp_tbl2[j], in floatx80_etoxm1()
2408 fp0 = floatx80_add(fp0, fp1, status); in floatx80_etoxm1()
2416 a = floatx80_mul(fp0, sc, status); in floatx80_etoxm1()
2423 fp0 = float32_to_floatx80(make_float32(0xBF800000), in floatx80_etoxm1()
2429 a = floatx80_add(fp0, float32_to_floatx80( in floatx80_etoxm1()
2445 fp0 = a; in floatx80_etoxm1()
2446 fp0 = floatx80_mul(fp0, fp0, status); /* S = X*X */ in floatx80_etoxm1()
2449 fp1 = floatx80_mul(fp1, fp0, status); /* S * B12 */ in floatx80_etoxm1()
2455 fp2 = floatx80_mul(fp2, fp0, status); in floatx80_etoxm1()
2456 fp1 = floatx80_mul(fp1, fp0, status); in floatx80_etoxm1()
2463 fp2 = floatx80_mul(fp2, fp0, status); in floatx80_etoxm1()
2464 fp1 = floatx80_mul(fp1, fp0, status); in floatx80_etoxm1()
2471 fp2 = floatx80_mul(fp2, fp0, status); in floatx80_etoxm1()
2472 fp1 = floatx80_mul(fp1, fp0, status); in floatx80_etoxm1()
2479 fp2 = floatx80_mul(fp2, fp0, status); in floatx80_etoxm1()
2480 fp1 = floatx80_mul(fp1, fp0, status); in floatx80_etoxm1()
2486 fp2 = floatx80_mul(fp2, fp0, status); in floatx80_etoxm1()
2487 fp1 = floatx80_mul(fp1, fp0, status); in floatx80_etoxm1()
2489 fp2 = floatx80_mul(fp2, fp0, status); in floatx80_etoxm1()
2492 fp0 = floatx80_mul(fp0, float32_to_floatx80( in floatx80_etoxm1()
2496 fp0 = floatx80_add(fp0, fp1, status); /* S*B1+Q */ in floatx80_etoxm1()
2501 a = floatx80_add(fp0, a, status); in floatx80_etoxm1()
2508 fp0 = a; in floatx80_etoxm1()
2511 fp0 = floatx80_mul(fp0, float64_to_floatx80( in floatx80_etoxm1()
2514 fp0 = floatx80_add(fp0, sc, status); in floatx80_etoxm1()
2519 a = floatx80_mul(fp0, float64_to_floatx80( in floatx80_etoxm1()
2526 a = floatx80_add(fp0, sc, status); in floatx80_etoxm1()
2550 floatx80 fp0, fp1; in floatx80_tanh() local
2592 fp0 = float32_to_floatx80(make_float32(sign), status); in floatx80_tanh()
2599 a = floatx80_add(fp0, float32_to_floatx80(make_float32(sign), in floatx80_tanh()
2606 fp0 = packFloatx80(0, aExp + 1, aSig); /* Y = 2|X| */ in floatx80_tanh()
2607 fp0 = floatx80_etox(fp0, status); /* FP0 IS EXP(Y) */ in floatx80_tanh()
2608 fp0 = floatx80_add(fp0, float32_to_floatx80( in floatx80_tanh()
2613 sign ^ 0xC0000000), status), fp0, in floatx80_tanh()
2615 fp0 = float32_to_floatx80(make_float32(sign | 0x3F800000), in floatx80_tanh()
2621 a = floatx80_add(fp1, fp0, status); in floatx80_tanh()
2629 fp0 = packFloatx80(0, aExp + 1, aSig); /* Y = 2|X| */ in floatx80_tanh()
2630 fp0 = floatx80_etoxm1(fp0, status); /* FP0 IS Z = EXPM1(Y) */ in floatx80_tanh()
2631 fp1 = floatx80_add(fp0, float32_to_floatx80(make_float32(0x40000000), in floatx80_tanh()
2644 a = floatx80_div(fp0, fp1, status); in floatx80_tanh()
2666 floatx80 fp0, fp1, fp2; in floatx80_sinh() local
2700 fp0 = floatx80_abs(a); /* Y = |X| */ in floatx80_sinh()
2701 fp0 = floatx80_sub(fp0, float64_to_floatx80( in floatx80_sinh()
2704 fp0 = floatx80_sub(fp0, float64_to_floatx80( in floatx80_sinh()
2707 fp0 = floatx80_etox(fp0, status); in floatx80_sinh()
2713 a = floatx80_mul(fp0, fp2, status); in floatx80_sinh()
2720 fp0 = floatx80_abs(a); /* Y = |X| */ in floatx80_sinh()
2721 fp0 = floatx80_etoxm1(fp0, status); /* FP0 IS Z = EXPM1(Y) */ in floatx80_sinh()
2722 fp1 = floatx80_add(fp0, float32_to_floatx80(make_float32(0x3F800000), in floatx80_sinh()
2724 fp2 = fp0; in floatx80_sinh()
2725 fp0 = floatx80_div(fp0, fp1, status); /* Z/(1+Z) */ in floatx80_sinh()
2726 fp0 = floatx80_add(fp0, fp2, status); in floatx80_sinh()
2733 a = floatx80_mul(fp0, float32_to_floatx80(fact, status), status); in floatx80_sinh()
2754 floatx80 fp0, fp1; in floatx80_cosh() local
2784 fp0 = packFloatx80(0, aExp, aSig); in floatx80_cosh()
2785 fp0 = floatx80_sub(fp0, float64_to_floatx80( in floatx80_cosh()
2788 fp0 = floatx80_sub(fp0, float64_to_floatx80( in floatx80_cosh()
2791 fp0 = floatx80_etox(fp0, status); in floatx80_cosh()
2797 a = floatx80_mul(fp0, fp1, status); in floatx80_cosh()
2805 fp0 = packFloatx80(0, aExp, aSig); /* |X| */ in floatx80_cosh()
2806 fp0 = floatx80_etox(fp0, status); /* EXP(|X|) */ in floatx80_cosh()
2807 fp0 = floatx80_mul(fp0, float32_to_floatx80(make_float32(0x3F000000), in floatx80_cosh()
2810 fp1 = floatx80_div(fp1, fp0, status); /* 1/(2*EXP(|X|)) */ in floatx80_cosh()
2815 a = floatx80_add(fp0, fp1, status); in floatx80_cosh()