Re: [问题] glibc 2.7 strlen.c

楼主: johnjohnlin (嗯?)   2015-04-28 23:24:06
试着练习一下
37 /* Handle the first few characters by reading one character at a time.
38 Do this until CHAR_PTR is aligned on a longword boundary. */
39 for (char_ptr = str; ((unsigned long int) char_ptr
40 & (sizeof (longword) - 1)) != 0;
41 ++char_ptr)
42 if (*char_ptr == '\0')
43 return char_ptr - str;
这边注解写的很清楚
把 char pointer 对齐到 long* 的位置
这样之后我们一次就可以比一个 long
59 magic_bits = 0x7efefeffL;
60 himagic = 0x80808080L;
61 lomagic = 0x01010101L;
62 if (sizeof (longword) > 4)
63 {
66 magic_bits = ((0x7efefefeL << 16) << 16) | 0xfefefeffL;
67 himagic = ((himagic << 16) << 16) | himagic;
68 lomagic = ((lomagic << 16) << 16) | lomagic;
69 }
70 if (sizeof (longword) > 8)
71 abort ();
magic number 先跳过
不过方便起见之后假设 long 是 4B
76 for (;;)
77 {
107 longword = *longword_ptr++;
这个跟你的写法差不多
不过他现在就是要一次比 4B
108
109 if (
122 ((longword - 0x01010101) & 0x80808080)
124 != 0)
125 {
整个程式里面神奇的只有这个
我把 magic number 展开了
如果 4B 里面有 0
那么每个 byte 减掉 1 的时候就会退位
导致该 byte 前面的 byte 的 MSB 是 1
128
129 const char *cp = (const char *) (longword_ptr - 1);
130
131 if (cp[0] == 0)
132 return cp - str;
133 if (cp[1] == 0)
134 return cp - str + 1;
(下略)
承上
会进来这边就是他认为这 4B 里面有可能 0x00
这边就是把 loop 展开
作者: EdisonX (卡卡兽)   2015-04-28 23:30:00
vc 也是用这种作法 , 只是它已包成了 strlen.asm ,所以..
作者: OPIV (Monitor)   2015-04-28 23:32:00
抱歉...J大 我不太了解loop展开的意思意思是,把原本的逐位检查改成if()if()if()吗?
作者: LiloHuang (十年一刻)   2015-04-28 23:34:00
楼主: johnjohnlin (嗯?)   2015-04-28 23:43:00
loop unrolling
作者: PkmX (阿猫)   2015-04-29 01:31:00
glibc的实作真变态 比我用AVX2一次比对32 bytes还要快...https://gist.github.com/PkmX/ea6ee1e0b56187006e80
楼主: johnjohnlin (嗯?)   2015-04-29 12:45:00
测试了一下,发现真的会快一点点,到底为什么啊
作者: PkmX (阿猫)   2015-04-29 13:08:00
可能glibc的版本我的cpu的micro-architecture比较喜欢吧XDrz不过基本上差距都在5%以内 所以可接受
作者: TobyH4cker (Toby (我要当好人))   2015-04-29 16:02:00
应该是有对齐比较符合CPU的运作方式这部分只要看编译后的asm就能比较了
作者: PkmX (阿猫)   2015-04-29 18:58:00
我的AVX2版本也有先对齐啊XDrz
作者: LiloHuang (十年一刻)   2015-04-29 20:41:00
因为 glibc 的实作里头有 SSE2 的版本也有 SSE4详见 http://goo.gl/kzQ9sX剩下就是得分析此情况下 AVX2 跟 SSE4 为啥会输了 :P这篇比较的表格也许可以看看 http://goo.gl/G04iPI说不定你的 glibc 刚好有 AVX2 实作也不一定 XD用 gdb 下个断点 step into 进去 strlen 看看应该可知
楼主: johnjohnlin (嗯?)   2015-04-29 22:39:00
我比较的对象是上一篇那个纯 C 的实作另外我有点好奇这样不会出现违规存取的问题吗?我电脑里面的是 SSE2 实作,速度跟纯 C 很接近
作者: PkmX (阿猫)   2015-04-30 08:25:00
基本上你的memory access都是align在4或8 bytes上不会有跨越page的问题 所以在x86上面ok

Links booklink

Contact Us: admin [ a t ] ucptt.com