![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to lshift.asm CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / gmp / mpn / m68k|
Return to lshift.asm CVS log | Up to [local] / OpenXM_contrib / gmp / mpn / m68k |
1.1 ohara 1: dnl mc68020 mpn_lshift -- mpn left shift.
2: dnl
3: dnl cycles/limb
4: dnl shift==1 shift>1
5: dnl 68040: 3.5 9.5
6:
7: dnl Copyright 1996, 1999, 2000, 2001, 2002 Free Software Foundation, Inc.
8: dnl
9: dnl This file is part of the GNU MP Library.
10: dnl
11: dnl The GNU MP Library is free software; you can redistribute it and/or
12: dnl modify it under the terms of the GNU Lesser General Public License as
13: dnl published by the Free Software Foundation; either version 2.1 of the
14: dnl License, or (at your option) any later version.
15: dnl
16: dnl The GNU MP Library is distributed in the hope that it will be useful,
17: dnl but WITHOUT ANY WARRANTY; without even the implied warranty of
18: dnl MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
19: dnl Lesser General Public License for more details.
20: dnl
21: dnl You should have received a copy of the GNU Lesser General Public
22: dnl License along with the GNU MP Library; see the file COPYING.LIB. If
23: dnl not, write to the Free Software Foundation, Inc., 59 Temple Place -
24: dnl Suite 330, Boston, MA 02111-1307, USA.
25:
26: include(`../config.m4')
27:
28:
29: C INPUT PARAMETERS
30: C res_ptr (sp + 4)
31: C s_ptr (sp + 8)
32: C s_size (sp + 16)
33: C cnt (sp + 12)
34:
35: define(res_ptr, `a1')
36: define(s_ptr, `a0')
37: define(s_size, `d6')
38: define(cnt, `d4')
39:
40: PROLOGUE(mpn_lshift)
41: C Save used registers on the stack.
42: moveml d2-d6/a2, M(-,sp)
43:
44: C Copy the arguments to registers.
45: movel M(sp,28), res_ptr
46: movel M(sp,32), s_ptr
47: movel M(sp,36), s_size
48: movel M(sp,40), cnt
49:
50: moveql #1, d5
51: cmpl d5, cnt
52: bne L(Lnormal)
53: cmpl s_ptr, res_ptr
54: bls L(Lspecial) C jump if s_ptr >= res_ptr
55:
56: ifelse(scale_available_p,1,`
57: lea M(s_ptr,s_size,l,4), a2
58: ',`
59: movel s_size, d0
60: asll #2, d0
61: lea M(s_ptr,d0,l), a2
62: ')
63: cmpl res_ptr, a2
64: bls L(Lspecial) C jump if res_ptr >= s_ptr + s_size
65:
66: L(Lnormal):
67: moveql #32, d5
68: subl cnt, d5
69:
70: ifelse(scale_available_p,1,`
71: lea M(s_ptr,s_size,l,4), s_ptr
72: lea M(res_ptr,s_size,l,4), res_ptr
73: ',`
74: movel s_size, d0
75: asll #2, d0
76: addl d0, s_ptr
77: addl d0, res_ptr
78: ')
79: movel M(-,s_ptr), d2
80: movel d2, d0
81: lsrl d5, d0 C compute carry limb
82:
83: lsll cnt, d2
84: movel d2, d1
85: subql #1, s_size
86: beq L(Lend)
87: lsrl #1, s_size
88: bcs L(L1)
89: subql #1, s_size
90:
91: L(Loop:)
92: movel M(-,s_ptr), d2
93: movel d2, d3
94: lsrl d5, d3
95: orl d3, d1
96: movel d1, M(-,res_ptr)
97: lsll cnt, d2
98: L(L1:)
99: movel M(-,s_ptr), d1
100: movel d1, d3
101: lsrl d5, d3
102: orl d3, d2
103: movel d2, M(-,res_ptr)
104: lsll cnt, d1
105:
106: dbf s_size, L(Loop)
107: subl #0x10000, s_size
108: bcc L(Loop)
109:
110: L(Lend:)
111: movel d1, M(-,res_ptr) C store least significant limb
112:
113: C Restore used registers from stack frame.
114: moveml M(sp,+), d2-d6/a2
115: rts
116:
117: C We loop from least significant end of the arrays, which is only
118: C permissable if the source and destination don't overlap, since the
119: C function is documented to work for overlapping source and destination.
120:
121: L(Lspecial):
122: clrl d0 C initialize carry
123: eorw #1, s_size
124: lsrl #1, s_size
125: bcc L(LL1)
126: subql #1, s_size
127:
128: L(LLoop):
129: movel M(s_ptr,+), d2
130: addxl d2, d2
131: movel d2, M(res_ptr,+)
132: L(LL1):
133: movel M(s_ptr,+), d2
134: addxl d2, d2
135: movel d2, M(res_ptr,+)
136:
137: dbf s_size, L(LLoop)
138: addxl d0, d0 C save cy in lsb
139: subl #0x10000, s_size
140: bcs L(LLend)
141: lsrl #1, d0 C restore cy
142: bra L(LLoop)
143:
144: L(LLend):
145: C Restore used registers from stack frame.
146: moveml M(sp,+), d2-d6/a2
147: rts
148:
149: EPILOGUE(mpn_lshift)