尝试写入矩阵后的 C++ 分段错误

Question

我将这个 3D 矩阵分配为一块内存，但是当我尝试写入该死的东西时，它给了我一个分段错误。这个东西在两个维度上都很好，但由于某种原因，我在第三个维度上遇到了麻烦……我不知道分配中的错误在哪里。对我来说它看起来很完美。

代码如下：

phi = new double**[xlength];
phi[0] = new double*[xlength*ylength];
phi[0][0] = new double[xlength*ylength*tlength];
for (int i=0;i<xlength;i++)
{
    phi[i] = phi[0] + ylength*i;
    for (int j=0;j<ylength;j++)
    {
        phi[i][j] = phi[i][0] + tlength*j;
    }
}

任何帮助将不胜感激。 （是的，我想要一个 3D 矩阵）

另外，如果重要的话，这就是我得到分段错误的地方：

for (int i = 0; i < xlength; i++)
    {
        for (int j = 0; j < ylength; j++)
        {
            phi[i][j][1] = 0.1*(4.0*i*h-i*i*h*h)
            *(2.0*j*h-j*j*h*h);
        }
    }

这确实适用于二维！

phi = new double*[xlength];
phi[0] = new double[xlength*ylength];
for (int i=0;i<xlength;i++)
{
    phi[i] = phi[0] + ylength*i;
}

Answer 1

您没有分配其他子矩阵，例如phi[1] 或 phi[0][1]

你至少需要

phi = new double**[xlength];
for (int i=0; i<xlength; i++) { 
    phi[i] = new double* [ylength];
    for (int j=0; j<ylength; j++) { 
       phi[i][j] = new double [zlength];
       for (k=0; k<zlength; k++) phi[i][j][k] = 0.0;
    }
}

你应该考虑使用std::vector（或者甚至，如果在C++ 2011中，std::array），即

std::vector<std::vector<double> > phi;

然后使用std::vector，您需要phi.resize(xlength) 和一个循环来调整每个子元素phi[i].resize(ylength) 等的大小。

---

如果你想一次分配所有内存，你可以有

double* phi = new double[xlength*ylength*zlength]

但是你不能使用phi[i][j][k] 符号，所以你应该

#define inphi(I,J,K) phi[(I)*xlength*ylength+(J)*xlength+(K)]

写inphi(i,j,k)而不是phi[i][j][k]

---

你的第二个代码不起作用：它是undefined behavior（它不会崩溃，因为你很幸运，它可能会在其他系统上崩溃......），只是一些内存泄漏还没有崩溃（但是以后可能会崩溃，甚至可能再次重新运行程序）。使用像valgrind这样的内存泄漏检测器