ifpsp060/src/fpsp.S

447 set ovfl_bit,		4			# overflow
454 set aovfl_bit,		6			# accrued overflow bit
475 set ovfl_mask,		0x00001000		# overflow exception mask
482 set aovfl_mask,		0x00000040		# accrued overflow
612 #	On the 060, if an FP overflow is present as the result of any	#
613 # instruction, the 060 will take an overflow exception whether the	#
624 # _real_ovfl() so that the operating system enabled overflow handler	#
626 #	Two other conditions exist. First, if overflow was disabled	#
631 # overflow was disabled and the trace exception was enabled, this	#
722 # the exceptional possibilities we have left ourselves with are ONLY overflow
723 # and inexact. and, the inexact is such that overflow occurred and was disabled
739 # overflow is enabled AND overflow, of course, occurred. so, we have the EXOP
756 # overflow occurred but is disabled. meanwhile, inexact is enabled. Therefore,
864 # _real_unfl() so that the operating system enabled overflow handler	#
970 # the exceptional possibilities we have left ourselves with are ONLY overflow
971 # and inexact. and, the inexact is such that overflow occurred and was disabled
989 # overflow is enabled AND overflow, of course, occurred. so, we have the EXOP
1361 	btst		&ovfl_bit,FPSR_EXCEPT(%a6) # was overflow set?
1367 	bra.w		fu_in_exc_ovfl		# go insert overflow frame
1392 	btst		&ovfl_bit,FPSR_EXCEPT(%a6) # did disabled overflow occur?
1621 # if a disabled overflow occurred and inexact was enabled but the result
1623 	btst		&ovfl_bit,FPSR_EXCEPT(%a6) # was overflow set?
1939 	btst		&ovfl_bit,FPSR_EXCEPT(%a6) # was overflow set?
1970 	btst		&ovfl_bit,FPSR_EXCEPT(%a6) # did disabled overflow occur?
2653 # no exception occurred. now, did a disabled, exact overflow occur with inexact
2654 # enabled? if so, then we have to stuff an overflow frame into the FPU.
2655 	btst		&ovfl_bit,FPSR_EXCEPT(%a6) # did overflow occur?
2670 # the enabled exception was inexact. so, if it occurs with an overflow
2671 # or underflow that was disabled, then we have to force an overflow or
2673 	btst		&ovfl_bit,FPSR_EXCEPT(%a6) # did overflow occur?
4185 # no exception that was enabled was set BUT if we got an exact overflow
4186 # and overflow wasn't enabled but inexact was (yech!) then this is
4188 	btst		&ovfl_bit,FPSR_EXCEPT(%a6) # did overflow occur?
4191 # the overflow w/ exact result happened but was inexact set in the FPCR?
4204 # the enabled exception was inexact. so, if it occurs with an overflow
4205 # or underflow that was disabled, then we have to force an overflow or
4206 # underflow frame. the eventual overflow or underflow handler will see that
4209 # a disabled overflow or underflow has also happened.
4210 	btst		&ovfl_bit,FPSR_EXCEPT(%a6) # did overflow occur?
4229 # insert an overflow frame
5577 #--there is a danger of unwanted overflow in first LOOP iteration.  In this
6008 #--there is a danger of unwanted overflow in first LOOP iteration.  In this
6925 #		neither overflow nor underflow. If AdjFlag = 1, that	#
6928 #		Hence, exp(X) may overflow or underflow or neither.	#
6931 #		over/underflow. Possible exception in 6.4 is overflow	#
6969 #	Notes:	Exp(X) will surely overflow or underflow, depending on	#
6973 #		inexact together with either overflow or underflow.	#
7702 #		overflow. Thus, we calculate sinh(X) as follows:	#
7709 #	5. (|X| > 16480 log2) sinh(X) must overflow. Return		#
7710 #		Huge*Huge to generate overflow and an infinity with	#
7815 #          However, invoking exp(|X|) may cause premature overflow.	#
7824 #       5. (|X| > 16480 log2) sinh(X) must overflow. Return		#
7825 #          sign(X)*Huge*Huge to generate overflow and an infinity with	#
8981 #	1. Generate overflow by Huge * Huge if X > 0; otherwise,	#
9469 # or overflow.
9566 #	    >= 2^14, an overflow or underflow is returned.		#
9679 	bra		t_ovfl_sc		# no; overflow
10343 #	  routine may have caused an underflow or overflow.	#
10381 # an underflow or overflow exception occurred.
11531 #	ovf_res() - return default overflow result			#
11550 # an exception. If so, return the default overflow/underflow result	#
11607 	beq.w		fmul_may_ovfl		# result may rnd to overflow
11608 	blt.w		fmul_ovfl		# result will overflow
11616 # - the result of the multiply operation will neither overflow nor underflow.
11651 # - the result of the multiply operation is an overflow.
11655 # - if overflow or inexact is enabled, we need a multiply result rounded to
11735 # - the result of the multiply operation MAY overflow.
11755 	fbge.w		fmul_ovfl_tst		# yes; overflow has occurred
11757 # no, it didn't overflow; we have correct result
11766 # - if overflow or inexact is enabled, we need a multiply result rounded to
12026 #	ovf_res() - return default overflow result			#
12078 # or overflow because of rounding to the correct precision. so...
12145 	cmpi.l		%d0,&0x3fff-0x407e	# will move in overflow?
12147 	blt.w		fin_sd_ovfl		# yes; go handle overflow
12150 # operand will NOT overflow or underflow when moved into the fp reg file
12188 	cmpi.l		%d0,&0x3fff-0x43fe	# will move in overflow?
12190 	blt.w		fin_sd_ovfl		# yes; go handle overflow
12240 # operand WILL overflow.
12294 # the move in MAY overflow. so...
12309 	fbge.w		fin_sd_ovfl_tst		# yes; overflow has occurred
12311 # no, it didn't overflow; we have correct result
12348 #	ovf_res() - return default overflow result			#
12367 # an exception. If so, return the default overflow/underflow result	#
12380 	long		0x3fff - 0x7ffe		# ext overflow exponent
12381 	long		0x3fff - 0x407e		# sgl overflow exponent
12382 	long		0x3fff - 0x43fe		# dbl overflow exponent
12429 	cmp.l		%d0,(tbl_fdiv_ovfl.b,%pc,%d1.w*4) # will result overflow?
12430 	ble.w		fdiv_may_ovfl		# result will overflow
12801 #	ovf_res() - return default overflow result			#
12818 # and an actual fneg performed to see if overflow/underflow would have	#
12819 # occurred. If so, return default underflow/overflow result. Else,	#
12851 # or overflow because of rounding to the correct precision. so...
12925 	cmpi.l		%d0,&0x3fff-0x407e	# will move in overflow?
12927 	blt.w		fneg_sd_ovfl		# yes; go handle overflow
12930 # operand will NOT overflow or underflow when moved in to the fp reg file
12968 	cmpi.l		%d0,&0x3fff-0x43fe	# will move in overflow?
12970 	blt.w		fneg_sd_ovfl		# yes; go handle overflow
13020 # operand WILL overflow.
13089 	fbge.w		fneg_sd_ovfl_tst	# yes; overflow has occurred
13091 # no, it didn't overflow; we have correct result
13421 #	ovf_res() - calculate overflow result				#
13437 #	Double and single precision can overflow and underflow. First,	#
13474 # or overflow because of rounding to the correct precision. so...
13543 	cmpi.l		%d0,&0x3fff-0x407e	# will move in overflow?
13545 	blt.w		fabs_sd_ovfl		# yes; go handle overflow
13548 # operand will NOT overflow or underflow when moved in to the fp reg file
13586 	cmpi.l		%d0,&0x3fff-0x43fe	# will move in overflow?
13588 	blt.w		fabs_sd_ovfl		# yes; go handle overflow
13635 # operand WILL overflow.
13704 	fbge.w		fabs_sd_ovfl_tst	# yes; overflow has occurred
13706 # no, it didn't overflow; we have correct result
13933 #	ovf_res() - return default overflow result			#
13952 # an exception. If so, return the default overflow/underflow result	#
13986 	beq.w		fsglmul_may_ovfl	# result may rnd to overflow
13987 	blt.w		fsglmul_ovfl		# result will overflow
14084 	fbge.w		fsglmul_ovfl_tst	# yes; overflow has occurred
14086 # no, it didn't overflow; we have correct result
14274 #	ovf_res() - return default overflow result			#
14293 # an exception. If so, return the default overflow/underflow result	#
14385 	cmp.l		%d1,&0x7fff		# did divide overflow?
14609 #	ovf_res() - return default overflow result			#
14689 	cmp.l		%d2,(tbl_fadd_ovfl.b,%pc,%d1.w*4) # is it an overflow?
15062 #	ovf_res() - return default overflow result			#
15142 	cmp.l		%d2,(tbl_fsub_ovfl.b,%pc,%d1.w*4) # is it an overflow?
15509 #	ovf_res() - return default overflow result			#
15527 # an exception. If so, return the default overflow/underflow result	#
15606 	cmpi.l		%d0,&0x3fff-0x407f	# will move in overflow?
15608 	blt.w		fsqrt_sd_ovfl		# yes; go handle overflow
15611 # operand will NOT overflow or underflow when moved in to the fp reg file
15651 	cmpi.l		%d0,&0x3fff-0x43ff	# will move in overflow?
15653 	blt.w		fsqrt_sd_ovfl		# yes; go handle overflow
15717 # operand WILL overflow.
15775 	bne.w		fsqrt_sd_ovfl		# yes, so overflow
15789 	fbge.w		fsqrt_sd_ovfl_tst	# yes; overflow has occurred
15791 # no, it didn't overflow; we have correct result
20155 #	ovf_res() - create default overflow result for sgl/dbl precision#
20172 #	fp0 : intermediate underflow or overflow result if		#
20183 #	For sgl or dbl precision, overflow or underflow can occur. If	#
20440 # would cause either an underflow or overflow. these cases are handled
20447 	cmpi.w		%d0,&SGL_HI		# will operand overflow?
20543 # it's definitely an overflow so call ovf_res to get the correct answer
20559 # overflow result. DON'T save the returned ccodes from ovf_res() since
20565 	fmovm.x		(%a0),&0x80		# load default overflow result
20593 # move out MAY overflow:
20616 	bra.w		fout_sgl_ovfl		# yes; go handle overflow
20678 # would cause either an underflow or overflow. these cases are handled
20685 	cmpi.w		%d0,&DBL_HI		# will operand overflow?
20764 # it's definitely an overflow so call ovf_res to get the correct answer
20780 # overflow result. DON'T save the returned ccodes from ovf_res() since
20786 	fmovm.x		(%a0),&0x80		# load default overflow result
20804 # move out MAY overflow:
20827 	bra.w		fout_dbl_ovfl		# yes; go handle overflow
22131 	bcc.b		scc_clr			# no mantissa overflow
22824 #	ovf_res(): routine to produce the default overflow result of	#
22845 #	The default overflow result can be determined by the sign of	#
23012 #	Note: this operation can never overflow.			#
23018 #	Note: this operation can never overflow.			#
23900 #     to overflow.  Only a negative iscale can cause this, so