about summary refs log tree commit diff
path: root/sysdeps/hppa/hppa1.1/udiv_qrnnd.S
blob: e6fb46ae6ebc8b02db7d2bc83aa2624a6bbc798a (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
;! HP-PA  __udiv_qrnnd division support, used from longlong.h.
;! This version runs fast on PA 7000 and later.

;! Copyright (C) 1993-2023 Free Software Foundation, Inc.

;! This file is part of the GNU MP Library.

;! The GNU MP Library is free software; you can redistribute it and/or modify
;! it under the terms of the GNU Lesser General Public License as published by
;! the Free Software Foundation; either version 2.1 of the License, or (at your
;! option) any later version.

;! The GNU MP Library is distributed in the hope that it will be useful, but
;! WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
;! or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
;! License for more details.

;! You should have received a copy of the GNU Lesser General Public License
;! along with the GNU MP Library.  If not, see
;! <https://www.gnu.org/licenses/>.


;! INPUT PARAMETERS
;! rem_ptr	gr26
;! n1		gr25
;! n0		gr24
;! d		gr23

	.text
L$0000:
	.word		0x43f00000
	.word		0x0
	.export		__udiv_qrnnd
__udiv_qrnnd:
	.proc
	.callinfo	frame=64,no_calls
	.entry
	ldo		64(%r30),%r30

	stws		%r25,-16(%r30)	;! n_hi
	stws		%r24,-12(%r30)	;! n_lo
	b,l		L$0,%r1
	ldo		L$0000-L$0(%r1),%r1
L$0:
	fldds		-16(%r30),%fr5
	stws		%r23,-12(%r30)
	comib,<=	0,%r25,L$1
	fcnvxf,dbl,dbl	%fr5,%fr5
	fldds		0(%r1),%fr4
	fadd,dbl	%fr4,%fr5,%fr5
L$1:
	fcpy,sgl	%fr0,%fr6L
	fldws		-12(%r30),%fr6R
	fcnvxf,dbl,dbl	%fr6,%fr4

	fdiv,dbl	%fr5,%fr4,%fr5

	fcnvfx,dbl,dbl	%fr5,%fr4
	fstws		%fr4R,-16(%r30)
	xmpyu		%fr4R,%fr6R,%fr6
	ldws		-16(%r30),%r28
	fstds		%fr6,-16(%r30)
	ldws		-12(%r30),%r21
	ldws		-16(%r30),%r20
	sub		%r24,%r21,%r22
	subb		%r25,%r20,%r1
	comib,=		0,%r1,L$2
	ldo		-64(%r30),%r30

	add		%r22,%r23,%r22
	ldo		-1(%r28),%r28
L$2:
	bv		0(%r2)
	stws		%r22,0(%r26)

	.exit
	.procend