[unix-history] / usr / src / sys / tahoe / math / Ksqrtf.s

/*	Ksqrtf.s	1.4	90/12/04	*/

#include "../tahoe/SYS.h"

	.text
ENTRY(Ksqrtf, 0)
	subl3	$76,fp,sp
	tstl	4(fp)			# if (arg <= 0.) {
	jgtr	L18
	jgeq	L19			# if (arg < 0.)
L19:	clrl	r1
	clrl	r0
	ret				# return(0.);}
L18:	
	pushl	20(fp)			# hfs
	subl3	$72,fp,-(sp)
	pushl	8(fp)
	pushl	4(fp)
	callf	$20,_Kfrexpf
	clrl	-60+4(fp)
	movl	r0,-60(fp)		# x = frexp(arg,&exp);
	jbr	L20
L2000001:
	pushl	20(fp)			# hfs
	ldd	two
	pushd
	ldd	-60(fp)
	pushd
	callf	$24,_Kmuld		# x * 2
	ldd	r0
	std	-60(fp)			# x *= 2;
	decl	-72(fp)			# exp--;}
L20:	cmpd2	-60(fp),half		# while(x < 0.5){
	jlss	L2000001
	bitl	$1,-72(fp)		# if(exp & 1) {
	jeql	L24
	pushl	20(fp)			# hfs
	ldd	two
	pushd
	ldd	-60(fp)
	pushd
	callf	$24,_Kmuld		# x * 2
	ldd	r0
	std	-60(fp)			# x *= 2;
	decl	-72(fp)			# exp--;}
L24:	pushl	20(fp)			# hfs
	ldd	-60(fp)
	pushd
	ldd	one
	pushd
	callf	$24,_Kaddd		# (1.0+x)
	pushl	20(fp)			# hfs
	ldd	r0
	pushd
	ldd	half
	pushd
	callf	$24,_Kmuld
	ldd	r0
	std	-68(fp)			 # temp=0.5*(1.0+x);
	jbr	L28
L2000003:
	pushl	20(fp)			# hfs
	ldd	big
	pushd
	ldd	-68(fp)
	pushd
	callf	$24,_Kmuld		# temp*(1L<<30)
	ldd	r0
	std	-68(fp)			# temp *=(1L<<30);
	subl2	$60,-72(fp)		# exp -= 60;}
L28:	cmpl	-72(fp),$60		# while(exp > 60){
	jleq	L31
	jbr	L2000003
L2000005:
	pushl	20(fp)			# hfs
	ldd	big
	pushd
	ldd	-68(fp)
	pushd
	callf	$24,_Kdivd
	ldd	r0
	std	-68(fp)			# temp /= (1L<<30);
	addl2	$60,-72(fp)		# exp +=60;}
L31:	cmpl	-72(fp),$-60		# while(exp < -60) {
	jlss	L2000005
	tstl	-72(fp)			# if(exp >=0)
	jlss	L34
	divl3	$2,-72(fp),r0
	shll	r0,$1,r0
	pushl	20(fp)			# hfs
	pushl	$0
	pushl	r0
	pushl	$0
	pushl	$0
	callf	$24,_Kcvtld
	pushl	20(fp)			# hfs
	ldd	r0
	pushd
	ldd	-68(fp)
	pushd
	callf	$24,_Kmuld		# temp * (1L <<(exp/2))
	jbr	L2000006		# temp *= 1L << (exp/2));
L34:	mnegl	-72(fp),r0		# -exp
	divl2	$2,r0			# -exp/2
	shll	r0,$1,r0
	pushl	20(fp)			# hfs
	pushl	$0
	pushl	r0
	pushl	$0
	pushl	$0
	callf	$24,_Kcvtld
	pushl	20(fp)			# hfs
	ldd	r0
	pushd
	ldd	-68(fp)
	pushd
	callf	$24,_Kdivd		# temp / (1L <<(-exp/2))
L2000006:
	ldd	r0
	std	-68(fp)			# temp (*/)= 1L << (exp/2));
	clrl	-76(fp)
L2000008:
	pushl	20(fp)			# hfs
	ldd	-68(fp)
	pushd
	ldd	4(fp)
	pushd
	callf	$24,_Kdivd
	pushl	20(fp)			# hfs
	ldd	r0
	pushd
	ldd	-68(fp)
	pushd
	callf	$24,_Kaddd
	pushl	20(fp)			# hfs
	ldd	r0	
	pushd
	ldd	half
	pushd
	callf	$24,_Kmuld
	ldd	r0
	std	-68(fp)			# temp=0.5*(temp+arg/temp);
	aobleq	$3,-76(fp),L2000008
	ldd	-68(fp)
	cvdf
	stf	r0
	ret

	.data
	.align	2
half:	.long	0x40000000, 0x00000000 # .double .5
two:	.long	0x41000000, 0x00000000 # .double 2
one:	.long	0x40800000, 0x00000000 # .double 1
big:	.long	0x4F800000, 0x00000000 # .double (1<<30)
Commit	Line	Data
a6730a00	1	/* Ksqrtf.s 1.4 90/12/04 */
2b252132	2
9d1834ba	3	#include "../tahoe/SYS.h"
2b252132	4
6edc04d6	5	.text
2b252132	6	ENTRY(Ksqrtf, 0)
6edc04d6 SL	7	subl3 $76,fp,sp
	8	tstl 4(fp) # if (arg <= 0.) {
	9	jgtr L18
	10	jgeq L19 # if (arg < 0.)
