~ubuntu-branches/ubuntu/precise/gnupg2/precise-proposed

« back to all changes in this revision

Viewing changes to mpi/hppa1.1/mpih-mul3.S

  • Committer: Bazaar Package Importer
  • Author(s): Matthias Urlichs
  • Date: 2006-01-24 04:31:42 UTC
  • mfrom: (1.1.2 upstream)
  • Revision ID: james.westby@ubuntu.com-20060124043142-pbg192or6qxv3yk2
Tags: 1.9.20-1
* New Upstream version. Closes:#306890,#344530
  * Closes:#320490: gpg-protect-tool fails to decrypt PKCS-12 files 
* Depend on libopensc2-dev, not -1-. Closes:#348106

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
/* hppa1.1   submul_1 -- Multiply a limb vector with a limb and add
 
2
 *                       the result to a second limb vector.
 
3
 *
 
4
 *      Copyright (C) 1992, 1993, 1994, 1998, 
 
5
 *                    2001 Free Software Foundation, Inc.
 
6
 *
 
7
 * This file is part of GnuPG.
 
8
 *
 
9
 * GnuPG is free software; you can redistribute it and/or modify
 
10
 * it under the terms of the GNU General Public License as published by
 
11
 * the Free Software Foundation; either version 2 of the License, or
 
12
 * (at your option) any later version.
 
13
 *
 
14
 * GnuPG 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 General Public License for more details.
 
18
 *
 
19
 * You should have received a copy of the GNU General Public License
 
20
 * 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
 *       The GNU MP Library itself is published under the LGPL;
 
29
 *       however I decided to publish this code under the plain GPL.
 
30
 */
 
31
 
 
32
 
 
33
/*******************
 
34
 * mpi_limb_t
 
35
 * mpihelp_submul_1( mpi_ptr_t res_ptr,      (r26)
 
36
 *                   mpi_ptr_t s1_ptr,       (r25)
 
37
 *                   mpi_size_t s1_size,     (r24)
 
38
 *                   mpi_limb_t s2_limb)     (r23)
 
39
 *
 
40
 *
 
41
 * This runs at 12 cycles/limb on a PA7000.  With the used instructions, it
 
42
 * can not become faster due to data cache contention after a store.  On the
 
43
 * PA7100 it runs at 11 cycles/limb, and that can not be improved either,
 
44
 * since only the xmpyu does not need the integer pipeline, so the only
 
45
 * dual-issue we will get are addc+xmpyu.  Unrolling could gain a cycle/limb
 
46
 * on the PA7100.
 
47
 *
 
48
 * There are some ideas described in mul1.S that applies to this code too.
 
49
 *
 
50
 * It seems possible to make this run as fast as addmul_1, if we use
 
51
 *       sub,>>= %r29,%r19,%r22
 
52
 *       addi    1,%r28,%r28
 
53
 * but that requires reworking the hairy software pipeline...
 
54
 */
 
55
 
 
56
 
 
57
        .code
 
58
        .export         mpihelp_submul_1
 
59
        .label          mpihelp_submul_1
 
60
        .proc
 
61
        .callinfo       frame=64,no_calls
 
62
        .entry
 
63
 
 
64
        ldo             64(%r30),%r30
 
65
        fldws,ma        4(%r25),%fr5
 
66
        stw             %r23,-16(%r30)          ; move s2_limb ...
 
67
        addib,=         -1,%r24,L$just_one_limb
 
68
         fldws          -16(%r30),%fr4          ; ... into fr4
 
69
        add             %r0,%r0,%r0             ; clear carry
 
70
        xmpyu           %fr4,%fr5,%fr6
 
71
        fldws,ma        4(%r25),%fr7
 
72
        fstds           %fr6,-16(%r30)
 
73
        xmpyu           %fr4,%fr7,%fr8
 
74
        ldw             -12(%r30),%r19          ; least significant limb in product
 
75
        ldw             -16(%r30),%r28
 
76
 
 
77
        fstds           %fr8,-16(%r30)
 
78
        addib,=         -1,%r24,L$end
 
79
         ldw            -12(%r30),%r1
 
80
 
 
81
; Main loop
 
82
        .label  L$loop
 
83
        ldws            0(%r26),%r29
 
84
        fldws,ma        4(%r25),%fr5
 
85
        sub             %r29,%r19,%r22
 
86
        add             %r22,%r19,%r0
 
87
        stws,ma         %r22,4(%r26)
 
88
        addc            %r28,%r1,%r19
 
89
        xmpyu           %fr4,%fr5,%fr6
 
90
        ldw             -16(%r30),%r28
 
91
        fstds           %fr6,-16(%r30)
 
92
        addc            %r0,%r28,%r28
 
93
        addib,<>        -1,%r24,L$loop
 
94
         ldw            -12(%r30),%r1
 
95
 
 
96
        .label  L$end
 
97
        ldw             0(%r26),%r29
 
98
        sub             %r29,%r19,%r22
 
99
        add             %r22,%r19,%r0
 
100
        stws,ma         %r22,4(%r26)
 
101
        addc            %r28,%r1,%r19
 
102
        ldw             -16(%r30),%r28
 
103
        ldws            0(%r26),%r29
 
104
        addc            %r0,%r28,%r28
 
105
        sub             %r29,%r19,%r22
 
106
        add             %r22,%r19,%r0
 
107
        stws,ma         %r22,4(%r26)
 
108
        addc            %r0,%r28,%r28
 
109
        bv              0(%r2)
 
110
         ldo            -64(%r30),%r30
 
111
 
 
112
        .label  L$just_one_limb
 
113
        xmpyu           %fr4,%fr5,%fr6
 
114
        ldw             0(%r26),%r29
 
115
        fstds           %fr6,-16(%r30)
 
116
        ldw             -12(%r30),%r1
 
117
        ldw             -16(%r30),%r28
 
118
        sub             %r29,%r1,%r22
 
119
        add             %r22,%r1,%r0
 
120
        stw             %r22,0(%r26)
 
121
        addc            %r0,%r28,%r28
 
122
        bv              0(%r2)
 
123
         ldo            -64(%r30),%r30
 
124
 
 
125
        .exit
 
126
        .procend
 
127