~ubuntu-branches/ubuntu/intrepid/ecl/intrepid

« back to all changes in this revision

Viewing changes to src/gmp/mpn/powerpc64/mode64/mul_1.asm

  • Committer: Bazaar Package Importer
  • Author(s): Peter Van Eynde
  • Date: 2007-04-09 11:51:51 UTC
  • mfrom: (1.1.3 upstream)
  • Revision ID: james.westby@ubuntu.com-20070409115151-ql8cr0kalzx1jmla
Tags: 0.9i-20070324-2
Upload to unstable. 

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
dnl  PowerPC-64 mpn_mul_1 -- Multiply a limb vector with a limb and store
 
2
dnl  the result in a second limb vector.
 
3
 
 
4
dnl  Copyright 1999, 2000, 2001 Free Software Foundation, Inc.
 
5
 
 
6
dnl  This file is part of the GNU MP Library.
 
7
 
 
8
dnl  The GNU MP Library is free software; you can redistribute it and/or modify
 
9
dnl  it under the terms of the GNU Lesser General Public License as published by
 
10
dnl  the Free Software Foundation; either version 2.1 of the License, or (at your
 
11
dnl  option) any later version.
 
12
 
 
13
dnl  The GNU MP Library is distributed in the hope that it will be useful, but
 
14
dnl  WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
 
15
dnl  or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
 
16
dnl  License for more details.
 
17
 
 
18
dnl  You should have received a copy of the GNU Lesser General Public License
 
19
dnl  along with the GNU MP Library; see the file COPYING.LIB.  If not, write
 
20
dnl  to the Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
 
21
dnl  Boston, MA 02110-1301, USA.
 
22
 
 
23
include(`../config.m4')
 
24
 
 
25
C INPUT PARAMETERS
 
26
C res_ptr       r3
 
27
C s1_ptr        r4
 
28
C size          r5
 
29
C s2_limb       r6
 
30
C cy_limb       r7
 
31
 
 
32
C PPC630: 6 to 18 cycles/limb, depending on multiplier.  This cannot be
 
33
C improved unless floating-point operations are used instead of the slow
 
34
C mulld/mulhdu.
 
35
 
 
36
ASM_START()
 
37
PROLOGUE(mpn_mul_1)
 
38
        li      r7,0            C cy_limb = 0
 
39
 
 
40
PROLOGUE(mpn_mul_1c)
 
41
        mtctr   r5
 
42
        addic   r0,r0,0
 
43
        addi    r3, r3, -8
 
44
        addi    r4, r4, -8
 
45
L(top): ldu     r0,8(r4)
 
46
        mulld   r9,r0,r6
 
47
        adde    r9,r9,r7
 
48
        mulhdu  r7,r0,r6
 
49
        stdu    r9,8(r3)
 
50
        bdnz    L(top)
 
51
 
 
52
        addze   r3,r7
 
53
        blr
 
54
EPILOGUE(mpn_mul_1)
 
55
EPILOGUE(mpn_mul_1c)