在使用矩阵进行计算时擦除 C++ 2d 向量行

Question

我正在尝试模拟一组神经元的平均行为。这意味着我需要对包含数十亿个元素的矩阵进行计算（steps~10^₆、neurons~10 ^₄）。

为了避免吃掉我的整个 RAM（并死去尝试），我决定在完成计算后立即从矩阵中删除行。我对 C++ 没有太多经验，但我的理解是 v.erase( v.begin()-i+1); 应该允许我这样做。

// Membrane potential matrix using st::vector
vector<vector<double>> v;
v.resize(steps + 1, vector<double>(neurons));
// Initialise v
for (size_t n = 0; n < neurons; n++) {
    v[0][n] = v0;
}

double v_avg[steps + 1] = {v0};

// Loop
for (size_t i = 1; i < steps + 1; i++) {
    for (size_t n = 0; n < neurons; n++) {
        if(v[i-1][n] >= vp) {
            v[i][n] = -vp;
        }
        else {
            v[i][n] = v[i-1][n] + h * ( pow(v[i-1][n], 2) + I[i] + eta[n] );
        }
        v_avg[i] += v[i][n]; // Sum of membrane potentials
    }
    cout << "step " << i << "/" << steps << " done\n";
    v.erase( v.begin()-i+1); // Erase row v[i-1]
    v_avg[i] = v_avg[i]/neurons; // Mean membrane potential
}
v.erase( v.begin()+steps+1 ); // Erase last row

我不确定为什么在 steps/2 步骤之后会出现分段错误（我正在使用较小的 steps 值进行测试）：

...    
step 10/20 done
[1]    1791 segmentation fault (core dumped)  ./qif_solve_vect

---

更新：

感谢@1201ProgramAlarm 我明白了我的问题所在。我的问题是：

1. 如何以从一开始就没有分配的方式使用矩阵。
2. 如何在保留索引的同时释放/释放行（与 v.erase( v.begin()) 不同）。这是必不可少的，因为稍后我将在每个神经元产生尖峰时为它们实现不同的不应时间 (v[i][n] = -vp;)。

Answer 1

在您的erase 语句中，您从v.begin() 中减去，这将导致无效的迭代器，因为它将指向向量开始之前。你的意思可能是v.erase( v.begin() + i - 1);。

但是，由于您已经分配了完整的矩阵，因此像这样擦除并不能节省任何空间。擦除会将所有剩余元素向下移动一个元素，并且您对下一个循环的索引将是错误的（因为您想一直使用v[0]）。