ifpsp060/src/fpsp.S

1442 	cmpi.w		%d0,&0x3f80		# is |exp| == $3f80?
1444 	cmpi.w		%d0,&0x407f		# no; is |exp| == $407f?
2099 #	OPERR	: if ((k_factor > +17) || (dec. exp exceeds 3 digits))
6733 	mov.l		(%a0),%d1		# pack exp w/ upper 16 fraction
6994 #		because EXPM1 is intended to evaluate exp(X)-1		#
7018 #	Step 4.	Approximate exp(R)-1 by a polynomial			#
7025 #			|p - (exp(R)-1)| <	|R| * 2^(-72.7)		#
7046 #	Step 6.	Reconstruction of exp(X)-1				#
7047 #			exp(X)-1 = 2^M * ( 2^(J/64) + p - 2^(-M) ).	#
7058 #	Step 7.	exp(X)-1 for |X| < 1/4.					#
7062 #	Step 8.	Calculate exp(X)-1, |X| < 2^(-65).			#
7076 #	Step 9.	Calculate exp(X)-1, |X| < 1/4, by a polynomial		#
7083 #			|p - (exp(X)-1)| < |X| 2^(-70.6)		#
7095 #	Step 10. Calculate exp(X)-1 for |X| >= 70 log 2.		#
7096 #		10.1 If X >= 70log2 , exp(X) - 1 = exp(X) for all	#
7098 #		10.2 If X <= -70log2, exp(X) - 1 = -1 for all practical	#
7637 	neg.w		%d0			# new exp = -(shft amt)
7639 	fmov.w		%d0,%fp0		# return exp in fp0
7645 	mov.w		SRC_EX(%a0),%d0		# get the exp
7646 	ori.w		&0x7fff,%d0		# clear old exp
7647 	bclr		&0xe,%d0		# make it the new exp +-3fff
7706 #		cosh(X) := Fact * exp(Y').				#
7821 #             sinh(X) := sgnFact * exp(Y').				#
9585 	andi.l		&0x00007fff,%d1		# strip sign from dst exp
11638 	mov.w		FP_SCR0_EX(%a6),%d1	# load {sgn,exp}
11643 	or.w		%d2,%d1			# concat old sign,new exp
12280 	mov.w		FP_SCR0_EX(%a6),%d1	# fetch {sgn,exp}
13060 	mov.w		FP_SCR0_EX(%a6),%d1	# fetch {sgn,exp}
13067 	or.w		%d2,%d1			# concat sign,exp
13377 	mov.w		SRC_EX(%a0),FP_SCR0_EX(%a6) # copy sign, zero exp
13517 	mov.w		FP_SCR0_EX(%a6),%d1	# fetch old sign,exp
13564 	mov.w		FP_SCR0_EX(%a6),%d1	# load sgn,exp
13569 	or.w		%d1,%d2			# concat old sign,new exp
13628 	or.w		%d2,%d1			# concat new sign,new exp
13629 	mov.w		%d1,FP_SCR1_EX(%a6)	# insert new exp
14009 	mov.w		FP_SCR0_EX(%a6),%d1	# load {sgn,exp}
14014 	or.w		%d2,%d1			# concat old sign,new exp
14056 	mov.w		FP_SCR0_EX(%a6),%d1	# fetch {sgn,exp}
14063 	or.w		%d2,%d1			# concat old sign,new exp
14356 	mov.w		FP_SCR0_EX(%a6),%d1	# load {sgn,exp}
14361 	or.w		%d2,%d1			# concat old sign,new exp
14409 	mov.w		FP_SCR0_EX(%a6),%d1	# fetch {sgn,exp}
14416 	or.w		%d2,%d1			# concat old sign,new exp
14466 	mov.w		FP_SCR0_EX(%a6),%d1	# fetch {sgn,exp}
14473 	or.w		%d2,%d1			# concat old sign, new exp
14685 	mov.w		(%sp),%d2		# fetch new sign, exp
14699 	or.w		%d2,%d1			# concat sign,new exp
14749 	or.w		%d2,%d1			# concat sign,new exp
14819 	mov.w		FP_SCR0_EX(%a6),%d1	# fetch {sgn,exp}
14826 	or.w		%d2,%d1			# concat sign,new exp
15198 	mov.w		(%sp),%d1		# fetch {sgn,exp}
15202 	or.w		%d2,%d1			# concat sign,exp
15272 	mov.w		FP_SCR0_EX(%a6),%d1	# fetch {sgn,exp}
15279 	or.w		%d2,%d1			# concat sgn,exp
15710 	or.w		%d2,%d1			# concat new sign,new exp
15711 	mov.w		%d1,FP_SCR1_EX(%a6)	# insert new exp
15882 	cmp.w		%d0, %d1		# is src exp >= dst exp?
15885 # dst exp is >  src exp; scale dst to exp = 0x3fff
15894 	bsr.l		norm			# normalize the denorm; result is new exp
15895 	neg.w		%d0			# new exp = -(shft val)
15896 	mov.w		%d0,L_SCR1(%a6)		# inset new exp
15900 	subi.w		&mantissalen+2,%d0	# subtract mantissalen+2 from larger exp
15909 	or.w		%d1,%d0			# concat {sgn,new exp}
15917 	bset		&0x0,1+FP_SCR0_EX(%a6)	# set exp = 1
15922 # src exp is >= dst exp; scale src to exp = 0x3fff
15930 	bsr.l		norm			# normalize the denorm; result is new exp
15931 	neg.w		%d0			# new exp = -(shft val)
15932 	mov.w		%d0,2+L_SCR1(%a6)	# inset new exp
15936 	subi.w		&mantissalen+2,%d0	# subtract mantissalen+2 from larger exp
15945 	or.w		%d1,%d0			# concat {sgn,new exp}
15953 	bset		&0x0,1+FP_SCR1_EX(%a6)	# set exp = 1
16000 	sub.l		%d1,%d0			# scale = BIAS + (-exp)
16042 	mov.w		FP_SCR0_EX(%a6),%d1	# extract operand's {sgn,exp}
16047 	btst		&0x0,%d1		# is exp even or odd?
16061 	sub.l		%d1,%d0			# scale = BIAS + (-exp)
16069 	btst		&0x0,%d0		# is exp even or odd?
16127 	sub.l		%d1,%d0			# scale = BIAS + (-exp)