Annotation of OpenXM_contrib/gmp/mpn/sparc32/rshift.S, Revision 1.1.1.1
1.1 maekawa 1: ! sparc __mpn_rshift --
2:
3: ! Copyright (C) 1995, 1996 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 %o0
25: ! src_ptr %o1
26: ! size %o2
27: ! cnt %o3
28:
29: #include "sysdep.h"
30:
31: .text
32: .align 4
33: .global C_SYMBOL_NAME(__mpn_rshift)
34: C_SYMBOL_NAME(__mpn_rshift):
35: ld [%o1],%g2 ! load first limb
36: sub %g0,%o3,%o5 ! negate shift count
37: add %o2,-1,%o2
38: andcc %o2,4-1,%g4 ! number of limbs in first loop
39: sll %g2,%o5,%g1 ! compute function result
40: be L0 ! if multiple of 4 limbs, skip first loop
41: st %g1,[%sp+80]
42:
43: sub %o2,%g4,%o2 ! adjust count for main loop
44:
45: Loop0: ld [%o1+4],%g3
46: add %o0,4,%o0
47: add %o1,4,%o1
48: addcc %g4,-1,%g4
49: srl %g2,%o3,%o4
50: sll %g3,%o5,%g1
51: mov %g3,%g2
52: or %o4,%g1,%o4
53: bne Loop0
54: st %o4,[%o0-4]
55:
56: L0: tst %o2
57: be Lend
58: nop
59:
60: Loop: ld [%o1+4],%g3
61: add %o0,16,%o0
62: addcc %o2,-4,%o2
63: srl %g2,%o3,%o4
64: sll %g3,%o5,%g1
65:
66: ld [%o1+8],%g2
67: srl %g3,%o3,%g4
68: or %o4,%g1,%o4
69: st %o4,[%o0-16]
70: sll %g2,%o5,%g1
71:
72: ld [%o1+12],%g3
73: srl %g2,%o3,%o4
74: or %g4,%g1,%g4
75: st %g4,[%o0-12]
76: sll %g3,%o5,%g1
77:
78: ld [%o1+16],%g2
79: srl %g3,%o3,%g4
80: or %o4,%g1,%o4
81: st %o4,[%o0-8]
82: sll %g2,%o5,%g1
83:
84: add %o1,16,%o1
85: or %g4,%g1,%g4
86: bne Loop
87: st %g4,[%o0-4]
88:
89: Lend: srl %g2,%o3,%g2
90: st %g2,[%o0-0]
91: retl
92: ld [%sp+80],%o0
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>