ifpsp060/src/fplsp.S

24 so long as this entire notice is retained without alteration in any modified and/or
26 No licenses are granted by implication, estoppel or otherwise under any patents
27 or trademarks of Motorola, Inc.
4906 #	fp0 = sin(X) or cos(X)						#
4923 #	2. If |X| >= 15Pi or |X| < 2**(-40), go to 7.			#
4926 #		k = N mod 4, so in particular, k = 0,1,2,or 3.		#
4951 #	1. If |X| >= 15Pi or |X| < 2**(-40), go to 6.			#
4954 #		k = N mod 4, so in particular, k = 0,1,2,or 3.		#
4959 #		j1 exclusive or with the l.s.b. of k.			#
5159 	or.l		&0x3F800000,%d1		# D0 IS SGN IN SINGLE
5203 # here, the operation may underflow iff the precision is sgl or dbl.
5381 	or.l		&0x3F800000,%d1
5490 	or.b		&0x80,FP_SCR0_EX(%a6)	# positive arg
5491 	or.b		&0x80,FP_SCR1_EX(%a6)
5543 	or.l		&0x5F000000,%d2		# d2 = SIGN(INARG)*2**63 IN SGL
5625 #	1. If |X| >= 15Pi or |X| < 2**(-40), go to 6.			#
5628 #		k = N mod 2, so in particular, k = 0 or 1.		#
5921 	or.b		&0x80,FP_SCR0_EX(%a6)	# positive arg
5922 	or.b		&0x80,FP_SCR1_EX(%a6)
5974 	or.l		&0x5F000000,%d2		# d2 = SIGN(INARG)*2**63 IN SGL
6054 #	Step 1. If |X| >= 16 or |X| < 1/16, go to Step 5.		#
6057 #		Note that k = -4, -3,..., or 3.				#
6282 	or.l		&0x04000000,XFRAC(%a6)	# SET 6-TH BIT TO 1
6310 	or.l		%d1,ATANF(%a6)		# ATANF IS NOW SIGN(F)*ATAN(|F|)
6569 	or.l		&0x3F800000,%d1		# +-1 IN SGL FORMAT
6660 	tst.b		(%a0)			# is X positive or negative?
6692 #	fp0 = exp(X) or exp(X)-1					#
6738 #			or 63.						#
6822 #		Hence, exp(X) may overflow or underflow or neither.	#
6826 #		or underflow. The inexact exception is not generated in	#
6863 #	Notes:	Exp(X) will surely overflow or underflow, depending on	#
6867 #		inexact together with either overflow or underflow.	#
6895 #			or 63.						#
7246 	tst.b		(%a0)			# is X positive or negative?
7343 	or.w		&0x8000,%d1		# signed/expo. of -2^(-M)
7511 #	fp0 = exponent(X) or mantissa(X)				#
7756 	or.l		&0x3F000000,%d1
7774 	or.l		&0x7FFB0000,%d1
7816 #	1. If |X| >= (5/2) log2 or |X| <= 2**(-40), go to 3.		#
7823 #	3. (|X| <= 2**(-40) or |X| >= (5/2) log2). If |X| < 1,		#
7928 	or.l		&0x3F800000,%d1		# SGN
7946 	or.l		&0x3F800000,%d1
7971 #	fp0 = log(X) or log(1+X)					#
8260 	or.l		&0x01000000,FFRAC(%a6)	# GET F: ATTACH A 1 AT THE EIGHTH BIT
8427 	or.l		%d7,%d4			# (D3,D4,D5) normalized
8496 	or.l		&0x01000000,FFRAC(%a6)	# F OBTAINED
8618 	or.l		&0x3F000000,%d1		# SIGN(X)*HALF
8656 #	fp0 = log_10(X) or log_2(X)					#
8820 #	fp0 = 2**X or 10**X						#
9076 	or.l		&0x00800001,%d1
9200 	or.l		&0x00800001,%d1
9207 #	    >= 2^14, an overflow or underflow is returned.		#
9292 # the src will force the dst to a DENORM value or worse. so, let's
9349 #	The input operands X and Y can be either normalized or		#
9353 #      fp0 = FREM(X,Y) or FMOD(X,Y)					#
9359 #                signQ := signX EOR signY. Record whether MOD or REM	#
9470 	or.l		%d7,%d4			# (D3,D4,D5) normalized
9515 	or.l		%d7,%d1			# (D0,D1,D2) normalized
9615 	or.l		%d7,%d1			# (D0,D1,D2) normalized
9686 	or.l		%d6,%d3			# sign and bits of Q
9690 #	or.l		%d3,%d6
9798 	bsr.l		unnorm_fix		# convert to norm,denorm,or zero
9848 	tst.b		SRC_EX(%a0)		# check sign for neg or pos
10104 #	t_catch(): Handle 060FPLSP OVFL,UNFL,or INEX2 exception during	#
10106 #	t_catch2(): Handle 060FPLSP OVFL,UNFL,or INEX2 exception during	#
10115 #	fp0 = default underflow or overflow result			#
10121 #	If an overflow or underflow occurred during the last		#
10131 	or.l		%d0,USER_FPSR(%a6)
10137 	or.l		%d0,USER_FPSR(%a6)
10223 	or.l		%d0,USER_FPSR(%a6)
10231 	or.l		%d0,USER_FPSR(%a6)
10529 # szr_inf(): Return +ZERO for a negative src operand or		#
10539 # sopr_inf(): Return +INF for a positive src operand or			#
10550 # setoxm1i(): Return minus one for a negative src operand or	#
10850 	or.l		%d3, %d0		# create hi(man)
10879 # unnorm_fix(): - changes an UNNORM to one of NORM, DENORM, or ZERO	#
10892 #	d0 = optype tag - is corrected to one of NORM, DENORM, or ZERO	#
10893 #	a0 = input operand has been converted to a norm, denorm, or	#
10929 	or.w		%d0, %d1		# {sgn,new exp}