6edc04d6 SL	11	L19: clrl r1
	12	clrl r0
	13	ret # return(0.);}
	14	L18:
	15	pushl 20(fp) # hfs
	16	subl3 $72,fp,-(sp)
	17	pushl 8(fp)
	18	pushl 4(fp)
	19	callf $20,_Kfrexpf
	20	clrl -60+4(fp)
	21	movl r0,-60(fp) # x = frexp(arg,&exp);
	22	jbr L20
	23	L2000001:
	24	pushl 20(fp) # hfs
2b252132	25	ldd two
6edc04d6 SL	26	pushd
	27	ldd -60(fp)
	28	pushd
	29	callf $24,_Kmuld # x * 2
	30	ldd r0
	31	std -60(fp) # x *= 2;
	32	decl -72(fp) # exp--;}
2b252132	33	L20: cmpd2 -60(fp),half # while(x < 0.5){
6edc04d6 SL	34	jlss L2000001
	35	bitl $1,-72(fp) # if(exp & 1) {
	36	jeql L24
	37	pushl 20(fp) # hfs
2b252132	38	ldd two
6edc04d6 SL	39	pushd
	40	ldd -60(fp)
	41	pushd
	42	callf $24,_Kmuld # x * 2
	43	ldd r0
	44	std -60(fp) # x *= 2;
	45	decl -72(fp) # exp--;}
	46	L24: pushl 20(fp) # hfs
	47	ldd -60(fp)
	48	pushd
2b252132	49	ldd one
6edc04d6 SL	50	pushd
	51	callf $24,_Kaddd # (1.0+x)
	52	pushl 20(fp) # hfs
	53	ldd r0
	54	pushd
2b252132	55	ldd half
6edc04d6 SL	56	pushd
	57	callf $24,_Kmuld
	58	ldd r0
	59	std -68(fp) # temp=0.5*(1.0+x);
	60	jbr L28
	61	L2000003:
	62	pushl 20(fp) # hfs
2b252132	63	ldd big
6edc04d6 SL	64	pushd
	65	ldd -68(fp)
	66	pushd
	67	callf $24,_Kmuld # temp*(1L<<30)
	68	ldd r0
	69	std -68(fp) # temp *=(1L<<30);
	70	subl2 $60,-72(fp) # exp -= 60;}
	71	L28: cmpl -72(fp),$60 # while(exp > 60){
	72	jleq L31
	73	jbr L2000003
	74	L2000005:
	75	pushl 20(fp) # hfs
2b252132	76	ldd big
6edc04d6 SL	77	pushd
	78	ldd -68(fp)
	79	pushd
	80	callf $24,_Kdivd
	81	ldd r0
	82	std -68(fp) # temp /= (1L<<30);
	83	addl2 $60,-72(fp) # exp +=60;}
	84	L31: cmpl -72(fp),$-60 # while(exp < -60) {
	85	jlss L2000005
	86	tstl -72(fp) # if(exp >=0)
	87	jlss L34
	88	divl3 $2,-72(fp),r0
	89	shll r0,$1,r0
	90	pushl 20(fp) # hfs
	91	pushl $0
	92	pushl r0
	93	pushl $0
	94	pushl $0
	95	callf $24,_Kcvtld
	96	pushl 20(fp) # hfs
	97	ldd r0
	98	pushd
	99	ldd -68(fp)
	100	pushd
	101	callf $24,_Kmuld # temp * (1L <<(exp/2))
	102	jbr L2000006 # temp *= 1L << (exp/2));
	103	L34: mnegl -72(fp),r0 # -exp
	104	divl2 $2,r0 # -exp/2
	105	shll r0,$1,r0
	106	pushl 20(fp) # hfs
	107	pushl $0
	108	pushl r0
	109	pushl $0
	110	pushl $0
	111	callf $24,_Kcvtld
	112	pushl 20(fp) # hfs
	113	ldd r0
	114	pushd
	115	ldd -68(fp)
	116	pushd
	117	callf $24,_Kdivd # temp / (1L <<(-exp/2))
	118	L2000006:
	119	ldd r0
	120	std -68(fp) # temp (*/)= 1L << (exp/2));
	121	clrl -76(fp)
	122	L2000008:
	123	pushl 20(fp) # hfs
	124	ldd -68(fp)
	125	pushd
	126	ldd 4(fp)
	127	pushd
	128	callf $24,_Kdivd
	129	pushl 20(fp) # hfs
	130	ldd r0
	131	pushd
	132	ldd -68(fp)
	133	pushd
	134	callf $24,_Kaddd
	135	pushl 20(fp) # hfs
	136	ldd r0
	137	pushd
2b252132	138	ldd half
6edc04d6 SL	139	pushd
	140	callf $24,_Kmuld
	141	ldd r0
	142	std -68(fp) # temp=0.5*(temp+arg/temp);
	143	aobleq $3,-76(fp),L2000008
	144	ldd -68(fp)
	145	cvdf
	146	stf r0
	147	ret
2b252132 SL	148
	149	.data
	150	.align 2
	151	half: .long 0x40000000, 0x00000000 # .double .5
	152	two: .long 0x41000000, 0x00000000 # .double 2
	153	one: .long 0x40800000, 0x00000000 # .double 1
	154	big: .long 0x4F800000, 0x00000000 # .double (1<<30)