020c31e3a5a2b8a1839a285ad288c1e88a6686f4
[libgcrypt.git] / mpi / hppa1.1 / udiv-qrnnd.S
1 /* HP-PA  __udiv_qrnnd division support, used from longlong.h.
2  *        This version runs fast on PA 7000 and later.
3  *
4  *      Copyright (C) 1993, 1994, 1998, 
5  *                    2001, 2002, 2004 Free Software Foundation, Inc.
6  *
7  * This file is part of Libgcrypt.
8  *
9  * Libgcrypt is free software; you can redistribute it and/or modify
10  * it under the terms of the GNU Lesser General Public License as
11  * published by the Free Software Foundation; either version 2.1 of
12  * the License, or (at your option) any later version.
13  *
14  * Libgcrypt is distributed in the hope that it will be useful,
15  * but WITHOUT ANY WARRANTY; without even the implied warranty of
16  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
17  * GNU Lesser General Public License for more details.
18  *
19  * You should have received a copy of the GNU Lesser General Public
20  * License along with this program; if not, write to the Free Software
21  * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA
22  *
23  * Note: This code is heavily based on the GNU MP Library.
24  *       Actually it's the same code with only minor changes in the
25  *       way the data is stored; this is to support the abstraction
26  *       of an optional secure memory allocation which may be used
27  *       to avoid revealing of sensitive data due to paging etc.
28  */
29
30
31
32 /* INPUT PARAMETERS
33  *   rem_ptr       gr26
34  *   n1            gr25
35  *   n0            gr24
36  *   d             gr23
37  */
38
39         .data
40         .align          8
41         .label L$0000
42         .word           0x43f00000
43         .word           0x0
44         .code
45         .export         __udiv_qrnnd
46         .label          __udiv_qrnnd
47         .proc
48         .callinfo       frame=64,no_calls
49         .entry
50         ldo             64(%r30),%r30
51
52         stws            %r25,-16(0,%r30)        ; n_hi
53         stws            %r24,-12(0,%r30)        ; n_lo
54         stw             %r19,-32(%r30)
55         addil           LT%L$0000,%r19
56         ldw             RT%L$0000(%r1),%r1
57         fldds           -16(0,%r30),%fr5
58         stws            %r23,-12(0,%r30)
59         comib,<=        0,%r25,L$1
60         fcnvxf,dbl,dbl  %fr5,%fr5
61         fldds           0(0,%r1),%fr4
62         fadd,dbl        %fr4,%fr5,%fr5
63         .label  L$1
64         fcpy,sgl        %fr0,%fr6L
65         fldws           -12(0,%r30),%fr6R
66         fcnvxf,dbl,dbl  %fr6,%fr4
67
68         fdiv,dbl        %fr5,%fr4,%fr5
69
70         fcnvfx,dbl,dbl  %fr5,%fr4
71         fstws           %fr4R,-16(%r30)
72         xmpyu           %fr4R,%fr6R,%fr6
73         ldws            -16(%r30),%r28
74         fstds           %fr6,-16(0,%r30)
75         ldws            -12(0,%r30),%r21
76         ldws            -16(0,%r30),%r20
77         sub             %r24,%r21,%r22
78         subb            %r25,%r20,%r1
79         comib,=         0,%r1,L$2
80         ldo             -64(%r30),%r30
81
82         add             %r22,%r23,%r22
83         ldo             -1(%r28),%r28
84         .label  L$2
85         bv              0(%r2)
86         stws            %r22,0(0,%r26)
87
88         .exit
89         .procend
90