代码优化-之-优化除法.docxVIP

出很多人的预料，它是一条很慢的操作，整数和浮点的除法都慢；我测试的英特尔P5 赛扬精度数除法指令 出很多人的预料，它是一条很慢的操作，整数和浮点的除法都慢；我测试的英特尔P5 赛扬 精度数除法指令 DIVSS 36 个 CPU 周期) (x86 求余运算和除法运算是用同一条 CPU 指 代码优化－之－优化除法 HouSisong@GM tag:代码优化,除法,牛顿迭代,减法代替除法,除法优化 说明：文章中的很多数据可能在不同的CPU 或不同的系统环境下有不同的结果，数据仅 供参考x86 供参考 x86 系列的 CPU 对于位运算、加、减等基本指令都能在1 个 CPU 周期内完成(现在的 CPU 还能乱序执行，从而使指令的平均 CPU 周期更小)；现在的 CPU,做乘法也是很快的(需要几个 CPU 周期，每个周期可能启动一个新的乘指令(x87))，但作为基本指令的除法却超 CP CPU 浮点数的除法差不多是 37 个 CPU 周期，整数的除法是 80 个 CPU 周期，AMD2200+ 浮点数的除法差不多是 21 个 CPU 周期，整数的除法是 40 个 CPU 周期。(改变 FPU 运算精度对于除法无效) (SSE 指令集的低路单精度数除法指令 DIVPS 18 个 CPU 周期,四路单 令实现的； 据说，很多 令实现的； 据说，很多 CPU 的整数除法都是用数学协处理器的浮点除法器完成的；有一个 推论就是，浮点除法和整数除法不能并行执行. (ps:intel 的 p4 imul 指令可能有 14 周期(或 15-18)的延迟才能得到结果) 本文将给出一些除法的优化方法或替代算法 本文将给出一些除法的优化方法或替代算法 (警告:某些替代算法并不能保证完全等价!) 1. 1.尽量少用除法 比如： if (x/yz) ... 改成： if ( ((y0)(xy*z))||((y0)(xy*z)) ) ... ( (少用求余) 比如： ++index; if (index=count) index=index % count; //assert(indexcount*2);改成： count; //assert(indexcount*2); 改成： ++index; if (index=count) index=index - count; 2. 2.用减法代替除法 如果知道被除数是除数的很小的倍数，那么可以用减法来代替除法比如： uint32 x=200; uint32 y=70; uint32 y=70; uint32 z=x/y;改成： uint32 z=x/y; 改成： uint z=0; while (x=y) { x-=y; ++z;} x-=y; ++z; } 一个用减法和移位完成的除法 一个用减法和移位完成的除法 (如果你没有除法指令可用:) uint32 div(uint64 u,uint32 z) //return u/z {uint32 { uint32 x=(uint32)(u32); uint32 y=(uint32)u; //y //y 保存商 x 保存余数 for (int for (int i=0;i32;++i) { uint32 t=((int32)x) 31; x=(x1)|(y31); y=y1; y=y1; if ((x|t)=z) if ((x|t)=z) { x-=z; ++y;} ++y; } } return y; }(该函数经过了测试； } (该函数经过了测试；z==0 需要自己处理；对于有符号除法，可以用取绝对值的方法(当 然也不是轻松就能 写出完全等价的有符号除法的:)； 如果不需s 的 64bit 长度，仅需要32bit,那么可以化简这个函数，但改进不多) 3. 3.用移位代替除法 (很多编译器能自动做好这个优化) 要求除数是 2 的次方的常量； (同理：对于某些应用，可以优先选取这样的数来做除数) 比如： uint32 x=213432575; uint32 y=x/8; 改成： y=x3; 对于有符号的整数； 对于有符号的整数； 比如： int32 x=213432575; □ int32 y=x/8; 改成： if 改成： if (x=0) y=x3; else y=(x+(13-1))3; 4.合并除法 4.合并除法 (替代方法不等价，很多编译器都不会帮你做这种优化) 适用于不能用其它方法避免除法的时候； 比如： double x=a/b/c; 改成： double x=a/(b*c); 比如： double x=a/b+c/b; 改成： 改成： double x=(a+c)/b; 比如： double x=a/b;