Switched to GPLv3.
[gnupg.git] / mpi / i386 / mpih-mul3.S
1 /* i80386 submul_1 -- Multiply a limb vector with a limb and add
2  *                    the result to a second limb vector.
3  *
4  *      Copyright (C) 1992, 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 3 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, see <http://www.gnu.org/licenses/>.
21  *
22  * Note: This code is heavily based on the GNU MP Library.
23  *       Actually it's the same code with only minor changes in the
24  *       way the data is stored; this is to support the abstraction
25  *       of an optional secure memory allocation which may be used
26  *       to avoid revealing of sensitive data due to paging etc.
27  *       The GNU MP Library itself is published under the LGPL;
28  *       however I decided to publish this code under the plain GPL.
29  */
30
31
32 #include "sysdep.h"
33 #include "asm-syntax.h"
34
35
36 /*******************
37  * mpi_limb_t
38  * mpihelp_submul_1( mpi_ptr_t res_ptr,      (sp + 4)
39  *                   mpi_ptr_t s1_ptr,       (sp + 8)
40  *                   mpi_size_t s1_size,     (sp + 12)
41  *                   mpi_limb_t s2_limb)     (sp + 16)
42  */
43
44 #define res_ptr edi
45 #define s1_ptr  esi
46 #define size    ecx
47 #define s2_limb ebp
48
49         TEXT
50         ALIGN (3)
51         GLOBL   C_SYMBOL_NAME(mpihelp_submul_1)
52 C_SYMBOL_NAME(mpihelp_submul_1:)
53
54         INSN1(push,l    ,R(edi))
55         INSN1(push,l    ,R(esi))
56         INSN1(push,l    ,R(ebx))
57         INSN1(push,l    ,R(ebp))
58
59         INSN2(mov,l     ,R(res_ptr),MEM_DISP(esp,20))
60         INSN2(mov,l     ,R(s1_ptr),MEM_DISP(esp,24))
61         INSN2(mov,l     ,R(size),MEM_DISP(esp,28))
62         INSN2(mov,l     ,R(s2_limb),MEM_DISP(esp,32))
63
64         INSN2(lea,l     ,R(res_ptr),MEM_INDEX(res_ptr,size,4))
65         INSN2(lea,l     ,R(s1_ptr),MEM_INDEX(s1_ptr,size,4))
66         INSN1(neg,l     ,R(size))
67         INSN2(xor,l     ,R(ebx),R(ebx))
68         ALIGN (3)
69 Loop:
70         INSN2(mov,l     ,R(eax),MEM_INDEX(s1_ptr,size,4))
71         INSN1(mul,l     ,R(s2_limb))
72         INSN2(add,l     ,R(eax),R(ebx))
73         INSN2(adc,l     ,R(edx),$0)
74         INSN2(sub,l     ,MEM_INDEX(res_ptr,size,4),R(eax))
75         INSN2(adc,l     ,R(edx),$0)
76         INSN2(mov,l     ,R(ebx),R(edx))
77
78         INSN1(inc,l     ,R(size))
79         INSN1(jnz,      ,Loop)
80         INSN2(mov,l     ,R(eax),R(ebx))
81
82         INSN1(pop,l     ,R(ebp))
83         INSN1(pop,l     ,R(ebx))
84         INSN1(pop,l     ,R(esi))
85         INSN1(pop,l     ,R(edi))
86         ret
87