<ins id="pjuwb"></ins>

<blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>

<noscript id="pjuwb"></noscript>

<sup id="pjuwb"><pre id="pjuwb"></pre></sup>

<dd id="pjuwb"></dd>

<abbr id="pjuwb"></abbr>

string

string

posts - 27, comments - 177, trackbacks - 0, articles - 0

C++博客 :: 首頁 :: 新隨筆 :: 聯系 :: 聚合

:: 管理

基于sse2的strstr函數

Posted on 2008-10-28 21:47 djx_zh 閱讀(2884) 評論(7) 編輯收藏引用

download the code
昨天實現了基于int類型的strstr函數，可以獲得1~2X左右的加速。今天按 lstrstr的流程實現了基于SSE2的STRSTR函數。可以得到2~4X左右的加速。

  1 char* lstrstrsse(char* text, char* pattern)
  2 {
  3     __m128i * sseiPtr = (__m128i *) text;
  4     unsigned char * chPtrAligned = (unsigned char*)text;
  5     __m128i sseiWord0 ;//= *sseiPtr ;
  6     __m128i sseiWord1 ;//= *sseiPtr ;
  7     __m128i sseiZero = _mm_set1_epi8(0);
  8     char chara = pattern[0];
  9     char charb = pattern[1];
10     register __m128i byte16a;
11     register __m128i byte16b;
12     char* bytePtr =text;
13     if(pattern ==NULL) return NULL;
14     if(pattern[0] == 0) return NULL;
15     if(pattern[1] == 0) return lstrchr(text,pattern[0]);
16     byte16a = _mm_set1_epi8(chara);
17     byte16b = _mm_set1_epi8(charb);
18 // process the unaligned bytes
19

20 // the aligned bytes
21 alignStart:
22     sseiWord0 = *sseiPtr;
23     sseiWord1 = *(sseiPtr+1);
24     while( haszeroByte(sseiWord0,sseiWord1,sseiZero) ==0)
25     {
26         unsigned int reta ;
27 searcha:
28         reta = hasByteC(sseiWord0,sseiWord1,  byte16a);
29         if(reta!=0 ) {
30             unsigned int retb ;
31 findouta:
32             retb = hasByteC(sseiWord0,sseiWord1,  byte16b);
33 findoutb:
34             if(((reta<<1) & retb)){
35                 // have ab
36                 int i=1;
37                 char * bytePtr0 = (char*) ( sseiPtr );
38                 int j;
39                 //printf("test::%0x,%d\n",reta ,bytePtr0 -text);
40                 bytePtr = (char*) ( sseiPtr );
41                 for(j =0;j<8;j++){
42                     if(reta & 0xff) {
43                         if(bytePtr0[0] == chara){
44                             i =1;
45                             bytePtr = bytePtr0 ;
46                             while((pattern[i] )&&(bytePtr[i] == pattern[i])) i++;
47                             if(pattern[i] == 0) return bytePtr;
48                         }
49                         if(bytePtr0[1] == chara){
50                             i =1;
51                             bytePtr = bytePtr0 + 1;
52                             while((pattern[i] )&&(bytePtr[i] == pattern[i])) i++;
53                             if(pattern[i] == 0) return bytePtr;
54                         }
55                         if(bytePtr0[2] == chara){
56                             i =1;
57                             bytePtr = bytePtr0 + 2;
58                             while((pattern[i] )&&(bytePtr[i] == pattern[i])) i++;
59                             if(pattern[i] == 0) return bytePtr;
60                         }
61                         if(bytePtr0[3] == chara){
62                             i =1;
63                             bytePtr = bytePtr0 + 3;
64                             while((pattern[i] )&&(bytePtr[i] == pattern[i])) i++;
65                             if(pattern[i] == 0) return bytePtr;
66                         }
67                     }
68                     reta = reta >> 4;
69                     bytePtr0 += 4;
70                 }
71             }
72             // search b
73             sseiPtr += 2;
74             sseiWord0 = *sseiPtr;
75             sseiWord1 = *(sseiPtr+1);
76
77             while( haszeroByte(sseiWord0,sseiWord1,sseiZero) ==0){
78                 retb = hasByteC(sseiWord0,sseiWord1,  byte16b);
79                 if(retb !=0){
80                     // findout b
81                     if((*((char*) sseiPtr)) == charb){
82                         //b000
83                         char * bytePtr = ((char*) ( sseiPtr )) -1;
84                         if(bytePtr[0] == chara){
85                             int i=1;
86                             while((pattern[i] )&&(bytePtr[i] == pattern[i])) i++;
87                             if(pattern[i] == 0) return bytePtr;
88                             if(bytePtr[i] == 0) return NULL;
89                         }
90
91                     }
92                     reta = hasByteC(sseiWord0,sseiWord1,  byte16a);
93                     if(reta !=0)
94                         goto findoutb;
95                     else{
96                         goto nextWord;
97                     }
98                 }
99                 sseiPtr += 2;
100                 sseiWord0 = *sseiPtr;
101                 sseiWord1 = *(sseiPtr+1);
102             }
103             // search  from (char*)sseiPtr
104             char * bytePtr = ((char*) ( sseiPtr )) -1;
105             if(bytePtr[0] == chara){
106                 int i=1;
107                 while((pattern[i] )&&(bytePtr[i] == pattern[i])) i++;
108                 if(pattern[i] == 0) return bytePtr;
109             }
110
111             goto prePareForEnd;
112         }
113 nextWord:
114         sseiPtr += 2;
115         sseiWord0 = *sseiPtr;
116         sseiWord1 = *(sseiPtr+1);
117     }
118 prePareForEnd:
119     {
120         unsigned int reta;
121         unsigned int retb;
122         reta =hasByteC(sseiWord0,sseiWord1,  byte16a);
123         retb =hasByteC(sseiWord0,sseiWord1,  byte16b);
124         if(((reta<<1) & retb)){
125             bytePtr = (char*)sseiPtr;
126             while(*bytePtr){
127                 if(*bytePtr == chara) {
128                     int i=1;
129                     while((pattern[i] )&&(bytePtr[i] == pattern[i])) i++;
130                     if(pattern[i] == 0) return bytePtr;
131                     if(bytePtr[i] == 0) return NULL;
132
133                 }
134                 bytePtr++;
135             }
136         }
137     }
138     return NULL;
139 }
140

