From 11da511c784eca003deb90c23570f0873954e0de Mon Sep 17 00:00:00 2001 From: Duncan Wilkie Date: Sat, 18 Nov 2023 06:11:09 -0600 Subject: Initial commit. --- gmp-6.3.0/mpn/power/lshift.asm | 61 ++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 61 insertions(+) create mode 100644 gmp-6.3.0/mpn/power/lshift.asm (limited to 'gmp-6.3.0/mpn/power/lshift.asm') diff --git a/gmp-6.3.0/mpn/power/lshift.asm b/gmp-6.3.0/mpn/power/lshift.asm new file mode 100644 index 0000000..efa2105 --- /dev/null +++ b/gmp-6.3.0/mpn/power/lshift.asm @@ -0,0 +1,61 @@ +dnl IBM POWER mpn_lshift -- Shift a number left. + +dnl Copyright 1992, 1994, 1999-2001 Free Software Foundation, Inc. + +dnl This file is part of the GNU MP Library. +dnl +dnl The GNU MP Library is free software; you can redistribute it and/or modify +dnl it under the terms of either: +dnl +dnl * the GNU Lesser General Public License as published by the Free +dnl Software Foundation; either version 3 of the License, or (at your +dnl option) any later version. +dnl +dnl or +dnl +dnl * the GNU General Public License as published by the Free Software +dnl Foundation; either version 2 of the License, or (at your option) any +dnl later version. +dnl +dnl or both in parallel, as here. +dnl +dnl The GNU MP Library is distributed in the hope that it will be useful, but +dnl WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY +dnl or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License +dnl for more details. +dnl +dnl You should have received copies of the GNU General Public License and the +dnl GNU Lesser General Public License along with the GNU MP Library. If not, +dnl see https://www.gnu.org/licenses/. + + +dnl INPUT PARAMETERS +dnl res_ptr r3 +dnl s_ptr r4 +dnl size r5 +dnl cnt r6 + +include(`../config.m4') + +ASM_START() +PROLOGUE(mpn_lshift) + sli 0,5,2 + cax 9,3,0 + cax 4,4,0 + sfi 8,6,32 + mtctr 5 C put limb count in CTR loop register + lu 0,-4(4) C read most significant limb + sre 3,0,8 C compute carry out limb, and init MQ register + bdz Lend2 C if just one limb, skip loop + lu 0,-4(4) C read 2:nd most significant limb + sreq 7,0,8 C compute most significant limb of result + bdz Lend C if just two limb, skip loop +Loop: lu 0,-4(4) C load next lower limb + stu 7,-4(9) C store previous result during read latency + sreq 7,0,8 C compute result limb + bdn Loop C loop back until CTR is zero +Lend: stu 7,-4(9) C store 2:nd least significant limb +Lend2: sle 7,0,6 C compute least significant limb + st 7,-4(9) C store it + br +EPILOGUE(mpn_lshift) -- cgit v1.2.3