source/libs/gmp/gmp-src/mpn/pa64/lshift.asm

   1 dnl  HP-PA 2.0 mpn_lshift -- Left shift.
   2
   3 dnl  Copyright 1997, 2000, 2002, 2003 Free Software Foundation, Inc.
   4
   5 dnl  This file is part of the GNU MP Library.
   6 dnl
   7 dnl  The GNU MP Library is free software; you can redistribute it and/or modify
   8 dnl  it under the terms of either:
   9 dnl
  10 dnl    * the GNU Lesser General Public License as published by the Free
  11 dnl      Software Foundation; either version 3 of the License, or (at your
  12 dnl      option) any later version.
  13 dnl
  14 dnl  or
  15 dnl
  16 dnl    * the GNU General Public License as published by the Free Software
  17 dnl      Foundation; either version 2 of the License, or (at your option) any
  18 dnl      later version.
  19 dnl
  20 dnl  or both in parallel, as here.
  21 dnl
  22 dnl  The GNU MP Library is distributed in the hope that it will be useful, but
  23 dnl  WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  24 dnl  or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  25 dnl  for more details.
  26 dnl
  27 dnl  You should have received copies of the GNU General Public License and the
  28 dnl  GNU Lesser General Public License along with the GNU MP Library.  If not,
  29 dnl  see https://www.gnu.org/licenses/.
  30
  31
  32 dnl  This runs at 1.5 cycles/limb on PA8000 and 1.0 cycles/limb on PA8500.
  33
  34 include(`../config.m4')
  35
  36 dnl  INPUT PARAMETERS
  37 define(`rp',`%r26')
  38 define(`up',`%r25')
  39 define(`n',`%r24')
  40 define(`cnt',`%r23')
  41
  42 ifdef(`HAVE_ABI_2_0w',
  43 `       .level  2.0w
  44 ',`     .level  2.0
  45 ')
  46 PROLOGUE(mpn_lshift)
  47         shladd          n, 3, up, up
  48         shladd          n, 3, rp, rp
  49         subi            64, cnt, cnt
  50         mtsar           cnt
  51         ldd             -8(up), %r21
  52         addib,=         -1, n, L(end)
  53         shrpd           %r0, %r21, %sar, %r29   C compute carry out limb
  54         depw,z          n, 31, 3, %r28          C r28 = (size & 7)
  55         sub             %r0, n, %r22
  56         depw,z          %r22, 28, 3, %r22       C r22 = 8 * (-size & 7)
  57         add             up, %r22, up            C offset up
  58         blr             %r28, %r0               C branch into jump table
  59         add             rp, %r22, rp            C offset rp
  60         b               L(0)
  61         nop
  62         b               L(1)
  63         copy            %r21, %r20
  64         b               L(2)
  65         nop
  66         b               L(3)
  67         copy            %r21, %r20
  68         b               L(4)
  69         nop
  70         b               L(5)
  71         copy            %r21, %r20
  72         b               L(6)
  73         nop
  74         b               L(7)
  75         copy            %r21, %r20
  76
  77 LDEF(loop)
  78 LDEF(0) ldd             -16(up), %r20
  79         shrpd           %r21, %r20, %sar, %r21
  80         std             %r21, -8(rp)
  81 LDEF(7) ldd             -24(up), %r21
  82         shrpd           %r20, %r21, %sar, %r20
  83         std             %r20, -16(rp)
  84 LDEF(6) ldd             -32(up), %r20
  85         shrpd           %r21, %r20, %sar, %r21
  86         std             %r21, -24(rp)
  87 LDEF(5) ldd             -40(up), %r21
  88         shrpd           %r20, %r21, %sar, %r20
  89         std             %r20, -32(rp)
  90 LDEF(4) ldd             -48(up), %r20
  91         shrpd           %r21, %r20, %sar, %r21
  92         std             %r21, -40(rp)
  93 LDEF(3) ldd             -56(up), %r21
  94         shrpd           %r20, %r21, %sar, %r20
  95         std             %r20, -48(rp)
  96 LDEF(2) ldd             -64(up), %r20
  97         shrpd           %r21, %r20, %sar, %r21
  98         std             %r21, -56(rp)
  99 LDEF(1) ldd             -72(up), %r21
 100         ldo             -64(up), up
 101         shrpd           %r20, %r21, %sar, %r20
 102         std             %r20, -64(rp)
 103         addib,>         -8, n, L(loop)
 104         ldo             -64(rp), rp
 105
 106 LDEF(end)
 107         shrpd           %r21, %r0, %sar, %r21
 108         std             %r21, -8(rp)
 109         bve             (%r2)
 110 ifdef(`HAVE_ABI_2_0w',
 111 `       copy            %r29,%r28
 112 ',`     extrd,u         %r29, 31, 32, %r28
 113 ')
 114 EPILOGUE(mpn_lshift)