sysdeps/mips/mips3/rshift.s

   1  # MIPS3 __mpn_rshift --
   2
   3  # Copyright (C) 1995 Free Software Foundation, Inc.
   4
   5  # This file is part of the GNU MP Library.
   6
   7  # The GNU MP Library is free software; you can redistribute it and/or modify
   8  # it under the terms of the GNU Library General Public License as published by
   9  # the Free Software Foundation; either version 2 of the License, or (at your
  10  # option) any later version.
  11
  12  # The GNU MP Library is distributed in the hope that it will be useful, but
  13  # WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  14  # or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Library General Public
  15  # License for more details.
  16
  17  # You should have received a copy of the GNU Library General Public License
  18  # along with the GNU MP Library; see the file COPYING.LIB.  If not, write to
  19  # the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
  20  # MA 02111-1307, USA.
  21
  22
  23  # INPUT PARAMETERS
  24  # res_ptr      $4
  25  # src_ptr      $5
  26  # size         $6
  27  # cnt          $7
  28
  29         .text
  30         .align  2
  31         .globl  __mpn_rshift
  32         .ent    __mpn_rshift
  33 __mpn_rshift:
  34         .set    noreorder
  35         .set    nomacro
  36
  37         ld      $10,0($5)       # load first limb
  38         dsubu   $13,$0,$7
  39         daddiu  $6,$6,-1
  40         and     $9,$6,4-1       # number of limbs in first loop
  41         beq     $9,$0,.L0       # if multiple of 4 limbs, skip first loop
  42          dsll   $2,$10,$13      # compute function result
  43
  44         dsubu   $6,$6,$9
  45
  46 .Loop0: ld      $3,8($5)
  47         daddiu  $4,$4,8
  48         daddiu  $5,$5,8
  49         daddiu  $9,$9,-1
  50         dsrl    $11,$10,$7
  51         dsll    $12,$3,$13
  52         move    $10,$3
  53         or      $8,$11,$12
  54         bne     $9,$0,.Loop0
  55          sd     $8,-8($4)
  56
  57 .L0:    beq     $6,$0,.Lend
  58          nop
  59
  60 .Loop:  ld      $3,8($5)
  61         daddiu  $4,$4,32
  62         daddiu  $6,$6,-4
  63         dsrl    $11,$10,$7
  64         dsll    $12,$3,$13
  65
  66         ld      $10,16($5)
  67         dsrl    $14,$3,$7
  68         or      $8,$11,$12
  69         sd      $8,-32($4)
  70         dsll    $9,$10,$13
  71
  72         ld      $3,24($5)
  73         dsrl    $11,$10,$7
  74         or      $8,$14,$9
  75         sd      $8,-24($4)
  76         dsll    $12,$3,$13
  77
  78         ld      $10,32($5)
  79         dsrl    $14,$3,$7
  80         or      $8,$11,$12
  81         sd      $8,-16($4)
  82         dsll    $9,$10,$13
  83
  84         daddiu  $5,$5,32
  85         or      $8,$14,$9
  86         bgtz    $6,.Loop
  87          sd     $8,-8($4)
  88
  89 .Lend:  dsrl    $8,$10,$7
  90         j       $31
  91         sd      $8,0($4)
  92         .end    __mpn_rshift