ifpsp060/src/fpsp.S

415 set SGL_LO,		0x3f81			# min sgl prec exponent
416 set SGL_HI,		0x407e			# max sgl prec exponent
516 set sgl_thresh,		0x3f81			# minimum sgl exponent
1400 	bsr.l		funimp_skew		# skew sgl or dbl inputs
1433 	cmpi.b		%d0,&0x11		# is class = 2 & fmt = sgl?
1977 # "non-skewed" operand for cases of sgl and dbl src INFs,NANs, and DENORMs.
2304 	cmpi.b		%d0,&0x1		# was src sgl?
3133 # denorm operand in the sgl or dbl format. NANs also become skewed, but can't
3197 	short		tbl_operr    - tbl_operr # sgl prec shouldn't happen
3329 # denorm operand in the sgl or dbl format. NANs also become skewed and must be
3367 	short		fsnan_out_s - tbl_snan # sgl prec shouldn't happen
3433 	lsr.l		&0x8,%d1		# shift mantissa for sgl
3434 	or.l		%d1,%d0			# create sgl SNAN
3448 	lsr.l		&0x8,%d1		# shift mantissa for sgl
3449 	or.l		%d1,%d0			# create sgl SNAN
3792 # in the sgl or dbl format.
5309 # here, the operation may underflow iff the precision is sgl or dbl.
10298 	bne.b		ovfl_sc_dbl		# no; sgl
10300 	tst.l		LOCAL_LO(%a0)		# is lo lw of sgl set?
10446 	long		0x3f810000, 0x00000000, 0x00000000, 0x0 # ZERO;sgl
10447 	long		0x3f810000, 0x00000000, 0x00000000, 0x0 # ZERO;sgl
10448 	long		0x3f810000, 0x00000000, 0x00000000, 0x0 # ZERO;sgl
10449 	long		0x3f810000, 0x00000100, 0x00000000, 0x0 # MIN; sgl
10466 	long		0xbf810000, 0x00000000, 0x00000000, 0x0 # ZERO;sgl
10467 	long		0xbf810000, 0x00000000, 0x00000000, 0x0 # ZERO;sgl
10468 	long		0xbf810000, 0x00000100, 0x00000000, 0x0 # MIN; sgl
10469 	long		0xbf810000, 0x00000000, 0x00000000, 0x0 # ZERO;sgl
11546 # norms/denorms into ext/sgl/dbl precision.				#
11569 	ori.b		&s_mode*0x10,%d0	# insert sgl prec
11702 	bne.b		fmul_ovfl_ena_sd	# it's sgl or dbl
11813 	bne.b		fmul_unfl_ena_sd	# no, sgl or dbl
12043 # sgl/dbl, must scale exponent and perform an "fmove". Check to see	#
12054 	ori.b		&s_mode*0x10,%d0	# insert sgl precision
12074 	bne.w		fin_not_ext		# no, so go handle dbl or sgl
12094 	bne.w		fin_not_ext		# no, so go handle dbl or sgl
12131 	cmpi.b		%d0,&s_mode*0x10	# separate sgl/dbl prec
12363 # norms/denorms into ext/sgl/dbl precision.				#
12381 	long		0x3fff - 0x407e		# sgl overflow exponent
12387 	ori.b		&s_mode*0x10,%d0	# insert sgl prec
12515 	bne.b		fdiv_ovfl_ena_sd	# no, do sgl or dbl
12583 	bne.b		fdiv_unfl_ena_sd	# no, sgl or dbl
12800 #	scale_to_zero_src() - scale sgl/dbl source exponent		#
12816 # norms/denorms into ext/sgl/dbl precisions. Extended precision can be	#
12828 	ori.b		&s_mode*0x10,%d0	# insert sgl precision
12847 	bne.w		fneg_not_ext		# no; go handle sgl or dbl
12871 	bne.b		fneg_not_ext		# no; go handle sgl or dbl
12911 	cmpi.b		%d0,&s_mode*0x10	# separate sgl/dbl prec
13451 	ori.b		&s_mode*0x10,%d0	# insert sgl precision
13470 	bne.b		fabs_not_ext		# no; go handle sgl or dbl
13529 	cmpi.b		%d0,&s_mode*0x10	# separate sgl/dbl prec
13948 # norms/denorms into ext/sgl/dbl precision.				#
13999 	fsglmul.x	FP_SCR0(%a6),%fp0	# execute sgl multiply
14026 	fsglmul.x	FP_SCR0(%a6),%fp0	# execute sgl multiply
14075 	fsglmul.x	FP_SCR0(%a6),%fp0	# execute sgl multiply
14097 	fsglmul.x	FP_SCR0(%a6),%fp0	# execute sgl multiply
14127 	fsglmul.x	FP_SCR0(%a6),%fp1	# execute sgl multiply
14152 	fsglmul.x	FP_SCR0(%a6),%fp0	# execute sgl multiply
14180 	fsglmul.x	FP_SCR0(%a6),%fp1	# execute sgl multiply
14289 # norms/denorms into ext/sgl/dbl precision.				#
14346 	fsgldiv.x	FP_SCR0(%a6),%fp0	# perform sgl divide
14430 	fsgldiv.x	FP_SCR0(%a6),%fp0	# execute sgl divide
14460 	fsgldiv.x	FP_SCR0(%a6),%fp1	# execute sgl divide
14488 	fsgldiv.x	FP_SCR0(%a6),%fp0	# execute sgl divide
14515 	fsgldiv.x	FP_SCR0(%a6),%fp1	# execute sgl divide
14638 	ori.b		&s_mode*0x10,%d0	# insert sgl prec
14713 	long		0x407f			# sgl ovfl
14718 	long		0x3f81			# sgl unfl
14742 	bne.b		fadd_ovfl_ena_sd	# no; prec = sgl or dbl
14807 	bne.b		fadd_unfl_ena_sd	# no; sgl or dbl
15091 	ori.b		&s_mode*0x10,%d0	# insert sgl prec
15166 	long		0x407f			# sgl ovfl
15171 	long		0x3f81			# sgl unfl
15523 # norms/denorms into ext/sgl/dbl precision.				#
15536 	ori.b		&s_mode*0x10,%d0	# insert sgl precision
15559 	bne.b		fsqrt_not_ext		# no; go handle sgl or dbl
15576 	bne.b		fsqrt_not_ext		# no; go handle sgl or dbl
15590 	cmpi.b		%d0,&s_mode*0x10	# separate sgl/dbl prec
19562 #	set_tag_s() - determine sgl prec optype tag			#
19766 	bsr.l		fetch_dreg		# fetch sgl in d0
19769 	lea		L_SCR1(%a6), %a0	# pass: ptr to the sgl
19770 	bsr.l		set_tag_s		# determine sgl type
19943 	lea		L_SCR1(%a6), %a0	# pass: ptr to sgl src op
19947 	cmpi.b		%d0, &DENORM		# is it a sgl DENORM?
19950 	cmpi.b		%d0, &SNAN		# is it a sgl SNAN?
19964 # must convert sgl denorm format to an Xprec denorm fmt suitable for
19966 # %a0 : points to sgl denorm
19969 	bfextu		(%a0){&9:&23}, %d0	# fetch sgl hi(_mantissa)
19989 # convert sgl to ext SNAN
19990 # %a0 : points to sgl SNAN
20153 #	_round() - needed to create EXOP for sgl/dbl precision		#
20155 #	ovf_res() - create default overflow result for sgl/dbl precision#
20156 #	unf_res() - create default underflow result for sgl/dbl prec.	#
20158 #	dst_sgl() - create rounded sgl precision result.		#
20173 #	      OVFL/UNFL occurred for a sgl or dbl operand		#
20183 #	For sgl or dbl precision, overflow or underflow can occur. If	#
20258 	ori.l		&0x00800000,%d1		# make smallest sgl
20304 	ori.l		&0x00800000,%d1		# make smallest sgl
20351 	ori.l		&0x00800000,%d1		# make smallest sgl
20435 	ori.b		&s_mode*0x10,%d0	# insert sgl prec
20558 # call ovf_res() w/ sgl prec and the correct rnd mode to create the default
20563 	mov.l		L_SCR3(%a6),%d0		# pass: sgl prec,rnd mode
20913 #	sgl_exp = ext_exp - $3fff(ext bias) + $7f(sgl bias)		#
22020 #			sgl = $0004xxxx					#
22083 	cmpi.b		%d1, &s_mode		# is prec = sgl?
22100 	cmpi.b		%d1, &s_mode		# is prec = sgl?
22117 	cmpi.b		%d1, &s_mode		# is prec = sgl?
22123 set	ad_1_sgl,	0x00000100	# constant to add 1 to l-bit in sgl prec
22140 	and.l		&0xffffff00, FTEMP_HI(%a0) # truncate bits beyond sgl limit
22193 	cmpi.b		%d1, &s_mode		# is prec sgl?
22223 	bne.b		ext_grs_not_ext		# no; go handle sgl or dbl
22236 	cmpi.b		%d1, &s_mode		# is rnd prec = sgl?
22240 # sgl:
22252 	bfextu		FTEMP_HI(%a0){&24:&2}, %d3 # sgl prec. g-r are 2 bits right
22253 	mov.l		&30, %d2		# of the sgl prec. limits
22608 #	set_tag_s(): return the optype of the input sgl fp number	#
22781 	mov.w		&s_mode,%d1		# force rnd prec = sgl