sysdeps/powerpc/mul_1.s

   1  # Multiply a limb vector by a limb, for PowerPC.
   2  # Copyright (C) 1993, 1994, 1995, 1997 Free Software Foundation, Inc.
   3  # This file is part of the GNU C Library.
   4  #
   5  # The GNU C Library is free software; you can redistribute it and/or
   6  # modify it under the terms of the GNU Library General Public License as
   7  # published by the Free Software Foundation; either version 2 of the
   8  # License, or (at your option) any later version.
   9  #
  10  # The GNU C Library is distributed in the hope that it will be useful,
  11  # but WITHOUT ANY WARRANTY; without even the implied warranty of
  12  # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13  # Library General Public License for more details.
  14  #
  15  # You should have received a copy of the GNU Library General Public
  16  # License along with the GNU C Library; see the file COPYING.LIB.  If not,
  17  # write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
  18  # Boston, MA 02111-1307, USA.
  19
  20  # mp_limb_t mpn_mul_1 (mp_ptr res_ptr, mp_srcptr s1_ptr,
  21  #                      mp_size_t s1_size, mp_limb_t s2_limb)
  22  # Calculate s1*s2 and put result in res_ptr; return carry.
  23
  24         .align 2
  25         .globl __mpn_mul_1
  26         .type    __mpn_mul_1,@function
  27
  28 __mpn_mul_1:
  29         mtctr   %r5
  30
  31         lwz     %r0,0(%r4)
  32         mullw   %r7,%r0,%r6
  33         mulhwu  %r10,%r0,%r6
  34         addi    %r3,%r3,-4              # adjust res_ptr
  35         addic   %r5,%r5,0               # clear cy with dummy insn
  36         bdz     Lend
  37
  38 Loop:   lwzu    %r0,4(%r4)
  39         stwu    %r7,4(%r3)
  40         mullw   %r8,%r0,%r6
  41         adde    %r7,%r8,%r10
  42         mulhwu  %r10,%r0,%r6
  43         bdnz    Loop
  44
  45 Lend:   stw     %r7,4(%r3)
  46         addze   %r3,%r10
  47         blr