ifpsp060/src/pfpsp.S

414 set SGL_LO,		0x3f81			# min sgl prec exponent
415 set SGL_HI,		0x407e			# max sgl prec exponent
515 set sgl_thresh,		0x3f81			# minimum sgl exponent
1399 	bsr.l		funimp_skew		# skew sgl or dbl inputs
1432 	cmpi.b		%d0,&0x11		# is class = 2 & fmt = sgl?
1976 # "non-skewed" operand for cases of sgl and dbl src INFs,NANs, and DENORMs.
2303 	cmpi.b		%d0,&0x1		# was src sgl?
3132 # denorm operand in the sgl or dbl format. NANs also become skewed, but can't
3196 	short		tbl_operr    - tbl_operr # sgl prec shouldn't happen
3328 # denorm operand in the sgl or dbl format. NANs also become skewed and must be
3366 	short		fsnan_out_s - tbl_snan # sgl prec shouldn't happen
3432 	lsr.l		&0x8,%d1		# shift mantissa for sgl
3433 	or.l		%d1,%d0			# create sgl SNAN
3447 	lsr.l		&0x8,%d1		# shift mantissa for sgl
3448 	or.l		%d1,%d0			# create sgl SNAN
3791 # in the sgl or dbl format.
6142 #			sgl = $0004xxxx					#
6205 	cmpi.b		%d1, &s_mode		# is prec = sgl?
6222 	cmpi.b		%d1, &s_mode		# is prec = sgl?
6239 	cmpi.b		%d1, &s_mode		# is prec = sgl?
6245 set	ad_1_sgl,	0x00000100	# constant to add 1 to l-bit in sgl prec
6262 	and.l		&0xffffff00, FTEMP_HI(%a0) # truncate bits beyond sgl limit
6315 	cmpi.b		%d1, &s_mode		# is prec sgl?
6345 	bne.b		ext_grs_not_ext		# no; go handle sgl or dbl
6358 	cmpi.b		%d1, &s_mode		# is rnd prec = sgl?
6362 # sgl:
6374 	bfextu		FTEMP_HI(%a0){&24:&2}, %d3 # sgl prec. g-r are 2 bits right
6375 	mov.l		&30, %d2		# of the sgl prec. limits
6730 #	set_tag_s(): return the optype of the input sgl fp number	#
6903 	mov.w		&s_mode,%d1		# force rnd prec = sgl
7054 #	_round() - needed to create EXOP for sgl/dbl precision		#
7056 #	ovf_res() - create default overflow result for sgl/dbl precision#
7057 #	unf_res() - create default underflow result for sgl/dbl prec.	#
7059 #	dst_sgl() - create rounded sgl precision result.		#
7074 #	      OVFL/UNFL occurred for a sgl or dbl operand		#
7084 #	For sgl or dbl precision, overflow or underflow can occur. If	#
7159 	ori.l		&0x00800000,%d1		# make smallest sgl
7205 	ori.l		&0x00800000,%d1		# make smallest sgl
7252 	ori.l		&0x00800000,%d1		# make smallest sgl
7336 	ori.b		&s_mode*0x10,%d0	# insert sgl prec
7459 # call ovf_res() w/ sgl prec and the correct rnd mode to create the default
7464 	mov.l		L_SCR3(%a6),%d0		# pass: sgl prec,rnd mode
7814 #	sgl_exp = ext_exp - $3fff(ext bias) + $7f(sgl bias)		#
7977 # norms/denorms into ext/sgl/dbl precision.				#
8000 	ori.b		&s_mode*0x10,%d0	# insert sgl prec
8133 	bne.b		fmul_ovfl_ena_sd	# it's sgl or dbl
8244 	bne.b		fmul_unfl_ena_sd	# no, sgl or dbl
8474 # sgl/dbl, must scale exponent and perform an "fmove". Check to see	#
8485 	ori.b		&s_mode*0x10,%d0	# insert sgl precision
8505 	bne.w		fin_not_ext		# no, so go handle dbl or sgl
8525 	bne.w		fin_not_ext		# no, so go handle dbl or sgl
8562 	cmpi.b		%d0,&s_mode*0x10	# separate sgl/dbl prec
8794 # norms/denorms into ext/sgl/dbl precision.				#
8812 	long		0x3fff - 0x407e		# sgl overflow exponent
8818 	ori.b		&s_mode*0x10,%d0	# insert sgl prec
8946 	bne.b		fdiv_ovfl_ena_sd	# no, do sgl or dbl
9014 	bne.b		fdiv_unfl_ena_sd	# no, sgl or dbl
9231 #	scale_to_zero_src() - scale sgl/dbl source exponent		#
9247 # norms/denorms into ext/sgl/dbl precisions. Extended precision can be	#
9259 	ori.b		&s_mode*0x10,%d0	# insert sgl precision
9278 	bne.w		fneg_not_ext		# no; go handle sgl or dbl
9302 	bne.b		fneg_not_ext		# no; go handle sgl or dbl
9342 	cmpi.b		%d0,&s_mode*0x10	# separate sgl/dbl prec
9882 	ori.b		&s_mode*0x10,%d0	# insert sgl precision
9901 	bne.b		fabs_not_ext		# no; go handle sgl or dbl
9960 	cmpi.b		%d0,&s_mode*0x10	# separate sgl/dbl prec
10379 # norms/denorms into ext/sgl/dbl precision.				#
10430 	fsglmul.x	FP_SCR0(%a6),%fp0	# execute sgl multiply
10457 	fsglmul.x	FP_SCR0(%a6),%fp0	# execute sgl multiply
10506 	fsglmul.x	FP_SCR0(%a6),%fp0	# execute sgl multiply
10528 	fsglmul.x	FP_SCR0(%a6),%fp0	# execute sgl multiply
10558 	fsglmul.x	FP_SCR0(%a6),%fp1	# execute sgl multiply
10583 	fsglmul.x	FP_SCR0(%a6),%fp0	# execute sgl multiply
10611 	fsglmul.x	FP_SCR0(%a6),%fp1	# execute sgl multiply
10720 # norms/denorms into ext/sgl/dbl precision.				#
10777 	fsgldiv.x	FP_SCR0(%a6),%fp0	# perform sgl divide
10861 	fsgldiv.x	FP_SCR0(%a6),%fp0	# execute sgl divide
10891 	fsgldiv.x	FP_SCR0(%a6),%fp1	# execute sgl divide
10919 	fsgldiv.x	FP_SCR0(%a6),%fp0	# execute sgl divide
10946 	fsgldiv.x	FP_SCR0(%a6),%fp1	# execute sgl divide
11069 	ori.b		&s_mode*0x10,%d0	# insert sgl prec
11144 	long		0x407f			# sgl ovfl
11149 	long		0x3f81			# sgl unfl
11173 	bne.b		fadd_ovfl_ena_sd	# no; prec = sgl or dbl
11238 	bne.b		fadd_unfl_ena_sd	# no; sgl or dbl
11522 	ori.b		&s_mode*0x10,%d0	# insert sgl prec
11597 	long		0x407f			# sgl ovfl
11602 	long		0x3f81			# sgl unfl
11954 # norms/denorms into ext/sgl/dbl precision.				#
11967 	ori.b		&s_mode*0x10,%d0	# insert sgl precision
11990 	bne.b		fsqrt_not_ext		# no; go handle sgl or dbl
12007 	bne.b		fsqrt_not_ext		# no; go handle sgl or dbl
12021 	cmpi.b		%d0,&s_mode*0x10	# separate sgl/dbl prec