使用 SSE 错误 __m128 到 *float 转换的矩阵乘法？

Question

我正在尝试使用 SSE Intrinsic 对矩阵乘法进行编程。我不确定我的代码是否正确，也无法编译它，因为我收到错误：

Error   1   error C2440: 'type cast' : cannot convert from 'float' to '__m128 *

有人可以仔细检查我的程序以确保我的矩阵乘法正确吗？另请注意，这是针对方阵的。

这是我的代码。

void Intrinsics (float * matrix_a, float * matrix_b, float * matrix_result, const int num_row, const int num_col) {
    __declspec(align(16)) float * a = matrix_a;
    __declspec(align(16)) float * b = matrix_b;
    __declspec(align(16)) float * c = matrix_result;

    for(int i = 0; i < num_row; ++i)
    {
       for(int j = 0; j < num_col; ++j)
        {
            __m128 *m3 = (__m128*)a[i];     // The error is here.
            __m128 *m4 = (__m128*)b[j];
            float* res;
            *(c + (j * num_col + i)) = 0;
          for(int k = 0; k < num_col; k += 4)
            {
                __m128 m5 = _mm_mul_ps(*m3,*m4);
                res = (float*)&m5;
                *(c + (j * num_col + i)) += res[0]+res[1]+res[2]+res[3];
                m3++;
                m4++;
            }
        }
    }
}

Answer 1

我假设它是__m128 *m3 = (__m128*)a[i]; 和之后的行会产生错误。您正在尝试将 float 类型转换为指向 __m128 的 指针，编译器可以抱怨该类型。

我不知道预期算法的细节。假设意图是访问四个浮点数 a[i]..a[i+3] 作为单个 __m128，你需要这样的东西：

__m128 *m3 = (__m128*)&a[i];
__m128 *m4 = (__m128*)&b[j];

或同等学历：

__m128 *m3 = (__m128*)(a + i);
__m128 *m4 = (__m128*)(b + j);