Annotation of OpenXM_contrib/gmp/mpn/power/rshift.s, Revision 1.1.1.2
1.1.1.2 ! maekawa 1: # IBM POWER __gmpn_rshift --
1.1 maekawa 2:
1.1.1.2 ! maekawa 3: # Copyright (C) 1992, 1994, 1999, 2000 Free Software Foundation, Inc.
1.1 maekawa 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
1.1.1.2 ! maekawa 8: # it under the terms of the GNU Lesser General Public License as published by
! 9: # the Free Software Foundation; either version 2.1 of the License, or (at your
1.1 maekawa 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
1.1.1.2 ! maekawa 14: # or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
1.1 maekawa 15: # License for more details.
16:
1.1.1.2 ! maekawa 17: # You should have received a copy of the GNU Lesser General Public License
1.1 maekawa 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 r3
25: # s_ptr r4
26: # size r5
27: # cnt r6
28:
29: .toc
1.1.1.2 ! maekawa 30: .globl __gmpn_rshift
! 31: .globl .__gmpn_rshift
! 32: .csect __gmpn_rshift[DS]
! 33: __gmpn_rshift:
! 34: .long .__gmpn_rshift, TOC[tc0], 0
! 35: .csect .text[PR]
! 36: .align 2
! 37: .__gmpn_rshift:
1.1 maekawa 38: sfi 8,6,32
39: mtctr 5 # put limb count in CTR loop register
40: l 0,0(4) # read least significant limb
41: ai 9,3,-4 # adjust res_ptr since it's offset in the stu:s
42: sle 3,0,8 # compute carry limb, and init MQ register
43: bdz Lend2 # if just one limb, skip loop
44: lu 0,4(4) # read 2:nd least significant limb
45: sleq 7,0,8 # compute least significant limb of result
46: bdz Lend # if just two limb, skip loop
47: Loop: lu 0,4(4) # load next higher limb
48: stu 7,4(9) # store previous result during read latency
49: sleq 7,0,8 # compute result limb
50: bdn Loop # loop back until CTR is zero
51: Lend: stu 7,4(9) # store 2:nd most significant limb
52: Lend2: sre 7,0,6 # compute most significant limb
53: st 7,4(9) # store it" \
54: br
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>