C++数据分析项目处理粒子帧数据时trimerCands内存访问错误求助

核心错误：数组索引越界覆盖栈内存 #

你遇到的内存错误是典型的栈缓冲区溢出，直接覆盖了栈上存储trimerCands对象内部指针的内存区域，导致后续访问trimerCands时触发非法内存访问。

具体问题点： #

• 你声明的trimerIDs数组第一个维度长度是当前帧的三聚物数量nTrimers，合法下标范围是0 ~ nTrimers -1，但你直接用全局索引i（取值范围是firstTrim ~ lastTrim，经过17帧迭代后i的数值已经远大于nTrimers）作为第一个维度的下标，直接越界写入栈内存。
• int trimerIDs[nTrimers][3] = {0}; 是C标准的变长数组（VLA），不属于C++标准规范，属于编译器扩展，且分配在栈上，当nTrimers较大时本身就容易触发栈溢出。

其他需要修复的问题： #

• 函数参数std::vector<trimerCoords> trimerCands是值传递，每次调用函数都会完整拷贝整个大体积的vector，不仅性能极差，如果你后续需要修改原vector的话也不会生效，建议改为const std::vector<trimerCoords>& trimerCands传递常量引用。
• 代码中用到的framesTBA看起来是全局变量，函数内没有定义也没有作为参数传入，逻辑上不安全，建议改为函数参数显式传入。
• 坐标匹配的双层循环时间复杂度是O(n²)，当单帧粒子数量较多时性能会很差，后续可以考虑用空间分区或者哈希表优化匹配逻辑。

最小修复代码示例： #

// 替换掉变长数组声明，用符合C++标准的vector+array
std::vector<std::array<int, 3>> trimerIDs(nTrimers, {0});

// 循环内索引计算修改
for (int i = firstTrim; i <= lastTrim; i++){
    // 计算当前三聚物在当前帧内的相对下标
    int local_idx = i - firstTrim;
    tempEntry.x = trimerCands[i].x1;
    tempEntry.y = trimerCands[i].y1;
    found = false;
    for (long unsigned int j = 0; j < currentFound.size(); j++){
        if (fabs(tempEntry.x - currentFound[j].x) + fabs(tempEntry.y - currentFound[j].y) < 0.001){
            trimerIDs[local_idx][0] = j; found = true; break;
        }
    }
    if (found == false) {
        currentFound.push_back(tempEntry);
        traceLengths.push_back(trimerCands[i].tLength1);
        trimerIDs[local_idx][0] = particleCounter;
        particleCounter++;
    }
}

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