Lines Matching refs:fp0
307 set EXC_FP0, EXC_FPREGS+(0*12) # offset of saved fp0
567 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
574 fmov.s 0x8(%a6),%fp0 # load sgl input
575 fmov.x %fp0,FP_SRC(%a6)
624 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
631 fmov.d 0x8(%a6),%fp0 # load dbl input
632 fmov.x %fp0,FP_SRC(%a6)
682 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
744 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
751 fmov.s 0x8(%a6),%fp0 # load sgl input
752 fmov.x %fp0,FP_SRC(%a6)
801 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
808 fmov.d 0x8(%a6),%fp0 # load dbl input
809 fmov.x %fp0,FP_SRC(%a6)
859 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
921 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
928 fmov.s 0x8(%a6),%fp0 # load sgl input
929 fmov.x %fp0,FP_SRC(%a6)
978 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
985 fmov.d 0x8(%a6),%fp0 # load dbl input
986 fmov.x %fp0,FP_SRC(%a6)
1036 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1098 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1105 fmov.s 0x8(%a6),%fp0 # load sgl input
1106 fmov.x %fp0,FP_SRC(%a6)
1155 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1162 fmov.d 0x8(%a6),%fp0 # load dbl input
1163 fmov.x %fp0,FP_SRC(%a6)
1213 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1275 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1282 fmov.s 0x8(%a6),%fp0 # load sgl input
1283 fmov.x %fp0,FP_SRC(%a6)
1332 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1339 fmov.d 0x8(%a6),%fp0 # load dbl input
1340 fmov.x %fp0,FP_SRC(%a6)
1390 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1452 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1459 fmov.s 0x8(%a6),%fp0 # load sgl input
1460 fmov.x %fp0,FP_SRC(%a6)
1509 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1516 fmov.d 0x8(%a6),%fp0 # load dbl input
1517 fmov.x %fp0,FP_SRC(%a6)
1567 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1629 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1636 fmov.s 0x8(%a6),%fp0 # load sgl input
1637 fmov.x %fp0,FP_SRC(%a6)
1686 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1693 fmov.d 0x8(%a6),%fp0 # load dbl input
1694 fmov.x %fp0,FP_SRC(%a6)
1744 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1806 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1813 fmov.s 0x8(%a6),%fp0 # load sgl input
1814 fmov.x %fp0,FP_SRC(%a6)
1863 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1870 fmov.d 0x8(%a6),%fp0 # load dbl input
1871 fmov.x %fp0,FP_SRC(%a6)
1921 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1983 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
1990 fmov.s 0x8(%a6),%fp0 # load sgl input
1991 fmov.x %fp0,FP_SRC(%a6)
2040 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2047 fmov.d 0x8(%a6),%fp0 # load dbl input
2048 fmov.x %fp0,FP_SRC(%a6)
2098 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2160 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2167 fmov.s 0x8(%a6),%fp0 # load sgl input
2168 fmov.x %fp0,FP_SRC(%a6)
2217 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2224 fmov.d 0x8(%a6),%fp0 # load dbl input
2225 fmov.x %fp0,FP_SRC(%a6)
2275 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2337 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2344 fmov.s 0x8(%a6),%fp0 # load sgl input
2345 fmov.x %fp0,FP_SRC(%a6)
2394 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2401 fmov.d 0x8(%a6),%fp0 # load dbl input
2402 fmov.x %fp0,FP_SRC(%a6)
2452 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2514 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2521 fmov.s 0x8(%a6),%fp0 # load sgl input
2522 fmov.x %fp0,FP_SRC(%a6)
2571 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2578 fmov.d 0x8(%a6),%fp0 # load dbl input
2579 fmov.x %fp0,FP_SRC(%a6)
2629 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2691 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2698 fmov.s 0x8(%a6),%fp0 # load sgl input
2699 fmov.x %fp0,FP_SRC(%a6)
2748 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2755 fmov.d 0x8(%a6),%fp0 # load dbl input
2756 fmov.x %fp0,FP_SRC(%a6)
2806 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2868 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2875 fmov.s 0x8(%a6),%fp0 # load sgl input
2876 fmov.x %fp0,FP_SRC(%a6)
2925 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
2932 fmov.d 0x8(%a6),%fp0 # load dbl input
2933 fmov.x %fp0,FP_SRC(%a6)
2983 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3045 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3052 fmov.s 0x8(%a6),%fp0 # load sgl input
3053 fmov.x %fp0,FP_SRC(%a6)
3102 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3109 fmov.d 0x8(%a6),%fp0 # load dbl input
3110 fmov.x %fp0,FP_SRC(%a6)
3160 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3222 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3229 fmov.s 0x8(%a6),%fp0 # load sgl input
3230 fmov.x %fp0,FP_SRC(%a6)
3279 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3286 fmov.d 0x8(%a6),%fp0 # load dbl input
3287 fmov.x %fp0,FP_SRC(%a6)
3337 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3399 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3406 fmov.s 0x8(%a6),%fp0 # load sgl input
3407 fmov.x %fp0,FP_SRC(%a6)
3456 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3463 fmov.d 0x8(%a6),%fp0 # load dbl input
3464 fmov.x %fp0,FP_SRC(%a6)
3514 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3576 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3583 fmov.s 0x8(%a6),%fp0 # load sgl input
3584 fmov.x %fp0,FP_SRC(%a6)
3633 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3640 fmov.d 0x8(%a6),%fp0 # load dbl input
3641 fmov.x %fp0,FP_SRC(%a6)
3691 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3753 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3760 fmov.s 0x8(%a6),%fp0 # load sgl input
3761 fmov.x %fp0,FP_SRC(%a6)
3810 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3817 fmov.d 0x8(%a6),%fp0 # load dbl input
3818 fmov.x %fp0,FP_SRC(%a6)
3868 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3930 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3937 fmov.s 0x8(%a6),%fp0 # load sgl input
3938 fmov.x %fp0,FP_SRC(%a6)
3987 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
3994 fmov.d 0x8(%a6),%fp0 # load dbl input
3995 fmov.x %fp0,FP_SRC(%a6)
4045 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
4107 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
4114 fmov.s 0x8(%a6),%fp0 # load sgl input
4115 fmov.x %fp0,FP_SRC(%a6)
4154 fmovm.x &0x03,-(%sp) # store off fp0/fp1
4155 fmovm.x (%sp)+,&0x40 # fp0 now in fp1
4156 fmovm.x (%sp)+,&0x80 # fp1 now in fp0
4166 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
4173 fmov.d 0x8(%a6),%fp0 # load dbl input
4174 fmov.x %fp0,FP_SRC(%a6)
4214 fmovm.x &0x03,-(%sp) # store off fp0/fp1
4215 fmovm.x (%sp)+,&0x40 # fp0 now in fp1
4216 fmovm.x (%sp)+,&0x80 # fp1 now in fp0
4226 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
4274 fmovm.x &0x03,-(%sp) # store off fp0/fp1
4275 fmovm.x (%sp)+,&0x40 # fp0 now in fp1
4276 fmovm.x (%sp)+,&0x80 # fp1 now in fp0
4290 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
4297 fmov.s 0x8(%a6),%fp0 # load sgl dst
4298 fmov.x %fp0,FP_DST(%a6)
4303 fmov.s 0xc(%a6),%fp0 # load sgl src
4304 fmov.x %fp0,FP_SRC(%a6)
4356 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
4363 fmov.d 0x8(%a6),%fp0 # load dbl dst
4364 fmov.x %fp0,FP_DST(%a6)
4369 fmov.d 0x10(%a6),%fp0 # load dbl src
4370 fmov.x %fp0,FP_SRC(%a6)
4422 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
4494 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
4501 fmov.s 0x8(%a6),%fp0 # load sgl dst
4502 fmov.x %fp0,FP_DST(%a6)
4507 fmov.s 0xc(%a6),%fp0 # load sgl src
4508 fmov.x %fp0,FP_SRC(%a6)
4560 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
4567 fmov.d 0x8(%a6),%fp0 # load dbl dst
4568 fmov.x %fp0,FP_DST(%a6)
4573 fmov.d 0x10(%a6),%fp0 # load dbl src
4574 fmov.x %fp0,FP_SRC(%a6)
4626 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
4698 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
4705 fmov.s 0x8(%a6),%fp0 # load sgl dst
4706 fmov.x %fp0,FP_DST(%a6)
4711 fmov.s 0xc(%a6),%fp0 # load sgl src
4712 fmov.x %fp0,FP_SRC(%a6)
4764 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
4771 fmov.d 0x8(%a6),%fp0 # load dbl dst
4772 fmov.x %fp0,FP_DST(%a6)
4777 fmov.d 0x10(%a6),%fp0 # load dbl src
4778 fmov.x %fp0,FP_SRC(%a6)
4830 fmovm.x &0xc0,EXC_FP0(%a6) # save fp0/fp1
4906 # fp0 = sin(X) or cos(X) #
4909 # fp0 = sin(X) #
5028 fmov.x (%a0),%fp0 # LOAD INPUT
5029 fmov.x %fp0,X(%a6) # save input at X
5048 fmov.x %fp0,%fp1
5061 fsub.x (%a1)+,%fp0 # X-Y1
5062 fsub.s (%a1),%fp0 # fp0 = R = (X-Y1)-Y2
5113 fmul.x %fp0,%fp2 # S(A2+T(A4+TA6))
5115 fmul.x X(%a6),%fp0 # R'*S
5139 fmul.x %fp0,%fp0 # FP0 IS S
5144 fmov.x %fp0,%fp1
5147 fmov.x %fp0,X(%a6) # X IS S
5345 fadd.s COSB1(%pc),%fp0 # B1+S(B2...)
5346 fmul.x SPRIME(%a6),%fp0 # S'(B1+S(B2+...))
5577 fsub.x %fp3,%fp0 # fp0 = A := R - P
5580 fmov.x %fp0,%fp3 # fp3 = A
5585 fadd.x %fp1,%fp0 # fp0 = R := A+a
5591 fsub.x %fp0,%fp3 # fp3 = A-R
5615 # fp0 = tan(X) #
5760 fmov.x (%a0),%fp0 # LOAD INPUT
5777 fmov.x %fp0,%fp1
5787 fsub.x (%a1)+,%fp0 # X-Y1
5789 fsub.s (%a1),%fp0 # FP0 IS R = (X-Y1)-Y2
5800 fmov.x %fp0,%fp1
5822 fmul.x %fp0,%fp2 # RS(P1+S(P2+SP3))
5826 fadd.x %fp2,%fp0 # R+RS(P1+S(P2+SP3))
5833 fdiv.x %fp1,%fp0 # last inst - possible exception set
5837 fmov.x %fp0,%fp1
5838 fmul.x %fp0,%fp0 # S = R*R
5843 fmul.x %fp0,%fp3 # SQ4
5844 fmul.x %fp0,%fp2 # SP3
5849 fmul.x %fp0,%fp3 # S(Q3+SQ4)
5850 fmul.x %fp0,%fp2 # S(P2+SP3)
5855 fmul.x %fp0,%fp3 # S(Q2+S(Q3+SQ4))
5856 fmul.x %fp0,%fp2 # S(P1+S(P2+SP3))
5861 fmul.x %fp3,%fp0 # S(Q1+S(Q2+S(Q3+SQ4)))
5864 fadd.s &0x3F800000,%fp0 # 1+S(Q1+...)
5872 fdiv.x (%sp)+,%fp0 # last inst - possible exception set
5882 fmov.x %fp0,-(%sp)
5885 fmov.x (%sp)+,%fp0 # last inst - posibble exception set
5918 ftest.x %fp0 # test sign of argument
5924 fadd.x FP_SCR0(%a6),%fp0 # high part of reduction is exact
5925 fmov.x %fp0,%fp1 # save high result in fp1
5926 fadd.x FP_SCR1(%a6),%fp0 # low part of reduction
5927 fsub.x %fp0,%fp1 # determine low component of result
5928 fadd.x FP_SCR1(%a6),%fp1 # fp0/fp1 are reduced argument.
5934 fmov.x %fp0,INARG(%a6) # +-2**K * F, 1 <= F < 2
5963 fmov.x %fp0,%fp2
6242 fmov.x (%a0),%fp0 # LOAD INPUT
6246 fmov.x %fp0,X(%a6)
6418 fdiv.x %fp0,%fp1 # FP1 IS -1/X
6422 fmov.x %fp1,%fp0 # FP0 IS X'
6426 fmov.x %fp0,%fp1
6439 fmul.x %fp0,%fp2 # Y*(C2+Z*C4)
6442 fmul.x X(%a6),%fp0 # X'*Y
6448 fadd.x X(%a6),%fp0
6457 fadd.x NPIBY2(%pc),%fp0
6461 fadd.x PPIBY2(%pc),%fp0
6470 fmov.x NPIBY2(%pc),%fp0
6472 fadd.x PTINY(%pc),%fp0
6476 fmov.x PPIBY2(%pc),%fp0
6478 fadd.x NTINY(%pc),%fp0
6495 # fp0 = arcsin(X) #
6524 fmov.x (%a0),%fp0 # LOAD INPUT
6678 fmov.x PIBY2(%pc),%fp0
6692 # fp0 = exp(X) or exp(X)-1 #
7121 fmov.x (%a0),%fp0 # load input from (a0)
7123 fmov.x %fp0,%fp1
7124 fmul.s &0x42B8AA3B,%fp0 # 64/log2 * X
7127 fmov.l %fp0,%d1 # N = int( X * 64/log2 )
7129 fmov.l %d1,%fp0 # convert to floating-format
7142 #--fp1,fp2 saved on the stack. fp0 is N, fp1 is X,
7144 fmov.x %fp0,%fp2
7145 fmul.s &0xBC317218,%fp0 # N * L1, L1 = lead(-log2/64)
7147 fadd.x %fp1,%fp0 # X + N*L1
7148 fadd.x %fp2,%fp0 # fp0 is R, reduced arg.
7156 fmov.x %fp0,%fp1
7176 fmul.x %fp0,%fp3 # fp3 IS R*S*(A2+S*A4)
7179 fadd.x %fp3,%fp0 # fp0 IS R+R*S*(A2+S*A4),
7182 fadd.x %fp2,%fp0 # fp0 is EXP(R) - 1
7188 fmul.x %fp1,%fp0 # 2^(J/64)*(Exp(R)-1)
7190 fadd.s (%a1),%fp0 # accurate 2^(J/64)
7192 fadd.x %fp1,%fp0 # 2^(J/64) + 2^(J/64)*...
7199 fmul.x ADJSCALE(%a6),%fp0
7203 fmul.x SCALE(%a6),%fp0 # multiply 2^(M)
7210 fadd.s &0x3F800000,%fp0 # 1+X in user mode
7218 fmov.x (%a0),%fp0 # load input from (a0)
7220 fmov.x %fp0,%fp1
7221 fmul.s &0x42B8AA3B,%fp0 # 64/log2 * X
7224 fmov.l %fp0,%d1 # N = int( X * 64/log2 )
7226 fmov.l %d1,%fp0 # convert to floating-format
7257 fmov.s &0x3F800000,%fp0
7260 fadd.s (%sp)+,%fp0
7286 fmov.x (%a0),%fp0 # load input from (a0)
7288 fmov.x %fp0,%fp1
7289 fmul.s &0x42B8AA3B,%fp0 # 64/log2 * X
7291 fmov.l %fp0,%d1 # N = int( X * 64/log2 )
7293 fmov.l %d1,%fp0 # convert to floating-format
7304 #--fp1,fp2 saved on the stack. fp0 is N, fp1 is X,
7306 fmov.x %fp0,%fp2
7307 fmul.s &0xBC317218,%fp0 # N * L1, L1 = lead(-log2/64)
7309 fadd.x %fp1,%fp0 # X + N*L1
7310 fadd.x %fp2,%fp0 # fp0 is R, reduced arg.
7319 fmov.x %fp0,%fp1
7349 fmul.x %fp0,%fp2 # fp2 IS R*S*(A2+S*(A4+S*A6))
7350 fadd.x %fp1,%fp0 # fp0 IS R+S*(A1+S*(A3+S*A5))
7352 fadd.x %fp2,%fp0 # fp0 IS EXP(R)-1
7359 fmul.x (%a1),%fp0 # 2^(J/64)*(Exp(R)-1)
7369 fadd.x %fp1,%fp0 # p+(t+OnebySc), fp1 released
7370 fadd.x (%a1),%fp0 # T+(p+(t+OnebySc))
7378 fadd.s 12(%a1),%fp0 # p+t
7379 fadd.x (%a1),%fp0 # T+(p+t)
7380 fadd.x ONEBYSC(%a6),%fp0 # OnebySc + (T+(p+t))
7385 fadd.s (%a1),%fp0 # fp0 is p+t
7387 fadd.x %fp1,%fp0 # (T+OnebySc)+(p+t)
7392 fmul.x SC(%a6),%fp0
7408 fmov.x (%a0),%fp0
7411 fadd.x SC(%a6),%fp0
7416 fmov.x (%a0),%fp0
7417 fmul.d TWO140(%pc),%fp0
7421 fadd.x SC(%a6),%fp0
7424 fmul.d TWON140(%pc),%fp0
7429 fmov.x (%a0),%fp0 # fp0 is X
7430 fmul.x %fp0,%fp0 # fp0 is S := X*X
7433 fmul.x %fp0,%fp1 # fp1 is S*B12
7437 fmul.x %fp0,%fp2 # fp2 is S*B11
7438 fmul.x %fp0,%fp1 # fp1 is S*(B10 + ...
7443 fmul.x %fp0,%fp2 # fp2 is S*(B9+...
7444 fmul.x %fp0,%fp1 # fp1 is S*(B8+...
7449 fmul.x %fp0,%fp2 # fp2 is S*(B7+...
7450 fmul.x %fp0,%fp1 # fp1 is S*(B6+...
7455 fmul.x %fp0,%fp2 # fp2 is S*(B5+...
7456 fmul.x %fp0,%fp1 # fp1 is S*(B4+...
7461 fmul.x %fp0,%fp2 # fp2 is S*(B3+...
7462 fmul.x %fp0,%fp1 # fp1 is S*(B2+...
7464 fmul.x %fp0,%fp2 # fp2 is S*S*(B3+...)
7467 fmul.s &0x3F000000,%fp0 # fp0 is S*B1
7472 fadd.x %fp1,%fp0 # fp0 is S*B1+Q
7475 fadd.x (%a0),%fp0
7484 fmov.s &0xBF800000,%fp0 # fp0 is -1
7486 fadd.s &0x00800000,%fp0 # -1 + 2^(-126)
7498 # returned as an extended precision number in fp0. #
7503 # an exponent of $3fff and is returned in fp0. The range of #
7511 # fp0 = exponent(X) or mantissa(X) #
7520 fmov.w %d0,%fp0 # return exp in fp0
7573 # fp0 = cosh(X) #
7615 fmov.x (%a0),%fp0 # LOAD INPUT
7626 fabs.x %fp0 # |X|
7634 fmul.s &0x3F000000,%fp0 # (1/2)EXP(|X|)
7638 fdiv.x %fp0,%fp1 # 1/(2 EXP(|X|))
7642 fadd.x %fp1,%fp0
7649 fabs.x %fp0
7650 fsub.d T1(%pc),%fp0 # (|X|-16381LOG2_LEAD)
7651 fsub.d T2(%pc),%fp0 # |X| - 16381 LOG2, ACCURATE
7655 fmovm.x &0x01,-(%sp) # save fp0 to stack
7656 lea (%sp),%a0 # pass ptr to fp0
7658 add.l &0xc,%sp # clear fp0 from stack
7663 fmul.x TWO16380(%pc),%fp0
7672 fmov.s &0x3F800000,%fp0
7675 fadd.s &0x00800000,%fp0
7687 # fp0 = sinh(X) #
7727 fmov.x (%a0),%fp0 # LOAD INPUT
7739 fabs.x %fp0 # Y = |X|
7750 fmov.x %fp0,%fp1
7752 fmov.x %fp0,-(%sp)
7753 fdiv.x %fp1,%fp0 # Z/(1+Z)
7757 fadd.x (%sp)+,%fp0
7762 fmul.s (%sp)+,%fp0 # last fp inst - possible exceptions set
7768 fabs.x %fp0
7769 fsub.d T1(%pc),%fp0 # (|X|-16381LOG2_LEAD)
7776 fsub.d T2(%pc),%fp0 # |X| - 16381 LOG2, ACCURATE
7780 fmovm.x &0x01,-(%sp) # save fp0 on stack
7781 lea (%sp),%a0 # pass ptr to fp0
7783 add.l &0xc,%sp # clear fp0 from stack
7788 fmul.x (%sp)+,%fp0 # possible exception
7805 # fp0 = tanh(X) #
7852 fmov.x (%a0),%fp0 # LOAD INPUT
7854 fmov.x %fp0,X(%a6)
7873 fmov.x X(%a6),%fp0 # FP0 IS Y = 2|X|
7883 fmov.x %fp0,%fp1
7890 fdiv.x V(%a6),%fp0
7911 fmov.x X(%a6),%fp0 # Y = 2|X|
7921 fadd.s &0x3F800000,%fp0 # EXP(Y)+1
7925 fdiv.x %fp0,%fp1 # -SIGN(X)2 / [EXP(Y)+1 ]
7929 fmov.s %d1,%fp0 # SGN IN SGL FMT
7933 fadd.x %fp1,%fp0
7939 fmov.x X(%a6),%fp0 # last inst - possible exception set
7947 fmov.s %d1,%fp0
7952 fadd.s %d1,%fp0
7971 # fp0 = log(X) or log(1+X) #
8214 fmov.x (%a0),%fp0 # LOAD INPUT
8268 fmov.x X(%a6),%fp0
8271 fsub.x F(%a6),%fp0 # Y-F
8278 fmul.x (%a0),%fp0 # FP0 IS U = (Y-F)/F
8280 fmov.x %fp0,%fp2
8307 fmul.x %fp0,%fp1 # U*V*(A2+V*(A4+V*A6))
8308 fadd.x %fp2,%fp0 # U+V*(A1+V*(A3+V*A5))
8312 fadd.x %fp1,%fp0 # FP0 IS LOG(F) + LOG(1+U)
8315 fadd.x KLOG2(%a6),%fp0 # FINAL ADD
8442 fmov.x (%a0),%fp0 # LOAD INPUT
8443 fabs.x %fp0 # test magnitude
8444 fcmp.x %fp0,LTHOLD(%pc) # compare with min threshold
8448 fmov.x (%a0),%fp0 # return signed argument
8452 fmov.x (%a0),%fp0 # LOAD INPUT
8454 fmov.x %fp0,%fp1 # FP1 IS INPUT Z
8455 fadd.s one(%pc),%fp0 # X := ROUND(1+Z)
8456 fmov.x %fp0,X(%a6)
8480 fadd.s one(%pc),%fp0 # FP0 IS 1+X
8501 fmov.s TWO(%pc),%fp0
8504 fsub.x F(%a6),%fp0 # 2-F
8512 fadd.x %fp1,%fp0 # FP0 IS Y-F = (2-F)+2Z
8519 fmov.s one(%pc),%fp0
8522 fsub.x F(%a6),%fp0 # 1-F
8528 fadd.x %fp1,%fp0 # FP0 IS Y-F
8540 fmov.s negone(%pc),%fp0
8546 fmov.s zero(%pc),%fp0
8566 # fp0 = arctanh(X) #
8610 fabs.x (%a0),%fp0 # Y = |X|
8611 fmov.x %fp0,%fp1
8613 fadd.x %fp0,%fp0 # 2Y
8615 fdiv.x %fp1,%fp0 # 2Y/(1-Y)
8631 fmul.s (%sp)+,%fp0
8635 fabs.x (%a0),%fp0 # |X|
8636 fcmp.s %fp0,&0x3F800000
8656 # fp0 = log_10(X) or log_2(X) #
8739 fmov.b &0x1,%fp0
8740 fcmp.x %fp0,(%a0) # if operand == 1,
8749 fmul.x INV_L10(%pc),%fp0
8761 fmul.x INV_L10(%pc),%fp0
8783 fmov.l %d1,%fp0
8791 fmul.x INV_L2(%pc),%fp0
8806 fmul.x INV_L2(%pc),%fp0
8820 # fp0 = 2**X or 10**X #
8990 fmov.x %fp0,X(%a6)
9005 fmov.x %fp0,%fp1
9035 fsub.x %fp1,%fp0 # X - (1/64)*INT(64 X)
9041 fmul.x LOG2(%pc),%fp0 # FP0 IS R
9054 fadd.s &0x3F800000,%fp0 # RETURN 1 + X
9087 fmov.x %fp0,X(%a6)
9102 fmov.x %fp0,%fp1
9134 fsub.x %fp1,%fp0 # X - N L_LEAD
9137 fsub.x %fp2,%fp0 # X - N L_TRAIL
9143 fmul.x LOG10(%pc),%fp0 # FP0 IS R
9153 fmov.x %fp0,%fp1
9169 fmul.x %fp0,%fp3 # FP3 IS R*S*(A2+S*A4)
9172 fadd.x %fp3,%fp0 # FP0 IS R+R*S*(A2+S*A4)
9173 fadd.x %fp2,%fp0 # FP0 IS EXP(R) - 1
9180 fmul.x FACT1(%a6),%fp0
9181 fadd.x FACT2(%a6),%fp0
9182 fadd.x FACT1(%a6),%fp0
9190 fmul.x ADJFACT(%a6),%fp0 # FINAL ADJUSTMENT
9295 fmovm.x DST(%a1),&0x80 # load fp0 with normalized src
9308 fmul.x (%sp)+,%fp0 # do the multiply
9333 fmov.x DST(%a1),%fp0 # simply return dest
9353 # fp0 = FREM(X,Y) or FMOD(X,Y) #
9629 fmov.x R(%a6),%fp0 # no exception
9640 fmov.x R(%a6),%fp0
9668 fsub.x Y(%a6),%fp0 # no exceptions
9676 fneg.x %fp0
9702 fmul.x Scale(%pc),%fp0 # may cause underflow
10073 fmul.x pls_huge(%pc),%fp0
10083 fmul.x pls_huge(%pc),%fp0
10095 fmul.x pls_huge(%pc),%fp0
10115 # fp0 = default underflow or overflow result #
10118 # fp0 = default result #
10149 # fp0 = default result #
10152 # fp0 = default result #