78 #define d6 REG (d6)
306 | DOUBLE_FLOAT or LONG_FLOAT) in d6, and the last operation code in d5.
328 	movew	d6,a0@(FORMT)	| and __format
333 	cmpw	IMM (SINGLE_FLOAT),d6
335 	cmpl	IMM (SINGLE_FLOAT),d6
350 	clrl	d6
351 	movew	a0@(TRAPE),d6
352 	andl	d6,d7
613 	moveq	IMM (DOUBLE_FLOAT),d6
623 	moveq	IMM (DOUBLE_FLOAT),d6
631 	moveq	IMM (DOUBLE_FLOAT),d6
639 	moveq	IMM (DOUBLE_FLOAT),d6
648 	moveq	IMM (DOUBLE_FLOAT),d6
702 	movel	d2,d6		| save sign in d6 
708         swap	d6		| and also b's sign bit '
710 	andw	IMM (0x8000),d6	|
711 	orw	d6,d7		| and combine them into d7, so that a's sign '
713 				| low word, so d6 is free to be used
715 	andl	IMM (0x8000),d6
716 	orl	d6,d7
723 	movel	IMM (0x001fffff),d6 | mask for the fraction
727 	andl	d6,d0		| get fraction in d0
728 	notl	d6		| make d6 into mask for the exponent
729 	andl	d6,d4		| get exponent in d4
731 	cmpl	d6,d4		| check for INFINITY or NaN
743 	andl	d6,d5		| get exponent in d5
745 	cmpl	d6,d5		| check for INFINITY or NaN
747 	notl	d6		| make d6 into mask for the fraction again
748 	andl	d6,d2		| and get fraction in d2
764 | and d4-d5-d6-d7 for the second. To do this we store (temporarily) the
779 	movel	d7,d6		|
845 	roxrl	IMM (1),d6
849 	btst	IMM (0),d6
852 10:	lsrl	IMM (1),d6
855 	bset	IMM (31),d6
873 	movel	d6,d7
874 	movel	d5,d6
884 	movew	d6,d7
886 	movew	d5,d6
887 	swap	d6
902 	movew	d4,d6
903 	subw	d5,d6		| keep d5 (largest exponent) in d4
907 	movel	d5,d6
909 	movel	d6,d4
910 	subl	d5,d6
920 	cmpw	IMM (DBL_MANT_DIG+2),d6
922 	cmpl	IMM (DBL_MANT_DIG+2),d6
926 	cmpw	IMM (32),d6	| if difference >= 32, shift by longs
928 	cmpl	IMM (32),d6	| if difference >= 32, shift by longs
933 	cmpw	IMM (16),d6	| if difference >= 16, shift by words	
935 	cmpl	IMM (16),d6	| if difference >= 16, shift by words	
963 	dbra	d6,4b
965 	subql	IMM (1),d6
969 	movel	d7,d6
977 	subw	IMM (32),d6
979 	subl	IMM (32),d6
992 	subw	IMM (16),d6
994 	subl	IMM (16),d6
1016 	exg	d6,a3		| a3 is free to be used
1021 	movel	d6,a4
1022 	movel	a3,d6
1025 	movel	d7,d6		|
1027 	swap	d6              |
1028 	movew	IMM (0),d6	| and b's sign in d6 '
1029 	eorl	d7,d6		| compare the signs
1034 	exg	d6,a3		|
1039 	movel	d6,a4
1040 	movel	a3,d6
1044 	addxl	d6,d2		|
1092 	clrl	d6
1094 	movew	a1@(6),d6	| rounding mode in d6
1097 	cmpw	IMM (ROUND_TO_PLUS),d6
1099 	cmpl	IMM (ROUND_TO_PLUS),d6
1134 	exg	d6,a3		|
1139 	movel	d6,a4
1140 	movel	a3,d6
1144 	subxl	d6,d2		|
1200 	clrl	d6
1202 	movew	a1@(6),d6	| rounding mode in d6
1205 	cmpw	IMM (ROUND_TO_PLUS),d6
1207 	cmpl	IMM (ROUND_TO_PLUS),d6
1282 	notl	d6
1364 	movel	d2,d6		| 
1383 	eorl	d7,d6		| to check sign bits
1424 	movel	d7,d6			| another (mask for fraction)
1425 	notl	d6			|
1446 	andl	d6,d0			| isolate fraction
1457 	andl	d6,d2			|
1500 	movew	d3,d6		| get the rightmost 11 bits of d3
1501 	andw	IMM (0x07ff),d6	|
1502 	orw	d6,d2		| and put them into d2
1507 	movel	d3,d6		| get a copy of d3
1509 	andl	IMM (0xffe00000),d6 | get the top 11 bits of d3
1511 	lsrl	d7,d6
1512 	orl	d6,d2		| stick them at the end of d2
1515 	movel	d2,d6		| move b into d6-d7
1547 	addxl	d6,d6		|
1613 	moveq	IMM (29),d6
1616 	lsll	d6,d7
1620 	lsll	d6,d7
1624 	lsll	d6,d7
1707 	andl	d6,d0
1721 	andl	d6,d2
1756 	movel	d7,d6		| another (mask for fraction)
1757 	notl	d6		|
1782 	andl	d6,d0		| and isolate fraction
1793 	andl	d6,d2		|
1825 	movel	IMM (0),d6	| d6-d7 will hold the result
1826 	movel	d6,d7		| 
1836 	bset	d5,d6		| set the corresponding bit in d6
1916 	movel	d6,d0
1921 | 2^21 <= d6 < 2^23. Thus bit 23 is not set, but bit 22 could be set.
2021 	andl	d6,d0
2035 	andl	d6,d2
2064 	movel	IMM (0),d6	| use d6-d7 to collect bits flushed right
2065 	movel	d6,d7		| use d6-d7 to collect bits flushed right
2079 	roxrl	IMM (1),d6	| 
2085 	btst	IMM (0),d6
2088 13:	lsrl	IMM (1),d6
2091 	bset	IMM (31),d6
2110 2:	orl	d6,d2		| this is a trick so we don't lose  '
2117 	clrl	d6
2119 	movew	a1@(6),d6	| rounding mode in d6
2122 	cmpw	IMM (ROUND_TO_PLUS),d6
2124 	cmpl	IMM (ROUND_TO_PLUS),d6
2243 	movel	d0,d6		| copy signs into d6 (a) and d7(b)
2262 	eorl	d6,d7
2265 	tstl	d6
2270 	tstl	d6
2329 	bclr	IMM (31),d6
2348 	moveq	IMM (DOUBLE_FLOAT),d6
2521 	moveq	IMM (SINGLE_FLOAT),d6
2530 	moveq	IMM (SINGLE_FLOAT),d6
2537 	moveq	IMM (SINGLE_FLOAT),d6
2544 	moveq	IMM (SINGLE_FLOAT),d6
2552 	moveq	IMM (SINGLE_FLOAT),d6
2610 	movel	d0,d6		| save a to get exponent
2613 	andl	d4,d6		| get exponent in d6
2615 	cmpl	d4,d6		| check for INFINITY or NaN
2617 	swap	d6		| put exponent into first word
2620 | Now we have a's exponent in d6 (second byte) and the mantissa in d0. '
2643 | same, and put the largest exponent in d6. Note that we are using two
2647 	cmpw	d6,d7		| compare exponents
2649 	cmpl	d6,d7		| compare exponents
2654 	subl	d6,d7		| keep the largest exponent
2709 	exg	d6,d7		| exchange the exponents
2711 	eorl	d6,d7
2712 	eorl	d7,d6
2713 	eorl	d6,d7
2715 	subl	d6,d7		| keep the largest exponent
2769 | Now we have a in d0-d1, b in d2-d3, and the largest exponent in d6 (the
2775 	exg	d6,a0		| get signs back
2778 	movel	d6,d4
2779 	movel	a0,d6
2785 	eorl	d6,d7		| combine sign bits
2792 	exg	d6,a0		| now we have the exponent in d6
2794 	movel	d6,d4
2795 	movel	a0,d6
2807 	movel	d6,d2
2835 	clrl	d6
2837 	movew	a1@(6),d6	| rounding mode in d6
2840 	cmpw	IMM (ROUND_TO_PLUS),d6
2842 	cmpl	IMM (ROUND_TO_PLUS),d6
2871 	movel	d6,d7		| put sign in d7
2900 	clrl	d6
2902 	movew	a1@(6),d6	| rounding mode in d6
2905 	cmpw	IMM (ROUND_TO_PLUS),d6
2907 	cmpl	IMM (ROUND_TO_PLUS),d6
2958 	movel	d5,d6		| d5 contains 0x01000000
2959 	swap	d6
3087 	movel	IMM (INFINITY),d6	| useful constant (+INFINITY)
3088 	movel	d6,d5			| another (mask for fraction)
3097 	cmpl	d6,d0		| is a big?
3100 	cmpl	d6,d1		| now compare b with INFINITY
3105 	andl	d6,d2		| and isolate exponent in d2
3116 	andl	d6,d3		|
3142 	movel	d1,d6		| second operand in d6
3149 	lsll	IMM (31-FLT_MANT_DIG+1),d6		
3155 	lsll	IMM (1),d6	| get bit bn
3225 	cmpl	d6,d1		| is b NaN?
3301 	movel	IMM (INFINITY),d6	| useful constant (+INFINITY)
3302 	movel	d6,d5			| another (mask for fraction)
3311 	cmpl	d6,d0			| is a big?
3314 	cmpl	d6,d1			| now compare b with INFINITY 
3321 	andl	d6,d2		| and isolate exponent in d2
3332 	andl	d6,d3		|
3358 | d4, d5, d6 hold some constants
3359 	movel	d7,a0		| d6-d7 will hold the ratio of the fractions
3360 	movel	IMM (0),d6	| 
3361 	movel	d6,d7
3366 	bset	d3,d6		| set a bit in d6
3400 	movel	d6,d0		| put the ratio in d0-d1
3511 	movel	IMM (0),d6	| d6 is used temporarily
3523 	roxrl	IMM (1),d6	| d6 collect bits we would lose otherwise
3527 	lsrl	IMM (1),d6
3530 	bset	IMM (31),d6
3541 2:	orl	d6,d1		| this is a trick so we don't lose  '
3547 	clrl	d6
3549 	movew	a1@(6),d6	| rounding mode in d6
3552 	cmpw	IMM (ROUND_TO_PLUS),d6
3554 	cmpl	IMM (ROUND_TO_PLUS),d6
3671 	movel	d0,d6
3684 	eorl	d6,d7
3687 	tstl	d6
3692 	tstl	d6
3741 	bclr	IMM (31),d6
3750 	moveq	IMM (SINGLE_FLOAT),d6

Indexes created Wed Jul 03 21:38:27 MDT 2024