rts/gmp/mpn/clipper/mul_1.s

   1 ; Clipper __gmpn_mul_1 -- Multiply a limb vector with a limb and store
   2 ; the result in a second limb vector.
   3
   4 ; Copyright (C) 1995, 2000 Free Software Foundation, Inc.
   5
   6 ; This file is part of the GNU MP Library.
   7
   8 ; The GNU MP Library is free software; you can redistribute it and/or modify
   9 ; it under the terms of the GNU Lesser General Public License as published by
  10 ; the Free Software Foundation; either version 2.1 of the License, or (at your
  11 ; option) any later version.
  12
  13 ; The GNU MP Library is distributed in the hope that it will be useful, but
  14 ; WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  15 ; or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  16 ; License for more details.
  17
  18 ; You should have received a copy of the GNU Lesser General Public License
  19 ; along with the GNU MP Library; see the file COPYING.LIB.  If not, write to
  20 ; the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
  21 ; MA 02111-1307, USA.
  22
  23 .text
  24         .align  16
  25 .globl  ___gmpn_mul_1
  26 ___gmpn_mul_1:
  27         subq    $8,sp
  28         storw   r6,(sp)
  29         loadw   12(sp),r2
  30         loadw   16(sp),r3
  31         loadq   $0,r6           ; clear carry limb
  32
  33 .Loop:  loadw   (r1),r4
  34         mulwux  r3,r4
  35         addw    r6,r4           ; add old carry limb into low product limb
  36         loadq   $0,r6
  37         addwc   r5,r6           ; propagate cy into high product limb
  38         storw   r4,(r0)
  39         addq    $4,r0
  40         addq    $4,r1
  41         subq    $1,r2
  42         brne    .Loop
  43
  44         movw    r6,r0
  45         loadw   0(sp),r6
  46         addq    $8,sp
  47         ret     sp