DAXPY优化代码触发Segmentation Fault，无法定位38-39行错误

错误根源：栈溢出 #

你遇到的Segmentation fault本质是栈空间不足导致的溢出。在main函数中，你声明了两个变长数组（VLA）：

double x[n];
double y[n];

当n=10000000时，每个double占8字节，两个数组总共需要约160MB内存。而操作系统给进程分配的栈空间通常很小（比如macOS默认栈大小仅8MB左右），这么大的数组直接分配在栈上会超出栈的容量，写入数组时就会触发EXC_BAD_ACCESS错误。

解决方案 #

方案1：改用堆内存分配（推荐，可移植性强） #

使用malloc在堆上分配内存，堆的可用空间远大于栈，适合存储大数据集。修改后的main函数代码如下：

int main() {
    int n = 10000000;
    double a = 2.0;
    // 在堆上分配内存
    double *x = malloc(n * sizeof(double));
    double *y = malloc(n * sizeof(double));
    
    // 检查内存分配是否成功
    if (x == NULL || y == NULL) {
        fprintf(stderr, "内存分配失败\n");
        return 1;
    }
    
    int i;
    for (i = 0; i < n; ++i) {
        x[i] = (double)i;
        y[i] = (double)i;
    }
    clock_t start = clock();
    daxpy(n, a, x, y);
    clock_t end = clock();
    double time_elapsed = (double)(end - start) / CLOCKS_PER_SEC;
    printf("Time elapsed: %f seconds\n", time_elapsed);
    
    // 释放堆内存
    free(x);
    free(y);
    return 0;
}

方案2：临时增大栈大小（不推荐，依赖平台） #

如果一定要使用栈上数组，可以通过编译器参数临时调整栈大小。比如在macOS下编译时：

gcc daxpy.c -o daxpy.x -Wl,-stack_size,0x10000000 # 设置栈大小为256MB

这种方法仅适用于特定平台，代码可移植性会大幅降低，不建议在生产环境使用。

额外优化补充 #

你的DAXPY循环展开逻辑存在小缺陷：当n不是4的整数倍时，最后剩余的1-3个元素不会被处理。可以在循环后添加一个小循环处理剩余元素：

void daxpy(int n, double a, double *x, double *y) {
    int i;
    double y0, y1, y2, y3, x0, x1, x2, x3;
    // 循环展开处理4个元素的批次
    for (i = 0; i < n - 3; i += 4) {
        x0 = x[i];
        x1 = x[i + 1];
        x2 = x[i + 2];
        x3 = x[i + 3];

        y0 = y[i];
        y1 = y[i + 1];
        y2 = y[i + 2];
        y3 = y[i + 3];

        y0 += a * x0;
        y1 += a * x1;
        y2 += a * x2;
        y3 += a * x3;

        y[i] = y0;
        y[i + 1] = y1;
        y[i + 2] = y2;
        y[i + 3] = y3;
    }
    // 处理剩余不足4个的元素
    for (; i < n; ++i) {
        y[i] += a * x[i];
    }
}

内容的提问来源于stack exchange，提问作者user3116936