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
|
dnl ******************************************************************************
dnl Copyright 2009 Paul Zimmermann and Alexander Kruppa.
dnl
dnl This file is part of the ECM Library.
dnl
dnl The ECM Library is free software; you can redistribute it and/or modify
dnl it under the terms of the GNU Lesser General Public License as published by
dnl the Free Software Foundation; either version 3 of the License, or (at your
dnl option) any later version.
dnl
dnl The ECM Library is distributed in the hope that it will be useful, but
dnl WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
dnl or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
dnl License for more details.
dnl
dnl You should have received a copy of the GNU Lesser General Public License
dnl along with the ECM Library; see the file COPYING.LIB. If not, write to
dnl the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston,
dnl MA 02110-1301, USA.
dnl ******************************************************************************
define(C, `
dnl')
C mp_limb_t mulredc1(mp_limb_t * z, const mp_limb_t x, const mp_limb_t y,
C const mp_limb_t m, mp_limb_t inv_m);
C
C arguments:
C r3 : ptr to result z
C r4 : input x
C r5 : input y
C r6 : modulus m'
C r7 = -1/m mod 2^64
C
C final carry returned in r3
include(`config.m4')
GLOBL GSYM_PREFIX`'mulredc1
GLOBL .GSYM_PREFIX`'mulredc1
.section ".opd", "aw"
.align 3
GSYM_PREFIX`'mulredc1:
.quad .GSYM_PREFIX`'mulredc1, .TOC.@tocbase, 0
.size GSYM_PREFIX`'mulredc1, 24
TEXT
.align 5 C powerPC 32 byte alignment
.GSYM_PREFIX`'mulredc1:
mulld r8, r4, r5 C x*y low half T0
mulhdu r9, r4, r5 C x*y high half T1
mulld r0, r7, r8 C u = t0 * invm
mulld r10, r0, r6 C u*m low
mulhdu r11, r0, r6 C u*m high
addc r8, r8, r10 C x*y + u*m low (= zero)
adde r9, r9, r11 C result
std r9, 0(r3) C store in z
addze r3, r8 C return carry
blr
.size .GSYM_PREFIX`'mulredc1, .-.GSYM_PREFIX`'mulredc1
|
