c - - assembly sse - 段错误 - 在另一个矩阵中移动一个矩阵行

Question

我有两个矩阵 nxn，A 和 B。

A=[1 2 3 4 ; 5 6 7 8 ; 9 10 11 12 ; 13 14 15 16]

B=[1 1 1 1 ; 1 1 1 1 ; 1 1 1 1 ; 1 1 1 1 ]

这两个矩阵在 C 中分配为 float*。在 main.c 中，我在程序集 'exportRow(A,B,n);' 中调用了一个过程

在 assembly 中，

    ...
    mov     eax,[ebp+a]
    mov     ebx,[ebp+b]
    mov     [xa],eax
    mov     [xb],ebx
    printps [xa],1
    printps [xb],1
    movaps  xmm0,[xa]
    movaps  [xb],xmm0        ;casuses segmentation fault
    ...

printps 是一个宏，它只用于打印 xa 和 xb 的一个位置内存(它们是双字)。我想使用 xa 和 xb 作为矩阵的指针。

输出:

[1,2,3,4] 
[1,1,1,1]
SEGMENTATION FAULT

为什么最后一条指令导致段错误？哪里有问题？如果我在最后一条指令中将 [xb] 替换为 [xa](仅用于测试)，我没有任何错误。

--编辑---

如果我修改 [xb]，修改应该不会影响 B？

不过，稍微修改一下:(实际上xa和xb是resd，并按32字节对齐。)

 movaps  xmm0,[xa+16]
 movaps [xb],xmm0
 printps [xb],1

段错误....

---编辑----

这是一个小代码示例:

int main(){
 int m=4;
 int n=4;
 float* A=malloc(sizeof(float)*m*n);
 float* B=malloc(sizeof(float)*m*n);
 for(int i=0;i<m*n;i++){
     A[i]=i;
     B[i]=1;
 }
 printMatrix(A,m,n);
 printMatrix(B,m,n);
 exportRow(A,B,n); //procedure assembly 
 printMatrix(B,m,n);

汇编代码:

section .data
   A equ 8
   B equ 12
   n equ 16
section .bss
   alignb 32
   xa resd 1
   alignb 32
   xb resd 1
section .text
global exportRow

exportRow: push ebp
           mov ebp,esp
           push ebx
           push esi
           push edi

           mov     eax,[ebp+a]
           mov     ebx,[ebp+b]
           mov     [xa],eax
           mov     [xb],ebx
           printps [xb],1
           printps [xa+16],1
           movaps  xmm0,[xa+16]
           movaps  [xb],xmm0        
           printps [xb],1

           pop edi
           pop esi
           pop ebx
           mov esp,ebp
           pop ebp
           ret

输出:

A=[1 2 3 4 ; 5 6 7 8 ; 9 10 11 12 ; 13 14 15 16]

B=[1 1 1 1 ; 1 1 1 1 ; 1 1 1 1 ; 1 1 1 1 ]

[1,1,1,1]
SEGMENTATION FAULT

大小写对齐 b 32:如果我得到 [xa+16] ，我有段错误。

大小写对齐b 16:如果我得到 [xa+16]，我有 [1 1 1 1]，为什么？我以为结果是 [5 6 7 8]。

还有一个问题(我认为这两个问题是相关的):使用帖子顶部的代码:

mov     eax,[ebp+a]
mov     ebx,[ebp+b]
mov     [xa],eax
mov     [xb],ebx
printps [xa],1
printps [xb],1
movaps  xmm0,[xa]
movaps  [xb],xmm0 
printps [xb],1

输出:

A=[1 2 3 4 ; 5 6 7 8 ; 9 10 11 12 ; 13 14 15 16]

B=[1 1 1 1 ; 1 1 1 1 ; 1 1 1 1 ; 1 1 1 1 ]

[ 0.000000 1.000000 2.000000 3.000000 ]
[ 1.000000 1.000000 1.000000 1.000000 ]
[ 0.000000 1.000000 2.000000 3.000000 ]

B=[1.000000 1.000000    1.000000    1.000000    
1.000000    1.000000    1.000000    1.000000    
1.000000    1.000000    1.000000    1.000000    
1.000000    1.000000    1.000000    1.000000]

我修改了xb，但是矩阵B完好无损。为什么？

Answer 1

xb 没有 16 字节对齐的内存地址。如果您检查十六进制的内存地址，则必须始终以 0 结尾。如果不是，则需要使用 movups 而不是 movaps。

但是，简单地以 64 字节对齐几乎总是更好。

在 Windows 结账时:_aligned_malloc

编辑:所以看看您的 .bss 部分，您似乎并没有按照您认为的那样去做。您正在将存储在“EAX”中的指针复制到由单字节“xa”指向的地址。这不可能是正确的。

您需要确保您的原始分配是 16 字节对齐的。从那里，您可以直接通过 movaps 将 A 指向 B 的矩阵的一行直接移动到 xmm0 寄存器，然后从该寄存器移动另一个 movaps!

编辑 2:总而言之，xa 和 xb 毫无意义。

movaps  xmm0,[eax+16]
movaps  [ebx],xmm0 

虽然 malloc 分配的内存必须是 16 字节对齐的......