编辑我删除了问题的一部分，因为它可能会令人困惑，感谢答案，我应该能够自己发现。

我试图写一个程序，使用频谱函数对信号进行重采样。我用C写代码，但我发现用C++编码的问题。
下面是我使用的部分代码。我有一个函数test()，它调用另一个函数test3(tmp)，传递一个数组作为指针。只要我合并组合+=运算符，并试图填充第一个函数的数组，我的代码就会变慢。
我初始化和计算使用零，但所有数组必须有双精度值。

#include <stdio.h>
#include <stdlib.h>

void test3(double *tmp);

void test() {
    double *tmp;
    int k;
    tmp = (double *) calloc(250 * 6, sizeof(double));
    for (k = 0; k < 1000; ++k)
        test3(tmp);
    free(tmp);
}

void test3(double *tmp) {
    int k;
    double *fx;
    double *srf;
    double *ptr;
    int p;
    int l;
    int c;
    double a;
    fx = (double *) calloc(2100 * 6, sizeof(double));
    srf = (double *) calloc(2100 * 250, sizeof(double));
    
    for (int k = 0; k < 2100 * 6; ++k)   fx[k] = 0;
    for (int k = 0; k < 2100 * 250; ++k) srf[k] = 0;
    
    ptr = (double *) calloc(250 * 6, sizeof(double));
    
    for (p = 0; p < 6; ++p) {
        for (l = 0; l < 250; ++l) {
            ptr[l + p*250] = 0.0;
            for (c = 0; c < 2100; ++c)
                ptr[l+p*250] += fx[c+p*2100] * srf[c+l*2100];
        }
    }
    for (k = 0; k < 250 * 6; ++k) tmp[k] = ptr[k];
    
    free(ptr);
    free(fx);
    free(srf);
}

int main(int argc, char *argv[]) {
    test();
    return 0;
}

字符串
我编译使用：

gcc -O2 -o test test.c -lm

型
当我运行代码time ./test时，

./test_lut  3.84s user 0.11s system 98% cpu 4.008 total

型