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 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 3 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.  If not, see http://www.gnu.org/licenses/.
  20
  21 .text
  22         .align  16
  23 .globl  ___gmpn_mul_1
  24 ___gmpn_mul_1:
  25         subq    $8,sp
  26         storw   r6,(sp)
  27         loadw   12(sp),r2
  28         loadw   16(sp),r3
  29         loadq   $0,r6           ; clear carry limb
  30
  31 .Loop:  loadw   (r1),r4
  32         mulwux  r3,r4
  33         addw    r6,r4           ; add old carry limb into low product limb
  34         loadq   $0,r6
  35         addwc   r5,r6           ; propagate cy into high product limb
  36         storw   r4,(r0)
  37         addq    $4,r0
  38         addq    $4,r1
  39         subq    $1,r2
  40         brne    .Loop
  41
  42         movw    r6,r0
  43         loadw   0(sp),r6
  44         addq    $8,sp
  45         ret     sp