Lines Matching refs:exponent
415 set SGL_LO, 0x3f81 # min sgl prec exponent
416 set SGL_HI, 0x407e # max sgl prec exponent
417 set DBL_LO, 0x3c01 # min dbl prec exponent
418 set DBL_HI, 0x43fe # max dbl prec exponent
419 set EXT_LO, 0x0 # min ext prec exponent
420 set EXT_HI, 0x7ffe # max ext prec exponent
516 set sgl_thresh, 0x3f81 # minimum sgl exponent
517 set dbl_thresh, 0x3c01 # minimum dbl exponent
1440 mov.w LOCAL_EX(%a0),%d0 # fetch src exponent
1455 addi.w &0x3f81,%d0 # adjust new exponent
1456 andi.w &0x8000,LOCAL_EX(%a0) # clear old exponent
1457 or.w %d0,LOCAL_EX(%a0) # insert new exponent
1461 andi.w &0x8000,LOCAL_EX(%a0) # clear bogus exponent
1466 ori.w &0x7fff,LOCAL_EX(%a0) # make exponent = $7fff
1470 mov.w LOCAL_EX(%a0),%d0 # fetch src exponent
1487 addi.w &0x3c01,%d0 # adjust new exponent
1488 andi.w &0x8000,LOCAL_EX(%a0) # clear old exponent
1489 or.w %d0,LOCAL_EX(%a0) # insert new exponent
1517 mov.w FP_SRC_EX(%a6),%d0 # get exponent
2311 mov.w FP_SRC_EX(%a6),%d0 # fetch DENORM exponent
2316 neg.w %d0 # make exponent negative
2322 andi.w &0x8000,FP_SRC_EX(%a6) # clear old exponent
2323 ori.w &0x3f80,FP_SRC_EX(%a6) # insert new "skewed" exponent
2328 mov.w FP_SRC_EX(%a6),%d0 # fetch DENORM exponent
2336 mov.w %d0,FP_SRC_EX(%a6) # insert exponent with cleared sign
2341 mov.w &0x3c00,%d0 # new exponent
2347 mov.w %d0,FP_SRC_EX(%a6) # insert new exponent
2529 # The packed operand is an INF or a NAN if the exponent field is all ones.
3167 mov.w FP_SRC_EX(%a6),%d1 # fetch exponent
3431 ori.l &0x7fc00000,%d0 # insert new exponent,SNAN bit
3445 ori.l &0x7fc00000,%d0 # insert new exponent,SNAN bit
3458 ori.l &0x7ff80000,%d0 # insert new exponent,SNAN bit
3621 # w/ an exponent value of 0x401e. we convert this to extended precision here.
3624 cmpi.w FP_SRC_EX(%a6),&0x401e # is exponent 0x401e?
6828 # are the sign and biased exponent field of |X|; the #
7602 # sgetexp(): returns the exponent portion of the input argument. #
7603 # The exponent bias is removed and the exponent value is #
7609 # an exponent of $3fff and is returned in fp0. The range of #
7617 # fp0 = exponent(X) or mantissa(X) #
7623 mov.w SRC_EX(%a0),%d0 # get the exponent
7652 mov.w %d0,FP_SCR0_EX(%a6) # insert new exponent
7663 # then load the exponent with +/1 $3fff.
7667 bsr.l norm # normalize exponent
8491 #----below), adjusting exponent and storing -k to ADJK
8496 mov.l (%a0),%d3 # D3 is exponent of smallest norm. #
9583 mov.w DST_EX(%a1),%d1 # get dst exponent
9641 clr.l -(%sp) # make zero exponent
9648 clr.l -(%sp) # make zero exponent
9659 swap %d0 # put exponent in high word
9660 clr.l -(%sp) # insert new exponent
10209 # so, normalize the mantissa, add 0x6000 to the new exponent,
10220 or.w %d0,FP_SCR0_EX(%a6) # insert new exponent
11528 # scale_to_zero_src() - scale src exponent to zero #
11529 # scale_to_zero_dst() - scale dst exponent to zero #
11552 # result operand to the proper exponent. #
11596 bsr.l scale_to_zero_src # scale src exponent
11599 bsr.l scale_to_zero_dst # scale dst exponent
11618 # - scale the result exponent using the scale factor. if both operands were
11644 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
11659 # of this operation then has its exponent scaled by -0x6000 to create the
11695 # - if precision is extended, then we have the EXOP. simply bias the exponent
11716 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
11770 # of this operation then has its exponent scaled by -0x6000 to create the
11836 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
12025 # scale_to_zero_src() - scale src exponent to zero #
12043 # sgl/dbl, must scale exponent and perform an "fmove". Check to see #
12109 # exponent and insert back into the operand.
12117 neg.w %d0 # new exponent = -(shft val)
12118 addi.w &0x6000,%d0 # add new bias to exponent
12123 mov.w %d0,FP_SCR0_EX(%a6) # insert new exponent
12171 or.w %d1,%d2 # concat old sign,new exponent
12172 mov.w %d2,FP_SCR0_EX(%a6) # insert new exponent
12224 mov.w FP_SCR0_EX(%a6),%d1 # load current exponent
12234 mov.w %d2,FP_SCR1_EX(%a6) # insert new exponent
12288 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
12345 # scale_to_zero_src() - scale src exponent to zero #
12346 # scale_to_zero_dst() - scale dst exponent to zero #
12369 # result operand to the proper exponent. #
12380 long 0x3fff - 0x7ffe # ext overflow exponent
12381 long 0x3fff - 0x407e # sgl overflow exponent
12382 long 0x3fff - 0x43fe # dbl overflow exponent
12418 bsr.l scale_to_zero_src # scale src exponent
12421 bsr.l scale_to_zero_dst # scale dst exponent
12458 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
12488 mov.w (%sp),%d0 # fetch new exponent
12529 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
12800 # scale_to_zero_src() - scale sgl/dbl source exponent #
12820 # scale the result exponent and return result. FPSR gets set based on #
12892 # exponent and insert back into the operand.
12897 neg.w %d0 # new exponent = -(shft val)
12898 addi.w &0x6000,%d0 # add new bias to exponent
12902 or.w %d1,%d0 # concat old sign, new exponent
12903 mov.w %d0,FP_SCR0_EX(%a6) # insert new exponent
12952 mov.w %d2,FP_SCR0_EX(%a6) # insert new exponent
13004 mov.w FP_SCR0_EX(%a6),%d1 # load current exponent
13068 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
13419 # scale_to_zero_src() - make exponent. = 0; get scale factor #
13438 # scale the operand such that the exponent is zero. Perform an "fabs" #
13440 # exponent would take an exception. If so, use unf_res() or ovf_res() #
13443 # result exponent and return. #
13481 mov.w %d1,FP_SCR0_EX(%a6) # insert exponent
13499 mov.w %d0,FP_SCR0_EX(%a6) # insert exponent
13510 # exponent and insert back into the operand.
13515 neg.w %d0 # new exponent = -(shft val)
13516 addi.w &0x6000,%d0 # add new bias to exponent
13520 or.w %d1,%d0 # concat old sign, new exponent
13521 mov.w %d0,FP_SCR0_EX(%a6) # insert new exponent
13570 mov.w %d2,FP_SCR0_EX(%a6) # insert new exponent
13619 mov.w FP_SCR0_EX(%a6),%d1 # load current exponent
13683 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
13930 # scale_to_zero_src() - scale src exponent to zero #
13931 # scale_to_zero_dst() - scale dst exponent to zero #
13954 # result operand to the proper exponent. #
13978 bsr.l scale_to_zero_src # scale exponent
13981 bsr.l scale_to_zero_dst # scale dst exponent
14015 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
14064 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
14141 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
14271 # scale_to_zero_src() - scale src exponent to zero #
14272 # scale_to_zero_dst() - scale dst exponent to zero #
14295 # result operand to the proper exponent. #
14362 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
14381 mov.w (%sp),%d1 # fetch new exponent
14417 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
14474 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
14614 # scale_to_zero_src() - set src operand exponent equal to zero #
14615 # scale_to_zero_dst() - set dst operand exponent equal to zero #
14629 # occur. Then, check result exponent to see if exception would have #
14631 # the correct result exponent and return. Set FPSR bits as appropriate. #
14700 mov.w %d1,(%sp) # insert new exponent
14750 mov.w %d1,(%sp) # insert new exponent
14827 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
14859 # ok, so now the result has a exponent equal to the smallest normalized
14860 # exponent for the selected precision. also, the mantissa is equal to
15067 # scale_to_zero_src() - set src operand exponent equal to zero #
15068 # scale_to_zero_dst() - set dst operand exponent equal to zero #
15082 # occur. Then, check result exponent to see if exception would have #
15084 # the correct result exponent and return. Set FPSR bits as appropriate. #
15138 mov.w (%sp),%d2 # fetch new exponent
15152 or.w %d2,%d1 # insert new exponent
15153 mov.w %d1,(%sp) # insert new exponent
15203 mov.w %d1,(%sp) # insert new exponent
15280 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
15312 # ok, so now the result has a exponent equal to the smallest normalized
15313 # exponent for the selected precision. also, the mantissa is equal to
15529 # result operand to the proper exponent. #
15633 mov.w %d2,FP_SCR0_EX(%a6) # insert new exponent
15657 # the exponent is 3fff or 3ffe. if it's 3ffe, then it's a safe number
15660 btst &0x0,1+FP_SCR0_EX(%a6) # is exponent 0x3fff?
15701 mov.w FP_SCR0_EX(%a6),%d1 # load current exponent
15765 mov.w %d1,FP_SCR0_EX(%a6) # insert new exponent
15774 btst &0x0,1+FP_SCR0_EX(%a6) # is exponent 0x3fff?
15844 # norm() - normalize mantissa after adjusting exponent #
15856 # If the DST exponent is > the SRC exponent, set the DST exponent #
15857 # equal to 0x3fff and scale the SRC exponent by the value that the #
15858 # DST exponent was scaled by. If the SRC exponent is greater or equal, #
15861 # plus two, then set the smallest exponent to a very small value as a #
15879 mov.w %d0,L_SCR1(%a6) # store src exponent
15880 mov.w %d1,2+L_SCR1(%a6) # store dst exponent
15906 add.w 0x2(%sp),%d0 # scale src exponent by scale factor
15910 mov.w %d0,FP_SCR0_EX(%a6) # insert new dst exponent
15916 andi.w &0x8000,FP_SCR0_EX(%a6) # zero src exponent
15942 add.w 0x2(%sp),%d0 # scale dst exponent by scale factor
15946 mov.w %d0,FP_SCR1_EX(%a6) # insert new dst exponent
15952 andi.w &0x8000,FP_SCR1_EX(%a6) # zero dst exponent
15962 # scale_to_zero_src(): scale the exponent of extended precision #
15976 # Set the exponent of the input operand to 0x3fff. Save the value #
15977 # of the difference between the original and new exponent. Then, #
15988 andi.l &0x7fff,%d1 # extract operand's exponent
15991 or.w &0x3fff,%d0 # insert new operand's exponent(=0)
15993 mov.w %d0,FP_SCR0_EX(%a6) # insert biased exponent
16007 neg.l %d0 # new exponent = -(shft val)
16015 # scale_sqrt(): scale the input operand exponent so a subsequent #
16030 # If the exponent of the input operand is even, set the exponent #
16032 # exponent of the input operand is off, set the exponent to ox3fff and #
16043 andi.l &0x7fff,%d1 # extract operand's exponent
16050 ori.w &0x3fff,FP_SCR0_EX(%a6) # insert new operand's exponent(=0)
16058 ori.w &0x3ffe,FP_SCR0_EX(%a6) # insert new operand's exponent(=0)
16072 ori.w &0x3fff,FP_SCR0_EX(%a6) # insert new operand's exponent(=0)
16079 ori.w &0x3ffe,FP_SCR0_EX(%a6) # insert new operand's exponent(=0)
16089 # scale_to_zero_dst(): scale the exponent of extended precision #
16103 # Set the exponent of the input operand to 0x3fff. Save the value #
16104 # of the difference between the original and new exponent. Then, #
16115 andi.l &0x7fff,%d1 # extract operand's exponent
16118 or.w &0x3fff,%d0 # insert new operand's exponent(=0)
16120 mov.w %d0,FP_SCR1_EX(%a6) # insert biased exponent
16133 neg.l %d0 # new exponent = -(shft val)
20422 or.w %d0,FP_SCR0_EX(%a6) # insert new exponent
20444 mov.w SRC_EX(%a0),%d0 # extract exponent
20596 # (3) if exp still equals zero, then insert original exponent
20614 fcmp.b %fp0,&0x2 # did exponent increase?
20682 mov.w SRC_EX(%a0),%d0 # extract exponent
20807 # (3) if exp still equals zero, then insert original exponent
20825 fcmp.b %fp0,&0x2 # did exponent increase?
20869 mov.w FTEMP_EX(%a0),%d0 # get exponent
20934 mov.w FTEMP_EX(%a0),%d0 # get exponent
20999 # we zero the exponent
21003 # algorithm allows the exponent to be non-zero. the 881/2 do not. Therefore,
21004 # if the mantissa is zero, I will zero the exponent, too.
21691 # According to the exponent underflow threshold for the given #
21693 # exponent of the operand to the threshold value. While shifting the #
21703 # table of exponent threshold values for each precision
21713 # Load the exponent threshold for the precision selected and check
21714 # to see if (threshold - exponent) is > 65 in which case we can
21746 # dnrm_lp(): normalize exponent/mantissa to specified threshold #
21773 # exponent is.
21776 sub.w FTEMP_EX(%a0), %d1 # d1 = threshold - uns exponent
21817 mov.w %d0, FTEMP_EX(%a0) # exponent = denorm threshold
21869 mov.w %d0, FTEMP_EX(%a0) # exponent = denorm threshold
22134 add.w &0x1, FTEMP_EX(%a0) # and incr exponent
22177 addq.w &0x1, FTEMP_EX(%a0) # incr exponent
22312 # a0 = the input operand's mantissa is normalized; the exponent #
22375 # zero; both the exponent and mantissa are changed. #
22398 mov.w FTEMP_EX(%a0), %d1 # extract exponent
22405 # exponent would not go < 0. Therefore, number stays normalized
22407 sub.w %d0, %d1 # shift exponent value
22408 mov.w FTEMP_EX(%a0), %d0 # load old exponent
22411 mov.w %d1, FTEMP_EX(%a0) # insert new exponent
22419 # exponent would go < 0, so only denormalize until exp = 0
22423 bgt.b unnorm_nrm_zero_lrg # no; go handle large exponent
22458 and.w &0x8000, FTEMP_EX(%a0) # force exponent to zero
22478 # Simply test the exponent, j-bit, and mantissa values to #
22487 mov.w FTEMP_EX(%a0), %d0 # extract exponent
22498 tst.w %d0 # is exponent = 0?
22519 andi.w &0x8000,FTEMP_EX(%a0) # clear exponent
22558 # Simply test the exponent, j-bit, and mantissa values to #
22621 # Simply test the exponent, j-bit, and mantissa values to #
22688 # exponent is extended to 16 bits and the sign is stored in the unused #
22706 mov.w FTEMP_EX(%a0), %d1 # extract exponent
22709 mov.w %d1, FTEMP_EX(%a0) # insert 16 bit exponent
22770 mov.w FTEMP_EX(%a0),%d1 # extract exponent
22773 mov.w %d1,FTEMP_EX(%a0) # insert 16 bit exponent
22970 # The packed operand is an INF or a NAN if the exponent field is all ones.
23008 # A1. Convert the bcd exponent to binary by successive adds and #
23023 # exponent equal to the exponent from A1 and the zero count #
23041 # the exponent factor. This is done by multiplying the #
23043 # exponent sign is positive, and dividing FP0 by FP1 if #
23081 # Calculate exponent:
23083 # 2. Calculate absolute value of exponent in d1 by mul and add.
23084 # 3. Correct for exponent sign.
23092 # (*) d1: accumulator for binary exponent
23099 # (*) L_SCR1: copy of original exponent word
23142 # ( ) L_SCR1: copy of original exponent word
23192 # for the adjusted exponent. That number is subtracted from the exp
23194 # of this is to reduce the value of the exponent and the possibility
23197 # 1. Branch on the sign of the adjusted exponent.
23212 # *Why 27? If the adjusted exponent is within -28 < expA < 28, than
23214 # exponent towards zero. Since all pwrten constants with a power
23216 # attempt to lessen the resultant exponent.
23229 # ( ) L_SCR1: copy of original exponent word
23232 # First check the absolute value of the exponent to see if this
23233 # routine is necessary. If so, then check the sign of the exponent
23235 # This section handles a positive adjusted exponent.
23292 # This section handles a negative adjusted exponent.
23340 # Calculate power-of-ten factor from adjusted and shifted exponent.
23346 # ( ) d1: exponent
23353 # ( ) d1: exponent
23359 # Pwrten calculates the exponent factor in the selected rounding mode
23434 btst &30,(%a0) # test the sign of the exponent
23553 # A15. Convert the exponent to bcd. #
23556 # Test the length of the final exponent string. If the #
23623 mov.l (%a0),L_SCR2(%a6) # save exponent for sign check
23672 # d0: k-factor/exponent
23686 # F_SCR2:Abs(X)/Abs(X) with $3fff exponent
23696 mov.w &0x3fff,FP_SCR1(%a6) # replace exponent with 0x3fff
23734 # d0: exponent/Unchanged
23748 # F_SCR2:Abs(X) with $3fff exponent/Unchanged
23806 # d0: exponent/scratch - final is 0
23820 # F_SCR2:Abs(X) with $3fff exponent/Unchanged
23923 # F_SCR2:Abs(X) with $3fff exponent/Unchanged
23947 mov.w (%sp),%d3 # grab exponent
23960 or.w %d3,(%sp) # insert new exponent
24053 # F_SCR2:Y adjusted for inex/Y with original exponent
24090 mov.l L_SCR2(%a6),FP_SCR1(%a6) # restore original exponent
24120 # F_SCR2:Y with original exponent/Unchanged
24229 # F_SCR2:Y with original exponent/Unchanged
24242 mov.l (%a0),%d0 # move exponent to d0
24243 swap %d0 # put exponent in lower word
24269 # A15. Convert the exponent to bcd.
24287 # d1: x/scratch (0);shift count for final exponent packing
24301 # F_SCR2:Y with original exponent/ILOG/10^4
24316 fmov.s F4933(%pc),%fp0 # force exponent to 4933
24324 fbneq.w not_zero # if zero, force exponent
24325 fmov.s FONE(%pc),%fp0 # force exponent to 1