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);

使用SSE错误__m128到浮点转换的矩阵乘法？

1 个答案: