内嵌汇编应用sse计算矩阵的问题 请教！！

zhangyiviva

我想利用sse指令计算一个4*4的float矩阵，然后写了一段内嵌汇编的c程序，结果编译没问题，运行时出了问题。

程序在输出“hello”之后就卡在那里不动了，请教各位是什么问题啊～～

1. #include <stdio.h>
   
2. 
   
3. 
   
4.     float a[4][4] __attribute__ ((aligned (16))) = {
   
5.         {1.0,2.0,3.0,4.0},
   
6.         {1.0,2.0,3.0,4.0},
   
7.         {1.0,2.0,3.0,4.0},
   
8.         {1.0,2.0,3.0,4.0}
   
9.     };
   
10. 
    
11.     float b[4][4] __attribute__ ((aligned (16))) = {
    
12.         {1.0,2.0,3.0,4.0},
    
13.         {1.0,2.0,3.0,4.0},
    
14.         {1.0,2.0,3.0,4.0},
    
15.         {1.0,2.0,3.0,4.0}
    
16.     };
    
17. 
    
18. 
    
19.     float c[4][4] __attribute__ ((aligned (16)));
    
20. 
    
21. void matMul(float result[4][4],float left[4][4],float right[4][4])
    
22. {
    
23.         printf("hello\n");
    
24.         __asm__ __volatile__
    
25.         (
    
26.                 "movaps (%2),%%xmm4        \n\t"
    
27.                 "movaps 16(%2),%%xmm5        \n\t"
    
28.                 "movaps 32(%2),%%xmm6        \n\t"
    
29.                 "movaps 48(%2),%%xmm7        \n\t"
    
30.                 "movl $0,%%edi                \n\t"
    
31.                 "movl $4,%%ecx                \n\t"
    
32. "lab:\n\t"        "movups (%1,%2,1),%%xmm0        \n\t"
    
33.                 "movaps %%xmm0,%%xmm1        \n\t"
    
34.                 "movaps %%xmm0,%%xmm2        \n\t"
    
35.                 "movaps %%xmm0,%%xmm3        \n\t"
    
36.                 "shufps $0x00,%%xmm0,%%xmm0        \n\t"
    
37.                 "shufps $0x55,%%xmm1,%%xmm1        \n\t"
    
38.                 "shufps $0xAA,%%xmm2,%%xmm2        \n\t"
    
39.                 "shufps $0xFF,%%xmm3,%%xmm3        \n\t"
    
40.                 "mulps %%xmm4,%%xmm0        \n\t"
    
41.                 "mulps %%xmm5,%%xmm1        \n\t"
    
42.                 "mulps %%xmm6,%%xmm2        \n\t"
    
43.                 "mulps %%xmm7,%%xmm3        \n\t"
    
44.                 "addps %%xmm0,%%xmm2        \n\t"
    
45.                 "addps %%xmm1,%%xmm3        \n\t"
    
46.                 "addps %%xmm2,%%xmm3        \n\t"
    
47.                 "movups %%xmm3,(%0,%2,1)        \n\t"
    
48.                 "addl $16,%%edi                \n\t"
    
49.                 "loop lab                \n\t"
    
50.                 :"=d"(result)
    
51.                 :"a"(left),"D"(right)
    
52.                 :"%ecx"
    
53.         );
    
54. }
    
55. 
    
56. int main()
    
57. {
    
58.     int i,j;
    
59. 
    
60.     matMul(c,a,b);
    
61.    
    
62.     for(i = 0; i<4; i++)
    
63.     {
    
64.         for(j = 0; j<4; j++)
    
65.             printf("%lf ",c[i][j]);
    
66.         printf("\n");
    
67.     }
    
68. }

复制代码

mirnshi

loop lab 到啥年月是个头呀？

zhangyiviva

回复 2# mirnshi

谢谢你，我是汇编初学者～～

我想用sse算完结果就跳出来，应该怎么改？


   

zhangyiviva

mirnshi 发表于 2012-05-31 12:57
loop lab 到啥年月是个头呀？

每次loop都要递减ecx不是么？

开始设置了4，应该循环四次就退出了啊。。。

mirnshi

zhangyiviva 发表于 2012-05-31 14:01
每次loop都要递减ecx不是么？

开始设置了4，应该循环四次就退出了啊。。。

晕，你总得跳出循环吧。

zhangyiviva

mirnshi 发表于 2012-05-31 14:02
晕，你总得跳出循环吧。

不是，自己跳出的么？

那应该怎么跳出啊，多谢了！

zhangyiviva

回复 5# mirnshi

懂了，谢谢你～～



   

zhangyiviva

回复 5# mirnshi

现在变成段错误了。。。