Lines Matching refs:s1
2126 fcvt.fx.s1 fXint = f8 // Convert arg to int (int repres. in FR)
2132 fcmp.lt.s1 p15, p14 = f8, f0
2137 fma.s1 fTwo = f1, f1, f1 // 2.0
2157 frcpa.s1 fRcpX, p0 = f1, f8 // initial approximation of 1/x
2181 fms.s1 fA5L = f1, f1, f8 // for 0.75 <= x < 1.3125 path
2196 fms.s1 fB4 = f8, f1, fTwo
2223 fnma.s1 fInvX = f8, fRcpX, f1 // start of 1st NR iteration
2239 fma.s1 fSix = fTwo, fTwo, fTwo // 6.0
2251 fma.s1 fRoot = f8, f1, f1 // x + 1
2282 fms.s1 fDx = f8, f1, fFltIntX // x - RTN(x)
2295 fma.s1 fRcpX = fInvX, fRcpX, fRcpX // end of 1st NR iteration
2303 fcmp.eq.s1 p9, p0 = f8,fFltIntX
2327 fcmp.lt.s1 p0, p15 = fAbsX, fSix // p15 is set when x < -6.0
2334 fma.s1 fThirteen = fSix, fTwo, f1 // 13.0
2343 fms.s1 FR_FracX = fSignifX, f1, FR_MHalf
2349 fma.s1 fDxSqr = fDx, fDx, f0 // deltaX^2
2364 fnma.s1 fInvX = f8, fRcpX, f1 // start of 2nd NR iteration
2372 fcmp.lt.s1 p10, p11 = fDx, f0
2407 fms.s1 fDxSqrL = fDx, fDx, fDxSqr // low part of deltaX^2
2421 fcmp.lt.s1 p14, p0 = fAbsX, fThirteen // set p14 if x > -13.0
2428 fma.s1 fXSqr = FR_FracX, FR_FracX, f0 // y^2
2435 fma.s1 fRcpX = fInvX, fRcpX, fRcpX // end of 2nd NR iteration
2442 fnma.s1 fMOne = f1, f1, f0 // -1.0
2464 fmpy.s1 FR_G = FR_G, FR_G2 // G = G_1 * G_2
2469 fadd.s1 FR_H = FR_H, FR_H2 // H = H_1 + H_2
2477 fma.s1 fDelX4 = fDxSqr, fDxSqr, f0 // deltaX^4
2483 fadd.s1 FR_h = FR_h, FR_h2 // h = h_1 + h_2
2489 fma.s1 fLnSin20 = fLnSin20, fDxSqr, fLnSin18
2495 fms.s1 fXSqrL = FR_FracX, FR_FracX, fXSqr // low part of y^2
2501 fma.s1 fX4 = fXSqr, fXSqr, f0 // y^4
2506 fma.s1 fLnSin24 = fLnSin24, fDxSqr, fLnSin22
2512 fma.s1 fLnSin28 = fLnSin28, fDxSqr, fLnSin26
2523 fma.s1 fLnSin32 = fLnSin32, fDxSqr, fLnSin30
2528 fma.s1 fB20 = f1, f1, FR_MHalf // 2.5
2534 fmpy.s1 FR_G = FR_G, FR_G3 // G = (G_1 * G_2) * G_3
2540 fadd.s1 FR_H = FR_H, FR_H3 // H = (H_1 + H_2) + H_3
2547 fma.s1 fLnSin36 = fLnSin36, fDxSqr, fLnSin34
2552 fma.s1 fLnSin20 = fLnSin20, fDxSqr, fLnSin16
2558 fma.s1 fX6 = fX4, fXSqr, f0 // y^6
2564 fma.s1 fDelX8 = fDelX4, fDelX4, f0 // deltaX^8
2571 fadd.s1 FR_h = FR_h, FR_h3 // h = (h_1 + h_2) + h_3
2576 fma.s1 fLnSin28 = fLnSin28, fDelX4, fLnSin24
2582 fma.s1 fRes1H = fA3, FR_FracX, f0 // (A3*y)hi
2587 fma.s1 fDelX6 = fDxSqr, fDelX4, f0 // DeltaX^6
2593 fma.s1 fA4 = fA5, FR_FracX, fA4 // A5*y + A4
2599 fms.s1 FR_r = FR_G, fSignifX, f1 // r = G * S_hi - 1
2610 fma.s1 fResH = fFloatN, FR_log2_hi, FR_H
2615 fma.s1 fRes6H = fA1, FR_FracX, f0 // (A1*y)hi
2621 fma.s1 fB14 = fLnSin6, fDxSqr, f0 // (LnSin6*deltaX^2)hi
2626 fma.s1 fA18 = fA19, FR_FracX, fA18
2632 fma.s1 fA23 = fA23, FR_FracX, fA22
2637 fma.s1 fA21 = fA21, FR_FracX, fA20
2643 fms.s1 fRes1L = fA3, FR_FracX, fRes1H // delta((A3*y)hi)
2648 fadd.s1 fPol = fRes1H, fA2 // (A3*y + A2)hi
2654 fma.s1 fRes2H = fA4, fXSqr, f0 // ((A5 + A4*y)*y^2)hi
2659 fma.s1 fA25 = fA25, FR_FracX, fA24
2666 fms.s1 fRes6L = fA1, FR_FracX, fRes6H // delta((A1*y)hi)
2672 fadd.s1 fRes3H = fRes6H, fA0 // (A1*y + A0)hi
2680 fadd.s1 fLnSinH = fB14, fLnSin4
2685 fms.s1 fB16 = fLnSin6, fDxSqr, fB14 // delta(LnSin6*deltaX^2)
2691 fma.s1 fA21 = fA21, fXSqr, fA18
2696 fma.s1 fLnSin36 = fLnSin36, fDelX4, fLnSin32
2702 fma.s1 fRes1L = fA3L, FR_FracX, fRes1L // (A3*y)lo
2708 fsub.s1 fPolL = fA2, fPol
2718 fms.s1 fRes2L = fA4, fXSqr, fRes2H
2724 fadd.s1 fRes4H = fRes2H, fPol
2731 fma.s1 fRes6L = fA1L, FR_FracX, fRes6L // (A1*y)lo
2737 fsub.s1 fRes3L = fA0, fRes3H
2743 fsub.s1 fLnSinL = fLnSin4, fLnSinH
2749 fma.s1 fB18 = fLnSinH, fDxSqr, f0
2755 fma.s1 fB16 = fLnSin6, fDxSqrL, fB16 // (LnSin6*deltaX^2)lo
2760 fma.s1 fA25 = fA25, fXSqr, fA23
2766 fadd.s1 fPolL = fPolL, fRes1H
2771 fadd.s1 fRes1L = fRes1L, fA2L // (A3*y)lo + A2lo
2777 fma.s1 fRes2L = fA4, fXSqrL, fRes2L // ((A5 + A4*y)*y^2)lo
2782 fsub.s1 fRes4L = fPol, fRes4H
2789 fma.s1 fRes7H = fRes4H, fXSqr, f0
2794 fma.s1 fA15 = fA15, FR_FracX, fA14
2800 fadd.s1 fRes3L = fRes3L, fRes6H
2805 fadd.s1 fRes6L = fRes6L, fA0L // (A1*y)lo + A0lo
2811 fadd.s1 fLnSinL = fLnSinL, fB14
2818 fms.s1 fB20 = fLnSinH, fDxSqr, fB18
2824 fadd.s1 fPolL = fPolL, fRes1L // (A3*y + A2)lo
2831 fadd.s1 fLnSin6 = fB18, fLnSin2
2837 fadd.s1 fRes4L = fRes4L, fRes2H
2842 fma.s1 fA17 = fA17, FR_FracX, fA16
2849 fms.s1 fRes7L = fRes4H, fXSqr, fRes7H
2854 fadd.s1 fPol = fRes7H, fRes3H
2860 fadd.s1 fRes3L = fRes3L, fRes6L // (A1*y + A0)lo
2865 fma.s1 fA25 = fA25, fX4, fA21
2872 fadd.s1 fLnSinL = fLnSinL, fB16
2877 fma.s1 fB20 = fLnSinH, fDxSqrL, fB20
2883 fsub.s1 fLnSin4 = fLnSin2, fLnSin6
2889 fma.s1 fLnSinH = fLnSin6, fDxSqr, f0
2896 fadd.s1 fRes2L = fRes2L, fPolL
2901 fma.s1 fA17 = fA17, fXSqr, fA15
2908 fma.s1 fRes7L = fRes4H, fXSqrL, fRes7L
2913 fsub.s1 fPolL = fRes3H, fPol
2919 fma.s1 fA13 = fA13, FR_FracX, fA12
2924 fma.s1 fA11 = fA11, FR_FracX, fA10
2931 fma.s1 fB20 = fLnSinL, fDxSqr, fB20
2936 fmpy.s1 FR_G = FR_G, FR_G2 // G = G_1 * G_2
2942 fadd.s1 fLnSin4 = fLnSin4, fB18
2947 fms.s1 fLnSinL = fLnSin6, fDxSqr, fLnSinH
2954 fadd.s1 fRes4L = fRes4L, fRes2L
2959 fadd.s1 fhDelX = fhDelX, FR_h2 // h = h_1 + h_2
2965 fadd.s1 fRes7L = fRes7L, fRes3L
2970 fadd.s1 fPolL = fPolL, fRes7H
2981 fadd.s1 FR_H = FR_H, FR_H2 // H = H_1 + H_2
2987 fmpy.s1 FR_G = FR_G, FR_G3 // G = (G_1 * G_2) * G_3
2993 fadd.s1 fLnSin2L = fLnSin2L, fB20
2999 fma.s1 fA25 = fA25, fX4, fA17
3004 fma.s1 fA13 = fA13, fXSqr, fA11
3010 fma.s1 fA9 = fA9, FR_FracX, fA8
3015 fma.s1 fA7 = fA7, FR_FracX, fA6
3021 fma.s1 fLnSin36 = fLnSin36, fDelX8, fLnSin28
3026 fma.s1 fLnSin14 = fLnSin14, fDxSqr, fLnSin12
3032 fma.s1 fLnSin10 = fLnSin10, fDxSqr, fLnSin8
3037 fadd.s1 FR_H = FR_H, FR_H3 // H = (H_1 + H_2) + H_3
3043 fms.s1 fRDx = FR_G, fNormDx, f1 // r = G * S_hi - 1
3049 fma.s1 FR_poly_lo = FR_r, FR_Q4, FR_Q3
3055 fmpy.s1 FR_rsq = FR_r, FR_r // rsq = r * r
3061 fma.s1 fRes7L = fRes4L, fXSqr, fRes7L
3067 fma.s1 fA25 = fA25, fX4, fA13
3072 fma.s1 fA9 = fA9, fXSqr, fA7
3079 fma.s1 FR_h = fFloatN, FR_log2_lo, FR_h
3084 fadd.s1 fhDelX = fhDelX, FR_h3 // h = (h_1 + h_2) + h_3
3090 fma.s1 fLnSin36 = fLnSin36, fDelX6, fLnSin20
3095 fma.s1 fLnSin14 = fLnSin14, fDelX4, fLnSin10
3102 fma.s1 fPolyLoDx = fRDx, FR_Q4, FR_Q3
3107 fmpy.s1 fRDxSq = fRDx, fRDx // rsq = r * r
3114 fma.s1 fResLnDxH = fFloatNDx, FR_log2_hi, FR_H
3119 fma.s1 FR_rcub = FR_rsq, FR_r, f0 // rcub = r^3
3125 fma.s1 fA9 = fA25, fX4, fA9
3130 fadd.s1 fPolL = fPolL, fRes7L
3136 fadd.s1 fLnSin4 = fLnSin4, fLnSin2L
3142 fma.s1 fhDelX = fFloatNDx, FR_log2_lo, fhDelX
3148 fma.s1 fLnSin36 = fLnSin36, fDelX8, fLnSin14
3154 fma.s1 fLnSinL = fLnSin6, fDxSqrL, fLnSinL
3161 fma.s1 fPolyLoDx = fPolyLoDx, fRDx, FR_Q2
3166 fma.s1 fRDxCub = fRDxSq, fRDx, f0 // rcub = r^3
3178 fadd.s1 fRes1H = fPol, fResH
3185 fma.s1 FR_poly_lo = FR_poly_lo, FR_r, FR_Q2
3190 fma.s1 fPolL = fA9, fX6, fPolL // P25lo
3203 fnma.s1 fRes2H = fResLnDxH, f1, fLnSinH
3211 fma.s1 fLnSinL = fLnSin4, fDxSqr, fLnSinL
3216 fma.s1 fLnSin36 = fLnSin36, fDelX6, f0
3223 fma.s1 fPolyHiDx = FR_Q1, fRDxSq, fRDx
3229 fma.s1 fPolyLoDx = fPolyLoDx, fRDxCub, fhDelX
3235 fsub.s1 fRes1L = fRes5H, fRes1H
3241 fnma.s1 fRes1H = fRes1H, f1, f0
3249 fma.s1 FR_poly_hi = FR_Q1, FR_rsq, FR_r
3255 fma.s1 FR_poly_lo = FR_poly_lo, FR_rcub, FR_h
3261 fms.s1 fRes2L = fResLnDxH, fMOne, fRes2H
3267 fma.s1 fLnSinL = fLnSin36, fDxSqr, fLnSinL
3273 fadd.s1 fResLnDxL = fPolyHiDx, fPolyLoDx
3279 fadd.s1 fRes1L = fRes1L, fRes5L
3285 fadd.s1 fYH = fRes2H, fRes1H
3292 fadd.s1 fResL = FR_poly_hi, FR_poly_lo
3311 fsub.s1 fLnSinL = fLnSinL, fResLnDxL
3316 fadd.s1 fRes2L = fRes2L, fLnSinH
3323 fadd.s1 fPolL = fResL, fPolL
3329 fsub.s1 fYL = fRes4H, fYH
3336 fadd.s1 fRes2L = fRes2L, fLnSinL
3342 fadd.s1 fRes1L = fRes1L, fPolL
3348 fadd.s1 fYL = fYL, fRes4L
3353 fsub.s1 fRes2L = fRes2L, fRes1L
3360 fadd.s1 fYL = fYL, fRes2L
3384 fma.s1 fA0L = fDxSqr, fDxSqr, f0 // deltaX^4
3391 fnma.s1 fResL = fDx, f1, f0 // -(x+1)
3543 fma.s1 fResH = fA20, fResL, fA19 //polynomial for log(|x|)
3550 fma.s1 fA2 = fA2, fDx, fA1 // polynomial for lgammal(|x|)
3559 fma.s1 fA18 = fA18, fResL, fA17 //polynomial for log(|x|)
3565 fma.s1 fA16 = fA16, fResL, fA15 //polynomial for log(|x|)
3570 fma.s1 fA4 = fA4, fDx, fA3 // polynomial for lgammal(|x|)
3576 fma.s1 fA14 = fA14, fResL, fA13 //polynomial for log(|x|)
3581 fma.s1 fA6 = fA6, fDx, fA5L // polynomial for lgammal(|x|)
3587 fma.s1 fPol = fA8, fDx, fA7 // polynomial for lgammal(|x|)
3594 fma.s1 fRes3L = fLnSin4, fDxSqr, fLnSin2
3605 fma.s1 fResH = fResH, fDxSqr, fA18 // High part of log(|x|)
3611 fma.s1 fA4 = fA4, fDxSqr, fA2 // Low part of lgammal(|x|)
3617 fma.s1 fRes3H = fLnSin8, fDxSqr, fLnSin6
3623 fmpy.s1 FR_G = FR_G, FR_G2 // G = G_1 * G_2
3628 fadd.s1 FR_H = FR_H, FR_H2 // H = H_1 + H_2
3634 fadd.s1 FR_h = FR_h, FR_h2 // h = h_1 + h_2
3639 fma.s1 fA16 = fA16, fDxSqr, fA14 // Low part of log(|x|)
3645 fma.s1 fPol = fPol, fDxSqr, fA6 // High part of lgammal(|x|)
3651 fma.s1 fResH = fResH, fA0L, fA16 // log(|x|)/deltaX^2 - deltaX
3657 fmpy.s1 FR_G = FR_G, FR_G3 // G = (G_1 * G_2) * G_3
3662 fadd.s1 FR_H = FR_H, FR_H3 // H = (H_1 + H_2) + H_3
3668 fadd.s1 FR_h = FR_h, FR_h3 // h = (h_1 + h_2) + h_3
3674 fma.s1 fResH = fResH, fDxSqr, fResL // log(|x|)
3679 fma.s1 fPol = fPol, fA0L, fA4 // lgammal(|x|)/|x|
3685 fms.s1 FR_r = FR_G, fAbsX, f1 // r = G * S_hi - 1
3691 fma.s1 fRes4H = fFloatN, FR_log2_hi, FR_H
3698 fma.s1 FR_h = fFloatN, FR_log2_lo, FR_h
3704 fma.s1 fResH = fPol, fDx, fResH // lgammal(|x|) + log(|x|)
3710 fma.s1 fRes3H = fRes3H, fA0L, fRes3L
3717 fma.s1 FR_poly_lo = FR_r, FR_Q4, FR_Q3
3722 fmpy.s1 FR_rsq = FR_r, FR_r // rsq = r * r
3729 fms.s1 fResH = fRes3H, fDxSqr, fResH
3737 fma.s1 FR_poly_lo = FR_poly_lo, FR_r, FR_Q2
3742 fma.s1 FR_rcub = FR_rsq, FR_r, f0 // rcub = r^3
3750 fma.s1 FR_poly_hi = FR_Q1, FR_rsq, FR_r
3758 fma.s1 FR_poly_lo = FR_poly_lo, FR_rcub, FR_h
3766 fadd.s1 fRes4L = FR_poly_hi, FR_poly_lo
3772 fsub.s1 fResH = fResH, fRes4L
3796 fnma.s1 fInvX = f8, fRcpX, f1 // start of 3rd NR iteration
3802 fadd.s1 FR_h = FR_h, FR_h2 // h = h_1 + h_2
3808 fmpy.s1 FR_G = FR_G, FR_G3 // G = (G_1 * G_2) * G_3
3814 fadd.s1 FR_H = FR_H, FR_H3 // H = (H_1 + H_2) + H_3
3821 fma.s1 fDx4 = fDxSqr, fDxSqr, f0 // deltaX^4
3826 fma.s1 fLnSin34 = fLnSin34, fDxSqr, fLnSin32
3834 fma.s1 fLnSin26 = fLnSin26, fDxSqr, fLnSin24
3839 fma.s1 fLnSin30 = fLnSin30, fDxSqr, fLnSin28
3845 fma.s1 fLnSin18 = fLnSin18, fDxSqr, fLnSin16
3851 fma.s1 fLnSin22 = fLnSin22, fDxSqr, fLnSin20
3857 fma.s1 fInvX = fInvX, fRcpX, fRcpX // end of 3rd NR iteration
3863 fadd.s1 FR_h = FR_h, FR_h3 // h = (h_1 + h_2) + h_3
3872 fms.s1 FR_r = FR_G, fSignifX, f1 // r = G * S_hi - 1
3878 fma.s1 fResH = fFloatN, FR_log2_hi, FR_H
3885 fma.s1 fLnSin36 = fLnSin36, fDx4, fLnSin34
3891 fma.s1 fRes5H = fLnSin4, fDxSqr, f0
3899 fma.s1 fLnSin30 = fLnSin30, fDx4, fLnSin26
3904 fma.s1 fDx8 = fDx4, fDx4, f0 // deltaX^8
3910 fma.s1 fDx6 = fDx4, fDxSqr, f0 // deltaX^6
3915 fma.s1 fLnSin22 = fLnSin22, fDx4, fLnSin18
3922 fma.s1 fRcpX = fInvX, fInvX, f0 // (1/x)^2
3927 fma.s1 FR_h = fFloatN, FR_log2_lo, FR_h// h = N * log2_lo + h
3933 fnma.s1 fInvXL = f8, fInvX, f1 // relative error of 1/x
3938 fma.s1 fLnSin8 = fLnSin8, fDxSqr, fLnSin6
3945 fma.s1 FR_poly_lo = FR_r, FR_Q4, FR_Q3
3950 fmpy.s1 FR_rsq = FR_r, FR_r // rsq = r * r
3956 fma.s1 fA1L = fB2, fInvX, f0 // (B2*(1/x))hi
3962 fms.s1 fRes4H = fResH, f1, f1 // ln(|x|)hi - 1
3968 fadd.s1 fRes2H = fRes5H, fLnSin2//(lnSin4*DeltaX^2 + lnSin2)hi
3974 fms.s1 fRes5L = fLnSin4, fDxSqr, fRes5H
3983 fma.s1 fInvX4 = fRcpX, fRcpX, f0 // (1/x)^4
3988 fma.s1 fB6 = fB6, fRcpX, fB4
3995 fma.s1 fB18 = fB18, fRcpX, fB16
4001 fma.s1 fInvXL = fInvXL, fInvX, f0 // low part of 1/x
4008 fma.s1 FR_poly_lo = FR_poly_lo, FR_r, FR_Q2
4013 fma.s1 FR_rcub = FR_rsq, FR_r, f0 // rcub = r^3
4019 fma.s1 fRes3H = fRes4H, f8, f0 // (-|x|*(ln(|x|)-1))hi
4025 fma.s1 FR_poly_hi = FR_Q1, FR_rsq, FR_r
4031 fms.s1 fA2L = fB2, fInvX, fA1L // delta(B2*(1/x))
4036 fnma.s1 fBrnH = fRes1H, f1, fA1L // (-C - S(1/x))hi
4042 fma.s1 fInvX8 = fInvX4, fInvX4, f0 // (1/x)^8
4047 fma.s1 fB10 = fB10, fRcpX, fB8
4054 fma.s1 fB20 = fB20, fInvX4, fB18
4059 fma.s1 fB14 = fB14, fRcpX, fB12
4065 fma.s1 fLnSin36 = fLnSin36, fDx8, fLnSin30
4070 fma.s1 fLnSin12 = fLnSin12, fDxSqr, fLnSin10
4076 fsub.s1 fRes2L = fLnSin2, fRes2H
4081 fma.s1 fPol = fRes2H, fDxSqr, f0 // high part of LnSin
4087 fnma.s1 fResH = fResH, FR_MHalf, fResH // -0.5*ln(|x|)hi
4092 fmpy.s1 fGDx = fGDx, FR_G2 // G = G_1 * G_2
4099 fma.s1 FR_poly_lo = FR_poly_lo, FR_rcub, FR_h
4105 fma.s1 fA2L = fB2L, fInvX, fA2L
4111 fma.s1 fB20 = fB20, fInvX4, fB14
4116 fma.s1 fB10 = fB10, fInvX4, fB6
4127 fma.s1 fLnSin14 = fLnSin14, fDx4, fLnSin12
4133 fma.s1 fLnSin36 = fLnSin36, fDx8, fLnSin22
4138 fms.s1 fRes3L = fRes4H, f8, fRes3H // delta(-|x|*(ln(|x|)-1))
4144 fmpy.s1 fGDx = fGDx, fG3Dx // G = (G_1 * G_2) * G_3
4150 fadd.s1 fRes4H = fRes3H, fResH
4156 fma.s1 fA2L = fInvXL, fB2, fA2L //(B2*(1/x))lo
4162 fadd.s1 fResL = FR_poly_hi, FR_poly_lo
4168 fma.s1 fB20 = fB20, fInvX8, fB10
4173 fma.s1 fInvX3 = fInvX, fRcpX, f0 // (1/x)^3
4179 fadd.s1 fHDx = fHDx, FR_H2 // H = H_1 + H_2
4184 fadd.s1 fRes5L = fRes5L, fLnSin2L
4190 fadd.s1 fRes2L = fRes2L, fRes5H
4195 fadd.s1 fhDx = fhDx, fh2Dx // h = h_1 + h_2
4201 fms.s1 fBrnL = fRes1H, fMOne, fBrnH
4206 fms.s1 FR_r = fGDx, fNormDx, f1 // r = G * S_hi - 1
4212 fma.s1 fRes3L = fResL, f8 , fRes3L // (-|x|*(ln(|x|)-1))lo
4217 fsub.s1 fRes4L = fRes3H, fRes4H
4224 fma.s1 fB20 = fB20, fInvX3, fA2L
4229 fnma.s1 fResL = fResL, FR_MHalf, fResL // -0.5*ln(|x|)lo
4235 fadd.s1 fHDx = fHDx, fH3Dx // H = (H_1 + H_2) + H_3
4240 fms.s1 fPolL = fRes2H, fDxSqr, fPol
4246 fadd.s1 fhDx = fhDx, fh3Dx // h = (h_1 + h_2) + h_3
4252 fadd.s1 fB14 = fRes4H, fBrnH
4259 fma.s1 FR_poly_lo = FR_r, FR_Q4, FR_Q3
4264 fmpy.s1 FR_rsq = FR_r, FR_r // rsq = r * r
4270 fadd.s1 fRes4L = fRes4L, fResH
4275 fadd.s1 fBrnL = fBrnL, fA1L
4282 fadd.s1 fRes3L = fRes3L, fResL
4287 fnma.s1 fB20 = fRes1L, f1, fB20 // -Clo - S(1/x)lo
4293 fadd.s1 fRes2L = fRes2L, fRes5L // (lnSin4*DeltaX^2 + lnSin2)lo
4298 fma.s1 fPolL = fDxSqrL, fRes2H, fPolL
4304 fma.s1 fLnSin14 = fLnSin14, fDx4, fLnSin8
4309 fma.s1 fLnSin36 = fLnSin36, fDx8, f0
4316 fma.s1 FR_poly_lo = FR_poly_lo, FR_r, FR_Q2
4321 fma.s1 FR_rcub = FR_rsq, FR_r, f0 // rcub = r^3
4328 fma.s1 FR_poly_hi = FR_Q1, FR_rsq, FR_r
4333 fsub.s1 fB12 = fRes4H, fB14
4340 fadd.s1 fRes4L = fRes4L, fRes3L
4345 fadd.s1 fBrnL = fBrnL, fB20 // (-C - S(1/x))lo
4352 fma.s1 fLnDeltaH = fFloatNDx, FR_log2_hi, fHDx
4358 fma.s1 fhDx = fFloatNDx, FR_log2_lo, fhDx
4364 fma.s1 fPolL = fRes2L, fDxSqr, fPolL
4369 fma.s1 fLnSin14 = fLnSin36, fDxSqr, fLnSin14
4376 fadd.s1 fBrnL = fBrnL, fRes4L
4381 fadd.s1 fB12 = fB12, fBrnH
4388 fma.s1 FR_poly_lo = FR_poly_lo, FR_rcub, fhDx
4393 fnma.s1 fRes1H = fLnDeltaH, f1, fPol//(-ln(|DeltaX|) + LnSin)hi
4399 fma.s1 fPolL = fDxSqrL, fRes2L, fPolL
4404 fma.s1 fLnSin36 = fLnSin14, fDx6, f0
4411 fadd.s1 fB12 = fB12, fBrnL
4418 fadd.s1 fLnDeltaL= FR_poly_hi, FR_poly_lo
4423 fms.s1 fRes1L = fLnDeltaH, fMOne, fRes1H
4429 fadd.s1 fPolL = fPolL, fLnSin36
4435 fadd.s1 f8 = fRes1H, fB14
4443 famax.s1 fMaxNegStir = fRes1H, fB14
4450 famin.s1 fMinNegStir = fRes1H, fB14
4456 fadd.s1 fRes1L = fRes1L, fPol
4462 fnma.s1 fPolL = fLnDeltaL, f1, fPolL
4468 fsub.s1 f9 = fMaxNegStir, f8 // delta1
4474 fadd.s1 fRes1L = fRes1L, fPolL // (-ln(|DeltaX|) + LnSin)lo
4480 fadd.s1 f9 = f9, fMinNegStir
4486 fadd.s1 fRes1L = fRes1L, fB12
4492 fadd.s1 f9 = f9, fRes1L
4634 fcmp.eq.s1 p6, p0 = f8, fRoot
4639 fms.s1 FR_FracX = f8, f1, fRoot
4666 fma.s1 fResH = fA2, FR_FracX, f0 // (A2*x)hi
4671 fma.s1 fA4L = FR_FracX, FR_FracX, f0 // x^2
4677 fma.s1 fA17 = fA17, FR_FracX, fA16
4682 fma.s1 fA13 = fA13, FR_FracX, fA12
4688 fma.s1 fA19 = fA19, FR_FracX, fA18
4693 fma.s1 fA15 = fA15, FR_FracX, fA14
4699 fma.s1 fPol = fA7, FR_FracX, fA6
4705 fma.s1 fA9 = fA9, FR_FracX, fA8
4711 fms.s1 fResL = fA2, FR_FracX, fResH // delta(A2*x)
4716 fadd.s1 fRes1H = fResH, fA1 // (A2*x + A1)hi
4722 fma.s1 fA11 = fA11, FR_FracX, fA10
4727 fma.s1 fA5L = fA4L, fA4L, f0 // x^4
4733 fma.s1 fA19 = fA19, fA4L, fA17
4738 fma.s1 fA15 = fA15, fA4L, fA13
4744 fma.s1 fPol = fPol, FR_FracX, fA5
4749 fma.s1 fA3L = fA4L, FR_FracX, f0 // x^3
4756 fma.s1 fResL = fA2L, FR_FracX, fResL
4761 fsub.s1 fRes1L = fA1, fRes1H
4767 fma.s1 fA11 = fA11, fA4L, fA9
4772 fma.s1 fA19 = fA19, fA5L, fA15
4778 fma.s1 fPol = fPol, FR_FracX, fA4
4784 fadd.s1 fResL = fResL, fA1L // (A2*x)lo + A1
4789 fadd.s1 fRes1L = fRes1L, fResH
4795 fma.s1 fRes2H = fRes1H, FR_FracX, f0 // ((A2*x + A1)*x)hi
4801 fma.s1 fA19 = fA19, fA5L, fA11
4807 fma.s1 fPol = fPol, FR_FracX, fA3
4813 fadd.s1 fRes1L = fRes1L, fResL // (A2*x + A1)lo
4820 fms.s1 fRes2L = fRes1H, FR_FracX, fRes2H
4825 fadd.s1 fRes3H = fRes2H, fA0 // ((A2*x + A1)*x + A0)hi
4831 fma.s1 fA19 = fA19, fA5L, f0
4838 fma.s1 fRes2L = fRes1L, FR_FracX, fRes2L // ((A2*x + A1)*x)lo
4843 fsub.s1 fRes3L = fRes2H, fRes3H
4849 fma.s1 fPol = fA19, FR_FracX, fPol
4855 fadd.s1 fRes3L = fRes3L, fA0
4860 fadd.s1 fRes2L = fRes2L, fA0L // ((A2*x + A1)*x)lo + A0L
4866 fadd.s1 fRes3L = fRes3L, fRes2L // (((A2*x + A1)*x) + A0)lo
4872 fma.s1 fRes3L = fPol, fA3L, fRes3L
4890 fma.s1 fA4L = f8, f8, f0 // x^2
4953 fma.s1 fA5L = fA4L, fA4L, f0 // x^4
4958 fms.s1 fB8 = f8, f8, fA4L // x^2 - <x^2>
5013 fma.s1 fB20 = fA5L, fA5L, f0 // x^8
5019 fms.s1 fB16 = fA4L, fA4L, fA5L
5040 fma.s1 fRes4H = fA3, fAbsX, f0
5052 fma.s1 fB18 = fB20, fB20, f0 // x^16
5063 fma.s1 fA21 = fA21, fAbsX, fA20 // v16
5076 fma.s1 fA19 = fA19, fAbsX, fA18 // v13
5082 fms.s1 fRes4L = fA3, fAbsX, fRes4H
5087 fadd.s1 fRes2H = fRes4H, fA2
5094 fma.s1 fA17 = fA17, fAbsX, fA16 // v12
5100 fmpy.s1 FR_G = FR_G, FR_G2 // G = G_1 * G_2
5105 fadd.s1 FR_H = FR_H, FR_H2 // H = H_1 + H_2
5111 fma.s1 fA15 = fA15, fAbsX, fA14 // v8
5116 fadd.s1 FR_h = FR_h, FR_h2 // h = h_1 + h_2
5123 fma.s1 fA13 = fA13, fAbsX, fA12 // v7
5129 fma.s1 fRes4L = fA3L, fAbsX, fRes4L
5134 fsub.s1 fRes2L = fA2, fRes2H
5140 fma.s1 fResH = fRes2H, fAbsX, f0
5145 fma.s1 fA22 = fA22, fA4L, fA21 // v15
5151 fma.s1 fA9 = fA9, fAbsX, fA8 // v4
5156 fma.s1 fA11 = fA11, fAbsX, fA10 // v5
5162 fma.s1 fA19 = fA19, fA4L, fA17 // v11
5167 fma.s1 fPolL = fA7, fAbsX, fA6
5173 fmpy.s1 FR_G = FR_G, FR_G3 // G = (G_1 * G_2) * G_3
5178 fadd.s1 FR_H = FR_H, FR_H3 // H = (H_1 + H_2) + H_3
5184 fadd.s1 fRes2L = fRes2L, fRes4H
5189 fadd.s1 fA2L = fA2L, fRes4L
5195 fms.s1 fResL = fRes2H, fAbsX, fResH
5200 fadd.s1 fRes2H = fResH, fA1
5206 fma.s1 fA11 = fA11, fA4L, fA9 // v3
5211 fma.s1 fA15 = fA15, fA4L, fA13 // v6
5219 fadd.s1 FR_h = FR_h, FR_h3 // h = (h_1 + h_2) + h_3
5225 fma.s1 fA5 = fA5, fAbsX, fA4
5231 fms.s1 FR_r = FR_G, fSignifX, f1 // r = G * S_hi - 1
5237 fms.s1 FR_log2_hi = fFloatN, FR_log2_hi, FR_H
5243 fadd.s1 fA3L = fRes2L, fA2L
5248 fma.s1 fA22 = fA22, fA5L, fA19
5254 fsub.s1 fRes2L = fA1, fRes2H
5259 fma.s1 fRes3H = fRes2H, f8, f0
5265 fma.s1 fA15 = fA15, fA5L, fA11 // v2
5270 fma.s1 fLnSin18 = fLnSin18, fA4L, fLnSin16
5277 fms.s1 FR_h = fFloatN, FR_log2_lo, FR_h
5282 fma.s1 fPolL = fPolL, fA4L, fA5
5289 fma.s1 FR_poly_lo = FR_r, FR_Q4, FR_Q3
5294 fmpy.s1 FR_rsq = FR_r, FR_r // rsq = r * r
5300 fma.s1 fResL = fA3L, fAbsX, fResL
5305 fma.s1 fLnSin30 = fLnSin30, fA4L, fLnSin28
5311 fadd.s1 fRes2L = fRes2L, fResH
5316 fms.s1 fRes3L = fRes2H, f8, fRes3H
5322 fadd.s1 fRes1H = fRes3H, FR_log2_hi
5327 fma.s1 fPol = fB20, fA22, fA15
5333 fma.s1 fLnSin34 = fLnSin34, fA4L, fLnSin32
5338 fma.s1 fLnSin14 = fLnSin14, fA4L, fLnSin12
5346 fma.s1 FR_poly_lo = FR_poly_lo, FR_r, FR_Q2
5351 fnma.s1 FR_rcub = FR_rsq, FR_r, f0 // rcub = r^3
5358 fma.s1 FR_poly_hi = FR_Q1, FR_rsq, FR_r
5363 fadd.s1 fA1L = fA1L, fResL
5370 fma.s1 fLnSin22 = fLnSin22, fA4L, fLnSin20
5375 fma.s1 fLnSin26 = fLnSin26, fA4L, fLnSin24
5382 fsub.s1 fRes1L = FR_log2_hi, fRes1H
5387 fma.s1 fPol = fPol, fA5L, fPolL
5393 fma.s1 fLnSin34 = fLnSin36, fA5L, fLnSin34
5398 fma.s1 fLnSin18 = fLnSin18, fA5L, fLnSin14
5404 fma.s1 fLnSin6 = fLnSin6, fA4L, fLnSin4
5409 fma.s1 fLnSin10 = fLnSin10, fA4L, fLnSin8
5416 fma.s1 FR_poly_hi = FR_Q1, FR_rsq, FR_r
5421 fadd.s1 fRes2L = fRes2L, fA1L
5428 fma.s1 FR_poly_lo = FR_poly_lo, FR_rcub, FR_h
5433 fma.s1 fB2 = fLnSin2, fA4L, f0
5439 fadd.s1 fRes1L = fRes1L, fRes3H
5444 fma.s1 fPol = fPol, fB10, f0
5450 fma.s1 fLnSin26 = fLnSin26, fA5L, fLnSin22
5455 fma.s1 fLnSin34 = fLnSin34, fA5L, fLnSin30
5461 fma.s1 fLnSin10 = fLnSin10, fA5L, fLnSin6
5466 fma.s1 fLnSin2L = fLnSin2L, fA4L, f0
5473 fma.s1 fRes3L = fRes2L, f8, fRes3L
5480 fsub.s1 FR_log2_lo = FR_poly_lo, FR_poly_hi
5485 fms.s1 fB4 = fLnSin2, fA4L, fB2
5491 fadd.s1 fRes2H = fRes1H, fPol
5497 fma.s1 fLnSin34 = fLnSin34, fB20, fLnSin26
5503 fma.s1 fLnSin18 = fLnSin18, fB20, fLnSin10
5508 fma.s1 fLnSin2L = fB8, fLnSin2, fLnSin2L
5515 fadd.s1 FR_log2_lo = FR_log2_lo, fRes3L
5521 fsub.s1 fRes2L = fRes1H, fRes2H
5527 fma.s1 fB6 = fLnSin34, fB18, fLnSin18
5532 fadd.s1 fB4 = fLnSin2L, fB4
5539 fadd.s1 fRes1L = fRes1L, FR_log2_lo
5545 fadd.s1 fRes2L = fRes2L, fPol
5551 fma.s1 fB12 = fB6, fA5L, f0
5557 fadd.s1 fRes2L = fRes2L, fRes1L
5564 fms.s1 fB14 = fB6, fA5L, fB12
5569 fadd.s1 fLnSin30 = fB2, fB12
5589 fma.s1 fB14 = fB16, fB6, fB14
5594 fsub.s1 fLnSin16 = fB2, fLnSin30
5600 fadd.s1 fResH = fLnSin30, fRes2H
5606 fadd.s1 fLnSin16 = fLnSin16, fB12
5611 fadd.s1 fB4 = fB14, fB4
5617 fadd.s1 fLnSin16 = fB4, fLnSin16
5622 fsub.s1 fResL = fRes2H, fResH
5628 fadd.s1 fResL = fResL, fLnSin30
5633 fadd.s1 fLnSin16 = fLnSin16, fRes2L
5639 fadd.s1 fResL = fResL, fLnSin16
5666 fnma.s1 fInvX = f8, fRcpX, f1 // start of 1st NR iteration
5721 fma.s1 fRcpX = fInvX, fRcpX, fRcpX // end of 1st NR iteration
5772 fnma.s1 fInvX = f8, fRcpX, f1 // start of 2nd NR iteration
5800 fcmp.gt.s1 p7,p0 = f8, fA1 // check if x > overflow threshold
5806 fma.s1 fRcpX = fInvX, fRcpX, fRcpX // end of 2nd NR iteration
5841 fmpy.s1 FR_G = FR_G, FR_G2 // G = G_1 * G_2
5846 fadd.s1 FR_H = FR_H, FR_H2 // H = H_1 + H_2
5853 fadd.s1 FR_h = FR_h, FR_h2 // h = h_1 + h_2
5858 fnma.s1 fInvX = f8, fRcpX, f1 // start of 3rd NR iteration
5885 fadd.s1 fB2 = fB2, fA1
5891 fma.s1 fInvX = fInvX, fRcpX, fRcpX // end of 3rd NR iteration
5904 fmpy.s1 FR_G = FR_G, FR_G3 // G = (G_1 * G_2) * G_3
5910 fadd.s1 FR_H = FR_H, FR_H3 // H = (H_1 + H_2) + H_3
5916 fadd.s1 FR_h = FR_h, FR_h3 // h = (h_1 + h_2) + h_3
5922 fma.s1 fRcpX = fInvX, fInvX, f0 // 1/x^2
5927 fma.s1 fA0L = fB2, fInvX, fA0L
5933 fms.s1 FR_r = fSignifX, FR_G, f1 // r = G * S_hi - 1
5939 fma.s1 fRes2H = fFloatN, FR_log2_hi, FR_H
5947 fma.s1 FR_h = fFloatN, FR_log2_lo, FR_h
5953 fma.s1 fRes1H = fFloatN, FR_log2_hi, FR_H
5959 fma.s1 fPol = fB18, fRcpX, fB16 // v9
5964 fma.s1 fA2L = fRcpX, fRcpX, f0 // v10
5970 fma.s1 fA3 = fB6, fRcpX, fB4 // v3
5975 fma.s1 fA4 = fB10, fRcpX, fB8 // v4
5981 fms.s1 fRes2H =fRes2H, f1, f1 // log_Hi(x) -1
5987 fma.s1 FR_poly_lo = FR_r, FR_Q4, FR_Q3
5993 fma.s1 fRes1H = fRes1H, FR_MHalf, f0 // -0.5*log_Hi(x)
5998 fmpy.s1 FR_rsq = FR_r, FR_r // rsq = r * r
6004 fma.s1 fA7 = fB14, fRcpX, fB12 // v7
6009 fma.s1 fA8 = fA2L, fB20, fPol // v8
6015 fma.s1 fA2 = fA4, fA2L, fA3 // v2
6020 fma.s1 fA4L = fA2L, fA2L, f0 // v5
6026 fma.s1 fResH = fRes2H, f8, f0 // (x*(ln(x)-1))hi
6032 fma.s1 FR_poly_lo = FR_poly_lo, FR_r, FR_Q2
6038 fma.s1 FR_rcub = FR_rsq, FR_r, f0 // rcub = r^3
6044 fma.s1 FR_poly_hi = FR_Q1, FR_rsq, FR_r
6050 fma.s1 fA11 = fRcpX, fInvX, f0 // 1/x^3
6055 fma.s1 fA6 = fA8, fA2L, fA7 // v6
6061 fms.s1 fResL = fRes2H, f8, fResH // d(x*(ln(x)-1))
6066 fadd.s1 fRes3H = fResH, fRes1H // (x*(ln(x)-1) -0.5ln(x))hi
6073 fma.s1 FR_poly_lo = FR_poly_lo, FR_rcub, FR_h
6079 fma.s1 fPol = fA4L, fA6, fA2 // v1
6091 fadd.s1 fRes4H = fRes3H, fA0 // (x*(ln(x)-1) -0.5ln(x))hi + Chi
6096 fsub.s1 fRes3L = fResH, fRes3H
6103 fadd.s1 fRes2L = FR_poly_hi, FR_poly_lo
6110 fma.s1 fA0L = fPol, fA11, fA0L // S(1/x) + Clo
6116 fadd.s1 fRes3L = fRes3L, fRes1H
6121 fsub.s1 fRes4L = fRes3H, fRes4H
6127 fma.s1 fResL = fRes2L, f8 , fResL // lo part of x*(ln(x)-1)
6134 fma.s1 fA0L = fRes2L, FR_MHalf, fA0L
6140 fadd.s1 fRes4L = fRes4L, fA0
6147 fadd.s1 fA0L = fA0L, fResL
6153 fadd.s1 fRes4L = fRes4L, fRes3L
6159 fadd.s1 fRes4L = fRes4L, fA0L
6176 fms.s1 FR_FracX = fSignifX, f1, FR_MHalf // y = x/8.0 - 1.5
6215 fma.s1 FR_rsq = FR_FracX, FR_FracX, f0 // y^2
6251 fma.s1 FR_Q4 = FR_rsq, FR_rsq, f0 // y^4
6262 fma.s1 fRes1H = FR_FracX, fA1, f0
6268 fma.s1 fA3 = fA3, FR_FracX, fA2 // v4
6274 fma.s1 fA5 = fA5, FR_FracX, fA4 // v5
6281 fma.s1 fA3L = FR_Q4, FR_Q4, f0 // v9 = y^8
6287 fma.s1 fA7 = fA7, FR_FracX, fA6 // v7
6293 fma.s1 fA9 = fA9, FR_FracX, fA8 // v8
6299 fms.s1 fRes1L = FR_FracX, fA1, fRes1H
6304 fma.s1 fA11 = fA11, FR_FracX, fA10 // v12
6310 fma.s1 fA13 = fA13, FR_FracX, fA12 // v13
6315 fma.s1 fRes2H = fRes1H, f1, fA0
6321 fma.s1 fA15 = fA15, FR_FracX, fA14 // v16
6326 fma.s1 fA5 = fA5, FR_rsq, fA3 // v3
6332 fma.s1 fA9 = fA9, FR_rsq, fA7 // v6
6338 fma.s1 fRes1L = FR_FracX, fA1L, fRes1L
6344 fms.s1 fRes2L = fA0, f1, fRes2H
6349 fma.s1 fA13 = fA13, FR_rsq, fA11 // v11
6355 fma.s1 fA9 = fA9, FR_Q4, fA5 // v2
6361 fma.s1 fRes1L = fRes1L, f1, fA0L
6367 fma.s1 fRes2L = fRes2L, f1, fRes1H
6372 fma.s1 fA15 = fA15, FR_Q4, fA13 // v10
6378 fma.s1 fRes2L = fRes1L, f1, fRes2L
6383 fma.s1 fPol = fA3L, fA15, fA9
6389 fma.s1 f8 = FR_rsq , fPol, fRes2H
6394 fma.s1 fPol = fPol, FR_rsq, f0
6400 fms.s1 fRes1L = fRes2H, f1, f8
6406 fma.s1 fRes1L = fRes1L, f1, fPol
6412 fma.s1 fRes1L = fRes1L, f1, fRes2L
6429 fms.s1 FR_FracX = fSignifX, f1, FR_MHalf
6453 fms.s1 FR_FracX = fSignifX, f1, FR_MHalf
6479 (p14) fms.s1 FR_FracX = f8, f1, FR_FracX
6485 (p15) fma.s1 FR_FracX = f8, f1, FR_FracX
6512 fms.s1 FR_FracX = f8, f1, fA5L
6524 fma.s1 FR_FracX = fB20, f1, f8 // 2.5 + x
6559 fma.s1 FR_FracX = FR_FracX, f1, f8
6615 fma.s1 fA4L = FR_FracX, FR_FracX, f0 // x^2
6659 fma.s1 fA5L = fA4L, fA4L, f0 // x^4
6664 fms.s1 fB2 = FR_FracX, FR_FracX, fA4L // x^2 - <x^2>
6671 fma.s1 fRes4H = fA3, FR_FracX, f0 // (A3*x)hi
6677 fma.s1 fA19 = fA19, FR_FracX, fA18 // D7*x + D6
6683 fma.s1 fResH = fA1, FR_FracX, f0 // (A1*x)hi
6688 fma.s1 fB6 = fA1L, FR_FracX, fA0L // A1L*x + A0L
6694 fma.s1 fA17 = fA17, FR_FracX, fA16 // D5*x + D4
6699 fma.s1 fA15 = fA15, FR_FracX, fA14 // D3*x + D2
6705 fma.s1 fA25 = fA25, FR_FracX, fA24 // C21*x + C20
6710 fma.s1 fA13 = fA13, FR_FracX, fA12 // D1*x + D0
6716 fma.s1 fA23 = fA23, FR_FracX, fA22 // C19*x + C18
6721 fma.s1 fA21 = fA21, FR_FracX, fA20 // C17*x + C16
6727 fms.s1 fRes4L = fA3, FR_FracX, fRes4H // delta((A3*x)hi)
6732 fadd.s1 fRes2H = fRes4H, fA2 // (A3*x + A2)hi
6738 fms.s1 fResL = fA1, FR_FracX, fResH // d(A1*x)
6743 fadd.s1 fRes1H = fResH, fA0 // (A1*x + A0)hi
6749 fma.s1 fA19 = fA19, fA4L, fA17 // Dhi
6754 fma.s1 fA11 = fA11, FR_FracX, fA10 // E7*x + E6
6767 fma.s1 fA15 = fA15, fA4L, fA13 // Dlo
6773 fma.s1 fA25 = fA25, fA4L, fA23
6779 fma.s1 fA9 = fA9, FR_FracX, fA8 // E5*x + E4
6784 fma.s1 fA7 = fA7, FR_FracX, fA6 // E3*x + E2
6790 fma.s1 fRes4L = fA3L, FR_FracX, fRes4L // (A3*x)lo
6795 fsub.s1 fRes2L = fA2, fRes2H
6801 fadd.s1 fResL = fResL, fB6 // (A1L*x + A0L) + d(A1*x)
6806 fsub.s1 fRes1L = fA0, fRes1H
6812 fma.s1 fA5 = fA5, FR_FracX, fA4 // E1*x + E0
6817 fma.s1 fB8 = fA5L, fA5L, f0 // x^8
6824 fma.s1 fA25 = fA25, fA4L, fA21
6829 fma.s1 fA19 = fA19, fA5L, fA15 // D
6835 fma.s1 fA11 = fA11, fA4L, fA9 // Ehi
6841 fadd.s1 fRes2L = fRes2L, fRes4H
6846 fadd.s1 fRes4L = fRes4L, fA2L // (A3*x)lo + A2L
6852 fma.s1 fRes3H = fRes2H, fA4L, f0 // ((A3*x + A2)*x^2)hi
6857 fadd.s1 fRes1L = fRes1L, fResH
6863 fma.s1 fRes3L = fRes2H, fB2, f0 // (A3*x + A2)hi*d(x^2)
6868 fma.s1 fA7 = fA7, fA4L, fA5 // Elo
6874 fma.s1 fA25 = fA25, fB8, fA19 // C*x^8 + D
6880 fadd.s1 fRes2L = fRes2L, fRes4L // (A3*x + A2)lo
6886 fms.s1 fB4 = fRes2H, fA4L, fRes3H // d((A3*x + A2)*x^2))
6891 fadd.s1 fRes1L = fRes1L, fResL // (A1*x + A0)lo
6897 fadd.s1 fB20 = fRes3H, fRes1H // Phi
6902 fma.s1 fA11 = fA11, fA5L, fA7 // E
6909 fma.s1 fRes3L = fRes2L, fA4L, fRes3L
6916 fadd.s1 fRes1L = fRes1L, fB4
6922 fsub.s1 fB18 = fRes1H, fB20
6927 fma.s1 fPol = fA25, fB8, fA11
6933 fadd.s1 fRes1L = fRes1L, fRes3L
6939 fadd.s1 fB18 = fB18, fRes3H
6944 fma.s1 fRes4H = fPol, fA5L, fB20
6950 fma.s1 fPolL = fPol, fA5L, f0
6956 fadd.s1 fB18 = fB18, fRes1L // Plo
6961 fsub.s1 fRes4L = fB20, fRes4H
6967 fadd.s1 fB18 = fB18, fPolL
6973 fadd.s1 fRes4L = fRes4L, fB18
6991 fma.s1 FR_FracX = fA5L, f1, f0 // x
6996 fma.s1 fB4 = fA5L, fA5L, f0 // x^2
7018 fma.s1 FR_FracX = fB4, f1, f0 // x
7023 fma.s1 fB4 = fB4, fB4, f0 // x^2
7058 fma.s1 fA5L = fB4, fB4, f0 // x^4
7063 fms.s1 fB2 = FR_FracX, FR_FracX, fB4 // x^2 - <x^2>
7103 fma.s1 fA0 = fA5L, fA5L, f0 // x^8
7108 fma.s1 fA0L = fA5L, FR_FracX, f0 // x^5
7126 fma.s1 fRes4H = fA4, fB4, f0 // A4*<x^2>
7131 fma.s1 fPol = fA24, FR_FracX, fA23 // C19*x + C18
7138 fma.s1 fRes1H = fA3, fB4, f0 // A3*<x^2>
7144 fma.s1 fA1L = fA3, fB2,fA1L // A3*d(x^2) + A1L
7150 fma.s1 fA20 = fA20, FR_FracX, fA19 // D7*x + D6
7155 fma.s1 fA18 = fA18, FR_FracX, fA17 // D5*x + D4
7161 fma.s1 fA16 = fA16, FR_FracX, fA15 // D3*x + D2
7166 fma.s1 fA14 = fA14, FR_FracX, fA13 // D1*x + D0
7172 fma.s1 fA2L = fA4, fB2,fA2L // A4*d(x^2) + A2L
7177 fma.s1 fA12 = fA12, FR_FracX, fA11 // E7*x + E6
7183 fms.s1 fRes2L = fA4, fB4, fRes4H // delta(A4*<x^2>)
7188 fadd.s1 fRes2H = fRes4H, fA2 // A4*<x^2> + A2
7194 fms.s1 fRes3L = fA3, fB4, fRes1H // delta(A3*<x^2>)
7199 fadd.s1 fRes3H = fRes1H, fA1 // A3*<x^2> + A1
7205 fma.s1 fA20 = fA20, fB4, fA18 // (D7*x + D6)*x^2 + D5*x + D4
7210 fma.s1 fA22 = fA22, FR_FracX, fA21 // C17*x + C16
7216 fma.s1 fA16 = fA16, fB4, fA14 // (D3*x + D2)*x^2 + D1*x + D0
7221 fma.s1 fPol = fA25, fB4, fPol // C20*x^2 + C19*x + C18
7227 fma.s1 fA2L = fA4L, fB4, fA2L // A4L*<x^2> + A4*d(x^2) + A2L
7232 fma.s1 fA1L = fA3L, fB4, fA1L // A3L*<x^2> + A3*d(x^2) + A1L
7238 fsub.s1 fRes4L = fA2, fRes2H // d1
7243 fma.s1 fResH = fRes2H, fB4, f0 // (A4*<x^2> + A2)*x^2
7249 fsub.s1 fRes1L = fA1, fRes3H // d1
7254 fma.s1 fB6 = fRes3H, FR_FracX, f0 // (A3*<x^2> + A1)*x
7260 fma.s1 fA10 = fA10, FR_FracX, fA9 // E5*x + E4
7265 fma.s1 fA8 = fA8, FR_FracX, fA7 // E3*x + E2
7272 fma.s1 fPol = fPol, fB4, fA22
7277 fma.s1 fA6 = fA6, FR_FracX, fA5 // E1*x + E0
7284 fadd.s1 fRes2L = fA2L, fRes2L
7290 fadd.s1 fRes3L = fA1L, fRes3L
7296 fadd.s1 fRes4L = fRes4L, fRes4H // d2
7301 fms.s1 fResL = fRes2H, fB4, fResH // d(A4*<x^2> + A2)*x^2)
7307 fadd.s1 fRes1L = fRes1L, fRes1H // d2
7312 fms.s1 fB8 = fRes3H, FR_FracX, fB6 // d((A3*<x^2> + A1)*x)
7318 fadd.s1 fB10 = fResH, fB6 // (A4*x^4 + .. + A1*x)hi
7323 fma.s1 fA12 = fA12, fB4, fA10 // Ehi
7330 fma.s1 fA20 = fA20, fA5L, fA16
7335 fma.s1 fA8 = fA8, fB4, fA6 // Elo
7341 fadd.s1 fRes2L = fRes2L, fRes4L // (A4*<x^2> + A2)lo
7347 fma.s1 fResL = fRes2H, fB2, fResL
7353 fadd.s1 fRes3L = fRes3L, fRes1L // (A4*<x^2> + A2)lo
7359 fsub.s1 fB12 = fB6, fB10
7365 fma.s1 fPol = fPol, fA0, fA20 // PolC*x^8 + PolD
7370 fma.s1 fPolL = fA12, fA5L, fA8 // E
7376 fma.s1 fResL = fB4, fRes2L, fResL // ((A4*<x^2> + A2)*x^2)lo
7382 fma.s1 fRes3L = fRes3L, FR_FracX, fB8 // ((A3*<x^2> + A1)*x)lo
7388 fadd.s1 fB12 = fB12, fResH
7394 fma.s1 fPol = fPol, fA0, fPolL
7400 fadd.s1 fRes3L = fRes3L, fResL
7406 fma.s1 fRes2H = fPol, fA0L, fB10
7412 fadd.s1 fRes3L = fB12, fRes3L
7418 fsub.s1 fRes4L = fB10, fRes2H
7424 fma.s1 fRes4L = fPol, fA0L, fRes4L
7430 fadd.s1 fRes4L = fRes4L, fRes3L
7525 fnorm.s1 f8 = f8 // normalize input value
7536 fcvt.fx.s1 fXint = f8 // Convert arg to int (int repres. in FR)
7542 fcmp.lt.s1 p15, p14 = f8, f0
7596 fma.s1 FR_X = f0,f0,f8