0

我将这个 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;
}
4

1 回答 1

0

您没有分配其他子矩阵，例如 egphi[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]

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

于 2013-05-15T04:48:32.870 回答