#include <stdio.h> #include <time.h> #include <stdlib.h> #include <math.h> #include <smmintrin.h> #include <emmintrin.h> #include <immintrin.h> #define MAXN 100000000 float a[MAXN]; float b[MAXN]; float c[MAXN]; float d[MAXN]; int main() { for (int i = 0; i < MAXN; ++i) { a[i] = 1.0 / (rand() + 1); b[i] = 1.0 / (rand() + 1); } for (int n = 0; n < 20; ++n) { for (int i = 0; i < MAXN; ++i) { d[i] += a[i] * b[i]; } } clock_t start, end; start = clock(); for (int i = 0; i < MAXN; i += 8) // 假设MAXN是8的倍数 { __m256d va = _mm256_load_pd(&a[i]); // 加载8个双精度浮点数到寄存器 __m256d vb = _mm256_load_pd(&b[i]); __m256d vc = _mm256_load_pd(&c[i]); __m256d result = _mm256_mul_pd(va, vb); // 向量乘法运算 result = _mm256_add_pd(result, vc); // 向量加法运算 _mm256_store_pd(&c[i], result); // 存储结果回内存 } end = clock(); printf("time=%f\n", (double)(end - start) / CLOCKS_PER_SEC); for (int i = 0; i < MAXN; ++i) { if (fabs(c[i] - d[i]) / d[i] > 0.0001) { printf("Check Failed at %d\n", i); return 0; } } printf("Check Passed"); } 34 38 C:\Users\Administrator\Desktop\add.cpp [Error] cannot convert 'float*' to 'const double*' for argument '1' to '__m256d _mm256_load_pd(const double*)'

这段代码出现了一个错误。错误的原因是在使用SSE指令集进行向量化计算时，加载数据时使用了错误的数据类型。在这段代码中，数组a、b、c、d的数据类型是float，但是在使用`_mm256_load_pd`函数加载数据时，需要传入`const double*`类型的参数。因此，需要将加载数据的部分修改为使用`_mm256_load_ps`函数来加载单精度浮点数。同时，由于AVX指令集对齐要求更高，需要确保数组a、b、c、d的地址是16字节对齐的。

以下是修改后的代码：

#include <stdio.h>
#include <time.h>
#include <stdlib.h>
#include <math.h>
#include <smmintrin.h>
#include <emmintrin.h>
#include <immintrin.h>

#define MAXN 100000000

float a[MAXN];
float b[MAXN];
float c[MAXN];
float d[MAXN];

int main()
{
    for (int i = 0; i < MAXN; ++i)
    {
        a[i] = 1.0 / (rand() + 1);
        b[i] = 1.0 / (rand() + 1);
    }
    for (int n = 0; n < 20; ++n)
    {
        for (int i = 0; i < MAXN; ++i)
        {
            d[i] += a[i] * b[i];
        }
    }
    clock_t start, end;
    start = clock();
    for (int i = 0; i < MAXN; i += 8) // 假设MAXN是8的倍数
    {
        __m256 va = _mm256_load_ps(&a[i]); // 加载8个单精度浮点数到寄存器
        __m256 vb = _mm256_load_ps(&b[i]);
        __m256 vc = _mm256_load_ps(&c[i]);

        __m256 result = _mm256_mul_ps(va, vb); // 向量乘法运算
        result = _mm256_add_ps(result, vc);     // 向量加法运算

        _mm256_store_ps(&c[i], result); // 存储结果回内存
    }
    end = clock();
    printf("time=%f\n", (double)(end - start) / CLOCKS_PER_SEC);
    for (int i = 0; i < MAXN; ++i)
    {
        if (fabs(c[i] - d[i]) / d[i] > 0.0001)
        {
            printf("Check Failed at %d\n", i);
            return 0;
        }
    }
    printf("Check Passed");
    return 0;
}

请注意，我对加载数据和存储结果的部分进行了修改，并将数据类型从`__m256d`改为了`__m256`，以适应单精度浮点数。另外，代码中的对齐问题需要注意，可以使用`posix_memalign`或`_aligned_malloc`等方法来确保数组的对齐。