快速无符号整数开方

最新推荐文章于 2024-01-29 21:49:01 发布

原创

最新推荐文章于 2024-01-29 21:49:01 发布 · 3k 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#测试

这篇博客展示了几个不同的无符号整数快速开方的实现，包括fast_sqrt1到fast_sqrt4以及iSqrt_FPU2_yaos和iSqrt_FPU1_lbc。代码中使用了FPU浮点运算进行优化，并通过测试用例验证了它们的性能和边界情况。

看了下liangbch 兄打包的Uint32sqrt测试程序，发现进一步改进FPU式的很难，不过我还是勉强做了些细微的改动，给出四个版本加上iSqrt_FPU2_yaos，iSqrt_FPU1_lbc。

代码： #include <stdio.h> #include <time.h> typedef unsigned int DWORD; double b32[] = {0.0,4294967296.0}; __declspec(naked) DWORD __fastcall iSqrt_FPU2_yaos(DWORD n) { __asm { push ecx mov eax, ecx and eax, 0x80000000 shr eax, 31 fld qword ptr [b32 + eax * 8] fild dword ptr [esp] faddp st(1), st fsqrt sub esp, 8 fstp qword ptr [esp] mov edx, dword ptr [esp + 4] mov eax, edx and edx,0x7ff00000 and eax,0xfffff shr edx, 20 or eax, 0x100000 xchg ecx, edx sub ecx, 1043 neg ecx shr eax, cl xchg edx, ecx add esp, 12 or ecx, ecx cmove eax, ecx ret } } _declspec(naked) DWORD fast_sqrt1(DWORD x) { _asm { sub esp,4 mov dword ptr [esp+12],0 fild qword ptr [esp+8]