mpn/i960/addmul_1.s

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