← Back to branch summary

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

« back to all changes in this revision

Viewing changes to src/gmp/mpn/powerpc64/submul_1.asm

  • Committer: Bazaar Package Importer
  • Author(s): Peter Van Eynde
  • Date: 2006-05-17 02:46:26 UTC
  • Revision ID: james.westby@ubuntu.com-20060517024626-lljr08ftv9g9vefl
Tags: upstream-0.9h-20060510
Import upstream version 0.9h-20060510

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/gmp/mpn/powerpc64/submul_1.asm
 
1
dnl  PowerPC-64 mpn_submul_1 -- Multiply a limb vector with a limb and subtract
 
2
dnl  the result from 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
 
10
dnl  by the Free Software Foundation; either version 2.1 of the License, or (at
 
11
dnl  your 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 to
 
20
dnl  the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
 
21
dnl  MA 02111-1307, 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_submul_1)
 
38
        li      r7,0                    C cy_limb = 0
 
39
 
 
40
PROLOGUE(mpn_submul_1c)
 
41
        mtctr   r5
 
42
        addic   r0,r0,0
 
43
        addi    r3, r3, -8
 
44
        addi    r4, r4, -8
 
45
.Loop:
 
46
        ldu     r0,8(r4)
 
47
        ld      r10,8(r3)
 
48
        mulld   r9,r0,r6
 
49
        adde    r9,r9,r7
 
50
        mulhdu  r7,r0,r6
 
51
        addze   r7,r7
 
52
        subfc   r9,r9,r10
 
53
        stdu    r9,8(r3)
 
54
        subfe   r11,r11,r11             C invert ...
 
55
        addic   r11,r11,1               C ... carry
 
56
        bdnz    .Loop
 
57
 
 
58
        addze   r3,r7
 
59
        blr
 
60
EPILOGUE(mpn_submul_1)
 
61
EPILOGUE(mpn_submul_1c)