OpenMP Task实现矩阵加法代码出现段错误，请求协助排查

以下代码执行时出现段错误：

void x(vector<vector<int>>&A,vector<vector<int>>&B,vector<vector<int>>&C,int n)
{
int i,j;
#pragma omp parallel
    #pragma omp single
    {
       for (i=0;i<n;++i)
         for(j=0;j<n;++j)
          #pragma omp task
           C[i][j] = A[i][j]+B[i][j];
    }
}

提问者认为parallel指令会创建多线程，single区域会被所有线程访问，为每个(i,j)元素创建n*n个独立任务，因此不存在竞态条件，但代码仍出现段错误，寻求解决。

• 纠正对single指令的误解：#pragma omp single并非让所有线程进入，而是仅单个线程执行该区域代码，其余线程会在区域末尾等待（除非添加nowait子句）。不过这不是段错误的直接诱因。
• 段错误的核心原因：循环变量作用域问题
  代码中i和j是函数级别的变量，当#pragma omp task创建任务时，任务不会立即执行，而是进入队列等待调度。而循环会持续修改i和j的值，当任务最终执行时，它访问的i、j早已不是创建任务时的迭代值——甚至可能超出n的范围，导致访问A、B、C时出现数组越界，触发段错误。
• 修复方案：缩小循环变量作用域
  将i和j的声明移到循环内部，让每个任务捕获当前迭代的变量副本：

void x(vector<vector<int>>&A,vector<vector<int>>&B,vector<vector<int>>&C,int n)
{
#pragma omp parallel
    #pragma omp single
    {
       for (int i=0;i<n;++i)
         for(int j=0;j<n;++j)
          #pragma omp task
           C[i][j] = A[i][j]+B[i][j];
    }
}

这样每个任务创建时都会绑定当前迭代的i、j值，避免了后续循环修改变量导致的访问越界。

• 额外检查：容器初始化
  需确保传入的A、B、C已正确初始化为n×n的二维vector。如果C的行/列数不足n，即使修复了变量作用域，仍会因越界访问触发段错误。示例初始化代码：

int n = 100;
vector<vector<int>> A(n, vector<int>(n));
vector<vector<int>> B(n, vector<int>(n));
vector<vector<int>> C(n, vector<int>(n));
x(A, B, C, n);

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