diff options
author | Werner Koch <wk@gnupg.org> | 1998-11-27 11:42:48 +0000 |
---|---|---|
committer | Werner Koch <wk@gnupg.org> | 1998-11-27 11:42:48 +0000 |
commit | 0b9ccc4150f5f6c5310267f54649612b905487b2 (patch) | |
tree | 8d1249d2b7a3202cf506046917fc599f463926b2 /mpi/powerpc32/mpih-mul3.S | |
parent | 37d536064fdd448553612b3a80de7058e9fd8756 (diff) | |
download | libgcrypt-0b9ccc4150f5f6c5310267f54649612b905487b2.tar.gz |
See ChangeLog: Fri Nov 27 12:39:29 CET 1998 Werner Koch
Diffstat (limited to 'mpi/powerpc32/mpih-mul3.S')
-rw-r--r-- | mpi/powerpc32/mpih-mul3.S | 56 |
1 files changed, 56 insertions, 0 deletions
diff --git a/mpi/powerpc32/mpih-mul3.S b/mpi/powerpc32/mpih-mul3.S index 46e85f68..a4a4d163 100644 --- a/mpi/powerpc32/mpih-mul3.S +++ b/mpi/powerpc32/mpih-mul3.S @@ -20,6 +20,11 @@ * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA */ +#include "sysdep.h" +#include "asm-syntax.h" + + +#ifndef USE_PPC_PATCHES /******************* * mpi_limb_t @@ -71,4 +76,55 @@ Lend: stw 8,4(3) addze 3,10 blr +#else +/* Multiply a limb vector by a single limb, for PowerPC. + Copyright (C) 1993, 1994, 1995, 1997 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Library General Public License as + published by the Free Software Foundation; either version 2 of the + License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Library General Public License for more details. + + You should have received a copy of the GNU Library General Public + License along with the GNU C Library; see the file COPYING.LIB. If not, + write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, + Boston, MA 02111-1307, USA. */ + +/* mp_limb_t mpn_submul_1 (mp_ptr res_ptr, mp_srcptr s1_ptr, + mp_size_t s1_size, mp_limb_t s2_limb) + Calculate res-s1*s2 and put result back in res; return carry. */ + +ENTRY(_mpihelp_submul_1) + mtctr %r5 + + lwz %r0,0(%r4) + mullw %r7,%r0,%r6 + mulhwu %r10,%r0,%r6 + lwz %r9,0(%r3) + subf %r8,%r7,%r9 + addc %r7,%r7,%r8 # invert cy (r7 is junk) + addi %r3,%r3,-4 # adjust res_ptr + bdz 1f + +0: lwzu %r0,4(%r4) + stwu %r8,4(%r3) + mullw %r8,%r0,%r6 + adde %r7,%r8,%r10 + mulhwu %r10,%r0,%r6 + lwz %r9,4(%r3) + addze %r10,%r10 + subf %r8,%r7,%r9 + addc %r7,%r7,%r8 # invert cy (r7 is junk) + bdnz 0b +1: stw %r8,4(%r3) + addze %r3,%r10 + blr +END(_mpihelp_submul_1) +#endif |