ifpsp060/src/fplsp.S

4933 #		even polynomial in r, 1 + r*r*(B1+s*(B2+ ... + s*B8)),	#
5129 #--SGN + S'*(B1 + S(B2 + S(B3 + S(B4 + ... + SB8)))), WHERE
5131 #--SGN + S'*([B1+T(B3+T(B5+TB7))] + [S(B2+T(B4+T(B6+TB8)))])
5134 #--WHILE B2 AND B3 ARE IN DOUBLE-EXTENDED FORMAT, B1 IS -1/2
5174 	fadd.x		COSB2(%pc),%fp2		# B2+T(B4+T(B6+TB8))
5177 	fmul.x		%fp2,%fp0		# S(B2+T(B4+T(B6+TB8)))
5339 	fadd.x		COSB2(%pc),%fp2		# B2+S(B3+...)
5342 	fmul.x		%fp2,%fp0		# S(B2+...)
5345 	fadd.s		COSB1(%pc),%fp0		# B1+S(B2...)
5346 	fmul.x		SPRIME(%a6),%fp0	# S'(B1+S(B2+...))
5410 	fadd.x		COSB2(%pc),%fp1		# B2+S(B3+...)
5413 	fmul.x		%fp0,%fp1		# S(B2+...)
5417 	fadd.s		COSB1(%pc),%fp1		# B1+S(B2...)
5419 	fmul.x		SPRIME(%a6),%fp1	# S'(B1+S(B2+...))
6352 #--ATAN(X) BY X + X*Y*(B1+Y*(B2+Y*(B3+Y*(B4+Y*(B5+Y*B6)))))
6353 #--WHICH IS X + X*Y*( [B1+Z*(B3+Z*B5)] + [Y*(B2+Z*(B4+Z*B6)] )
6379 	fadd.d		ATANB2(%pc),%fp2	# B2+Z*(B4+Z*B6)
6382 	fmul.x		%fp0,%fp2		# Y*(B2+Z*(B4+Z*B6))
6385 	fadd.x		%fp2,%fp1		# [B1+Z*(B3+Z*B5)]+[Y*(B2+Z*(B4+Z*B6))]
6387 	fmul.x		%fp1,%fp0		# X*Y*([B1+Z*(B3+Z*B5)]+[Y*(B2+Z*(B4+Z*B6))])
6447 #					...	+[Y*(B2+Z*(B4+Z*B6))])
6971 #			p = X + X*X*(B1 + X*(B2 + ... + X*B12))		#
6975 #		precision; and B2 is double extended.			#
6983 #			Q	=	X*S*(B2 + X*(B3 + ... + X*B12))	#
6986 #			Q = [ X*S*(B2 + S*(B4 + ... + S*B12)) ] +	#
7459 	fadd.x		EM1B2(%pc),%fp1		# fp1 is B2+S*...
7462 	fmul.x		%fp0,%fp1		# fp1 is S*(B2+...
7465 	fmul.x		(%a0),%fp1		# fp1 is X*S*(B2...
8341 #--U + U*V*(B1 + V*(B2 + V*(B3 + V*(B4 + V*B5)))) BY
8342 #--U + U*V*(  [B1 + W*(B3 + W*B5)]  +  [V*(B2 + W*B4)]  )
8356 	fadd.d		LOGB2(%pc),%fp2		# B2+W*B4
8360 	fmul.x		%fp0,%fp2		# V*(B2+W*B4)
8365 	fadd.x		%fp2,%fp1		# B1+W*(B3+W*B5) + V*(B2+W*B4), FP2 RELEASED
8368 	fmul.x		%fp1,%fp0		# U*V*( [B1+W*(B3+W*B5)] + [V*(B2+W*B4)] )