1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
|
/* Multiple versions of strrchr
All versions must be listed in ifunc-impl-list.c.
Copyright (C) 2009-2013 Free Software Foundation, Inc.
This file is part of the GNU C Library.
The GNU C Library is free software; you can redistribute it and/or
modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.
The GNU C Library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with the GNU C Library; if not, see
<http://www.gnu.org/licenses/>. */
#include <sysdep.h>
#include <init-arch.h>
/* Define multiple versions only for the definition in libc and for
the DSO. In static binaries we need strrchr before the initialization
happened. */
#if defined SHARED && !defined NOT_IN_libc
.text
ENTRY(strrchr)
.type strrchr, @gnu_indirect_function
cmpl $0, __cpu_features+KIND_OFFSET(%rip)
jne 1f
call __init_cpu_features
1: leaq __strrchr_sse2(%rip), %rax
testl $bit_Slow_SSE4_2, __cpu_features+CPUID_OFFSET+index_Slow_SSE4_2(%rip)
jnz 2f
testl $bit_SSE4_2, __cpu_features+CPUID_OFFSET+index_SSE4_2(%rip)
jz 2f
leaq __strrchr_sse42(%rip), %rax
ret
2: testl $bit_Slow_BSF, __cpu_features+FEATURE_OFFSET+index_Slow_BSF(%rip)
jz 3f
leaq __strrchr_sse2_no_bsf(%rip), %rax
3: ret
END(strrchr)
/*
This implementation uses SSE4 instructions to compare up to 16 bytes
at a time looking for the last occurrence of the character c in the
string s:
char *strrchr (const char *s, int c);
We use 0x4a:
_SIDD_SBYTE_OPS
| _SIDD_CMP_EQUAL_EACH
| _SIDD_MOST_SIGNIFICANT
on pcmpistri to compare xmm/mem128
0 1 2 3 4 5 6 7 8 9 A B C D E F
X X X X X X X X X X X X X X X X
against xmm
0 1 2 3 4 5 6 7 8 9 A B C D E F
C C C C C C C C C C C C C C C C
to find out if the first 16byte data element has a byte C and the
last offset. There are 4 cases:
1. The first 16byte data element has EOS and has the byte C at the
last offset X.
2. The first 16byte data element is valid and has the byte C at the
last offset X.
3. The first 16byte data element has EOS and doesn't have the byte C.
4. The first 16byte data element is valid and doesn't have the byte C.
Here is the table of ECX, CFlag, ZFlag and SFlag for 3 cases:
case ECX CFlag ZFlag SFlag
1 X 1 1 0
2 X 1 0 0
3 16 0 1 0
4 16 0 0 0
We exit from the loop for cases 1 and 3 with jz which branches
when ZFlag is 1. If CFlag == 1, ECX has the offset X for case 1. */
.section .text.sse4.2,"ax",@progbits
.align 16
.type __strrchr_sse42, @function
.globl __strrchr_sse42
.hidden __strrchr_sse42
__strrchr_sse42:
cfi_startproc
CALL_MCOUNT
testb %sil, %sil
je __strend_sse4
# ifdef __CHKP__
bndcl (%rdi), %bnd0
bndcu (%rdi), %bnd0
# endif
xor %eax,%eax /* RAX has the last occurrence of s. */
movd %esi, %xmm1
punpcklbw %xmm1, %xmm1
movl %edi, %esi
punpcklbw %xmm1, %xmm1
andl $15, %esi
pshufd $0, %xmm1, %xmm1
movq %rdi, %r8
je L(loop)
/* Handle unaligned string using psrldq. */
leaq L(psrldq_table)(%rip), %rdx
andq $-16, %r8
movslq (%rdx,%rsi,4),%r9
movdqa (%r8), %xmm0
addq %rdx, %r9
jmp *%r9
/* Handle unaligned string with offset 1 using psrldq. */
.p2align 4
L(psrldq_1):
psrldq $1, %xmm0
.p2align 4
L(unaligned_pcmpistri):
pcmpistri $0x4a, %xmm1, %xmm0
jnc L(unaligned_no_byte)
leaq (%rdi,%rcx), %rax
L(unaligned_no_byte):
/* Find the length of the unaligned string. */
pcmpistri $0x3a, %xmm0, %xmm0
movl $16, %edx
subl %esi, %edx
cmpl %ecx, %edx
/* Return RAX if the unaligned fragment to next 16B already
contain the NULL terminator. */
jg L(exit)
addq $16, %r8
# ifdef __CHKP__
bndcu (%r8), %bnd0
# endif
/* Loop start on aligned string. */
.p2align 4
L(loop):
pcmpistri $0x4a, (%r8), %xmm1
jbe L(match_or_eos)
addq $16, %r8
# ifdef __CHKP__
bndcu (%r8), %bnd0
# endif
jmp L(loop)
.p2align 4
L(match_or_eos):
je L(had_eos)
L(match_no_eos):
leaq (%r8,%rcx), %rax
addq $16, %r8
# ifdef __CHKP__
bndcu (%r8), %bnd0
# endif
jmp L(loop)
.p2align 4
L(had_eos):
jnc L(exit)
leaq (%r8,%rcx), %rax
# ifdef __CHKP__
bndcu (%rax), %bnd0
# endif
.p2align 4
L(exit):
ret
/* Handle unaligned string with offset 15 using psrldq. */
.p2align 4
L(psrldq_15):
psrldq $15, %xmm0
jmp L(unaligned_pcmpistri)
/* Handle unaligned string with offset 14 using psrldq. */
.p2align 4
L(psrldq_14):
psrldq $14, %xmm0
jmp L(unaligned_pcmpistri)
/* Handle unaligned string with offset 13 using psrldq. */
.p2align 4
L(psrldq_13):
psrldq $13, %xmm0
jmp L(unaligned_pcmpistri)
/* Handle unaligned string with offset 12 using psrldq. */
.p2align 4
L(psrldq_12):
psrldq $12, %xmm0
jmp L(unaligned_pcmpistri)
/* Handle unaligned string with offset 11 using psrldq. */
.p2align 4
L(psrldq_11):
psrldq $11, %xmm0
jmp L(unaligned_pcmpistri)
/* Handle unaligned string with offset 10 using psrldq. */
.p2align 4
L(psrldq_10):
psrldq $10, %xmm0
jmp L(unaligned_pcmpistri)
/* Handle unaligned string with offset 9 using psrldq. */
.p2align 4
L(psrldq_9):
psrldq $9, %xmm0
jmp L(unaligned_pcmpistri)
/* Handle unaligned string with offset 8 using psrldq. */
.p2align 4
L(psrldq_8):
psrldq $8, %xmm0
jmp L(unaligned_pcmpistri)
/* Handle unaligned string with offset 7 using psrldq. */
.p2align 4
L(psrldq_7):
psrldq $7, %xmm0
jmp L(unaligned_pcmpistri)
/* Handle unaligned string with offset 6 using psrldq. */
.p2align 4
L(psrldq_6):
psrldq $6, %xmm0
jmp L(unaligned_pcmpistri)
/* Handle unaligned string with offset 5 using psrldq. */
.p2align 4
L(psrldq_5):
psrldq $5, %xmm0
jmp L(unaligned_pcmpistri)
/* Handle unaligned string with offset 4 using psrldq. */
.p2align 4
L(psrldq_4):
psrldq $4, %xmm0
jmp L(unaligned_pcmpistri)
/* Handle unaligned string with offset 3 using psrldq. */
.p2align 4
L(psrldq_3):
psrldq $3, %xmm0
jmp L(unaligned_pcmpistri)
/* Handle unaligned string with offset 2 using psrldq. */
.p2align 4
L(psrldq_2):
psrldq $2, %xmm0
jmp L(unaligned_pcmpistri)
cfi_endproc
.size __strrchr_sse42, .-__strrchr_sse42
.section .rodata.sse4.2,"a",@progbits
.p2align 4
L(psrldq_table):
.int L(loop) - L(psrldq_table)
.int L(psrldq_1) - L(psrldq_table)
.int L(psrldq_2) - L(psrldq_table)
.int L(psrldq_3) - L(psrldq_table)
.int L(psrldq_4) - L(psrldq_table)
.int L(psrldq_5) - L(psrldq_table)
.int L(psrldq_6) - L(psrldq_table)
.int L(psrldq_7) - L(psrldq_table)
.int L(psrldq_8) - L(psrldq_table)
.int L(psrldq_9) - L(psrldq_table)
.int L(psrldq_10) - L(psrldq_table)
.int L(psrldq_11) - L(psrldq_table)
.int L(psrldq_12) - L(psrldq_table)
.int L(psrldq_13) - L(psrldq_table)
.int L(psrldq_14) - L(psrldq_table)
.int L(psrldq_15) - L(psrldq_table)
# undef ENTRY
# define ENTRY(name) \
.type __strrchr_sse2, @function; \
.align 16; \
.globl __strrchr_sse2; \
.hidden __strrchr_sse2; \
__strrchr_sse2: cfi_startproc; \
CALL_MCOUNT
# undef END
# define END(name) \
cfi_endproc; .size __strrchr_sse2, .-__strrchr_sse2
# undef libc_hidden_builtin_def
/* It doesn't make sense to send libc-internal strrchr calls through a PLT.
The speedup we get from using SSE4.2 instruction is likely eaten away
by the indirect call in the PLT. */
# define libc_hidden_builtin_def(name) \
.globl __GI_strrchr; __GI_strrchr = __strrchr_sse2
#endif
#include "../strrchr.S"
|
