rts/gmp/mpn/i960/mul_1.s

   1 # I960 __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  4
  25         .globl  ___gmpn_mul_1
  26 ___gmpn_mul_1:
  27         subo    g2,0,g2
  28         shlo    2,g2,g4
  29         subo    g4,g1,g1
  30         subo    g4,g0,g13
  31         mov     0,g0
  32
  33         cmpo    1,0             # clear C bit on AC.cc
  34
  35 Loop:   ld      (g1)[g2*4],g5
  36         emul    g3,g5,g6
  37
  38         addc    g0,g6,g6        # relies on that C bit is clear
  39         st      g6,(g13)[g2*4]
  40         addc    0,g7,g0
  41
  42         addo    g2,1,g2
  43         cmpobne 0,g2,Loop       # when branch is taken, clears C bit
  44
  45         ret