大数运算基选择测试

无心人

两个二进制数字， bit位长度是28*30*32
分别用28，30，32个bit为一个处理单位
做加法和乘法
测试在不同基选择下的效果
===============================
这么写如果不妥， 你们可提出意见
咱一起商量测试方案

liangbch

应该加上一些限制，如不得使用高级算法，如分治法，toom-cook 算法等。

无心人

28*30*32=26880
除非有贼快的FFT处理器，否则没人敢用FFT
同样Toom-n在这么小，优势也不大
就那个K算法还能加速
如果位数太大
可以考虑做多次测试
28*32， 28*30，30*32各一次

无心人

//我先提交我的函数，由旧代码修改成
//终于接近完美了，不过还有几条指令似乎可以省掉
//怎么省想不出来 ：）
void  lMul_yaos(uint32  *result, uint32  *left, uint32  *right, uint32 sLeft, uint32 sRight)
{
__asm          {
             mov edx, sLeft
             mov  esi, dword ptr [left]
             mov  edi, dword ptr [right]
             mov  ebx, dword ptr [result]
mbinmul2:
             mov eax, ebx
             pxor mm0, mm0
             mov ecx, sRight                        
             movd mm1, dword ptr [esi]
             mov ebp, edi
mbinmul3:
             movd mm2, dword ptr [edi]
             lea edi, [edi+4]
             movd mm3, dword ptr [ebx]                        
             pmuludq mm2, mm1
             paddq mm0, mm3
             paddq mm0, mm2
             movd dword ptr [ebx], mm0
             psrlq mm0, 32
             lea ebx, [ebx+4]
             sub ecx, 1
             jnz mbinmul3
             mov edi, ebp
             movd dword ptr [ebx], mm0
             mov ebx, eax
             lea esi, [esi+4]
             lea ebx, [ebx+4]
             sub edx, 1
             jnz mbinmul2
             emms
             }
}

gxqcn

建议将测试的bit长度缩小为：7*15*32 = 3360
在 28*30*32bits 时早该用更好的算法了。

无心人

好
同意你的意见
我只管提交我的标准32位乘法
具体测试你们来
：）

无心人

我的代码必须配合一个清零函数
等我写来
：）

liangbch

我打算 完全用C 编写一个程序，看看你的那个使用 SSE2 指令编写的程序能比我的C版的程序快多少。

无心人

有限的很吧
假设我的时间是100
C写的应该在120-160之间

liangbch

原帖由 无心人 于 2008-4-2 11:24 发表


有限的很吧
假设我的时间是100
C写的应该在120-160之间

  由那个128bit * 128 bit 帖子，得出一个结论，使用MMX寄存器的版本(96楼)最高可比使用ALU指令的版本(102楼)快2.33倍(在我的电脑上,前者运行时间为400.665,或者运行时间为935.273),也就是说 使用 MMX的版本 加速比 可达 2.33倍.
  退一步,完全使用c的版本性能可以达到102楼的性能, 我的C版程序运行时间应为 你的4#楼那个版本的2.33倍,你说预计可达你的你的版本的120-160%, 这等于承认了基为$2^30$的版本速度为基为$2^32$版本的 146%- 194%.