使用CPU 专用指令提高float运算速度的一个小例子

system888net

若有一批float要进行y= sqrt(x1*x1+x2*x2)+0.5 的计算
在32bit的CPU上用常规模式进行计算则明显耗时.
而用INTEL CPU自己专用的SSE指令进行运算,则有明显的优势:
windows2000+ VC

1. 
   
2. int i;
   
3.              float x1_float_data_block[]={...}; // x1数据块
   
4.              float x2_float_data_block[]={...}; // x2数据块
   
5.              float y_float_data_block[sizeof(x1_float_data_block)/sizeof(float)];    // y数据块,结果
   
6.              int num_float=sizeof(x1_float_data_block)/sizeof(float);                // 数据块中float的数量
   
7. 
   
8.              int count = num_float/4;
   
9.              float f = 0.5f;
   
10.              _asm    {
    
11.                          movss   xmm2, f         
    
12.                          shufps  xmm2, xmm2, 0   
    
13.                          mov         esi, x1_float_data_block               
    
14.                          mov         edx, x2_float_data_block               
    
15.                          mov         edi, y_float_data_block               
    
16.                          mov         ecx, count                             
    
17.                  sse_loop:        
    
18.                          movaps      xmm0, [esi]
    
19.                          mulps       xmm0, xmm0  
    
20.                          movaps      xmm1, [edx]
    
21.                          mulps       xmm1, xmm1  
    
22.                          addps       xmm0, xmm1  
    
23.                          sqrtps      xmm0, xmm0  
    
24.                          addps       xmm0, xmm2  
    
25.                          movaps      [edi],xmm0  
    
26.                          add         esi, 16     
    
27.                          add         edx, 16     
    
28.                          add         edi, 16     
    
29.                          dec         ecx         
    
30.                          jnz         sse_loop                    
    
31.                      }
    
32.                for(i=0;i<num_float;i++)  
    
33.                          printf("...",y_float_data_block[i],...
    
34.                .....
    

复制代码

[ 本帖最后由 system888net 于 2008-5-18 10:51 编辑 ]

system888net

这条printf("...",y_float_data_block,...
在上面edit了几次总是显示不出来[ i ],不知为何

[ 本帖最后由 system888net 于 2008-4-27 20:23 编辑 ]

chzCPU

运算速度能提高一点还是许多?

system888net

原帖由 jamesr 于 2008-4-28 19:16 发表
使用 代把码括起来。
另外，应该说明这段代码的运行环境。这种内联汇编不是gcc格式的。

谢谢提醒,已加了code 和/code,平台 windows2000+vc

[ 本帖最后由 system888net 于 2008-5-18 10:51 编辑 ]

chuhongze

速度是明显快了,顶...

chuhongze

此程序调用的是FPU内部的指令