求随机产生0-2^32-1内整数算法

无心人

晕， 那样Call似乎不行
int rand(void)
汇编如何call?

无心人

要不手工写个小随机算法插入里面做位随机

无心人

找到个比较复杂的
给汇编化加位随机化


DWORD BitX[5];

void BitRandomInit(DWORD seed)
{
  DWORD i;
  DWORD s = seed;
  // make random numbers and put them into the buffer
  for (i = 0; i < 5; i++) {
    s = s * 29943829 - 1;
    BitX[ i ] = s;
  }
}

DWORD BitRandom(void)
{
        __asm
        {
        mov ecx, 31
        mov eax, 0
loop1:
        mov  ebx, 2111111111
    movd xmm5, ebx        
        movd xmm1, dword ptr [BitX + 12]
        pmuludq xmm5, xmm1
        mov  edx, 1492
        movd xmm6, edx
        movd xmm2, dword ptr [BitX + 8]
        pmuludq xmm6, xmm2
        paddq xmm5, xmm6
        movd dword ptr [BitX + 12], xmm2
        mov  esi, 1776
        movd xmm7, esi
        movd xmm3, dword ptr [BitX + 4]
        pmuludq xmm7, xmm3
        paddq xmm5, xmm7
        movd dword ptr [BitX + 8], xmm3
        mov  edi, 5115
    movd xmm0, edi
        movd xmm4, dword ptr [BitX]
        pmuludq xmm0, xmm4
        paddq xmm5, xmm0
        movd dword ptr [BitX + 4], xmm4
        movd xmm0, dword ptr [BitX + 16]
    paddq xmm0, xmm5
        movd edx, xmm0
        psrlq xmm0, 32
        movd dword ptr [BitX + 16], xmm0
        mov dword ptr [BitX], edx
        and  edx, 1
        shl edx, cl
        add eax, edx
        sub ecx, 1
        jnc loop1
        }
}

mathe

这样应该会好一些。不过这个要看计算随机数部分数据占用整个计算的比例是多少。通常，如果前面部分计算占用比列过大，很难衡量后面部分代码的时间，因为我们测量到的大部分时间是前面计算的时间。
还有，同样一段代码，最好多运行几次，看看运行时间的稳定性如何，如果多次运行发现时间变化相对比较大（比如已经超出两种不同算法得到的时间差），那么测量的意义也不是很大。

无心人

终于通过调试了
10万的时间是109ms

无心人

封存
处理完128位乘法如果有兴趣可能写个BBS位算法

无心人

1. 
   
2. DWORD BitX[5];
   
3. void BitRandomInit(DWORD seed)
   
4. {
   
5.   DWORD i;
   
6.   DWORD s = seed;
   
7.   // make random numbers and put them into the buffer
   
8.   for (i = 0; i < 5; i++)
   
9.   {
   
10.     s = s * 29943829 - 1;
    
11.     BitX[ i ] = s;
    
12.   }
    
13. }
    
14. 
    
15. 
    
16. DWORD BitRandom(void)
    
17. {
    
18.    __asm
    
19.   {
    
20.   mov ecx, 31
    
21.   mov eax, 0
    
22. loop1:
    
23.   mov ebx, 2111111111
    
24.   movd xmm5, ebx
    
25.   movd xmm1, dword ptr [BitX + 12] //BitX[3]
    
26.   pmuludq xmm5, xmm1
    
27.   mov edx, 1492
    
28.   movd xmm6, edx
    
29.   movd xmm2, dword ptr [BitX + 8] //BitX[2]
    
30.   pmuludq xmm6, xmm2
    
31.   paddq xmm5, xmm6
    
32.   movd dword ptr [BitX + 12], xmm2 //BitX[3]
    
33.   mov esi, 1776
    
34.   movd xmm7, esi
    
35.   movd xmm3, dword ptr [BitX + 4] //BitX[1]
    
36.   pmuludq xmm7, xmm3
    
37.   paddq xmm5, xmm7
    
38.   movd dword ptr [BitX + 8], xmm3 //BitX[2]
    
39.   mov edi, 5115
    
40.   movd xmm0, edi
    
41.   movd xmm4, dword ptr [BitX]    //BitX[0]
    
42.   pmuludq xmm0, xmm4
    
43.   paddq xmm5, xmm0
    
44.   movd dword ptr [BitX + 4], xmm4 //BitX[1]
    
45.   movd xmm0, dword ptr [BitX + 16] //BitX[4]
    
46.   paddq xmm0, xmm5
    
47.   movd edx, xmm0
    
48.   psrlq xmm0, 32
    
49.   movd dword ptr [BitX + 16], xmm0 //BitX[4]
    
50.   mov dword ptr [BitX], edx        //BitX[0]
    
51.   and edx, 1
    
52.   shl edx, cl
    
53.   add eax, edx
    
54.   sub ecx, 1
    
55.   jnc loop1
    
56.   }
    
57. }

复制代码

这个代码很有价值，重整理下