Feedback

# re: 基于sse2的strstr函數回復 更多評論

2008-10-30 00:53 by 肥仔

超過了C的strstr?

# re: 基于sse2的strstr函數回復 更多評論

2008-10-30 09:33 by djxzh

@肥仔
就目前的測試結果，是這樣。還沒有測試最壞情況下會是什么結果。

# re: 基于sse2的strstr函數[未登錄] 回復 更多評論

2008-10-30 10:35 by megax

做一個從后面開始查找的試試？

# re: 基于sse2的strstr函數回復 更多評論

2008-10-30 10:41 by djxzh

@megax
你是說BM之類的算法嗎？那些算法需要對模式串預處理。

# re: 基于sse2的strstr函數回復 更多評論

2008-10-30 10:54 by vczh

用了SSE的指令集就可以同時計算一小部分內容了。

# re: 基于sse2的strstr函數[未登錄] 回復 更多評論

2008-10-31 12:47 by megax

不是，我說的是從一個字符串后面開始查找想要查找的內容。不是說具體的算法

# re: 基于sse2的strstr函數 回復 更多評論

2008-11-01 11:02 by 金山詞霸2008

沒想到strstr函數的設計還這么復雜

刷新評論列表

只有注冊用戶登錄后才能發表評論。
【推薦】100%開源！大型工業跨平臺軟件C++源碼提供，建模，組態！



網站導航: 博客園 IT新聞 BlogJava 博問 Chat2DB 管理

狠狠色丁香久久综合五月| 亚洲人成精品久久久久| 久久国产乱子伦免费精品| 国产精品乱码久久久久久软件| 日韩一区二区久久久久久| 国产精品久久久久AV福利动漫| 亚洲熟妇无码另类久久久| 中文成人无码精品久久久不卡| 很黄很污的网站久久mimi色| 久久精品中文字幕久久| 久久成人国产精品二三区| 精品国产一区二区三区久久| 久久免费精品一区二区| 精品久久久久久久久久中文字幕| 精品人妻伦九区久久AAA片69| 欧美性大战久久久久久| 久久天天躁狠狠躁夜夜2020一| 亚洲国产一成人久久精品| 国产精品久久自在自线观看| 国产精品九九久久免费视频 | 亚洲美日韩Av中文字幕无码久久久妻妇| 久久久噜噜噜久久中文字幕色伊伊| 99久久精品国产综合一区| 久久久久97国产精华液好用吗| 伊人久久大香线蕉AV一区二区| 亚洲欧美日韩久久精品第一区| 国产亚洲欧美精品久久久| 香蕉久久夜色精品国产小说| 日韩影院久久| AV无码久久久久不卡网站下载| 美女写真久久影院| 精品久久久久久中文字幕大豆网| 国产精品久久久久影视不卡| 亚洲国产精品综合久久一线| 国产亚洲精久久久久久无码| 美女久久久久久| 国产精品久久久久久福利漫画| 伊人色综合久久天天网| 亚洲伊人久久大香线蕉苏妲己| 2020国产成人久久精品| 7国产欧美日韩综合天堂中文久久久久 |