Question

我用C ++编写一个简单的代码，然后在linux ubuntu 11.04上用g ++编译它，我没有得到任何错误但是当我运行可执行文件时，我得到了这个错误“seg fault”。

我知道我的代码没有问题，而且这个错误与编译器有关。

有人能帮助我吗？

我的代码是：

#include <math.h>
int main()
{
    double a[200][200][200],b[200][200][200],c[200][200][200];
    int i,j,k;

    double const pi=3.14;

    for(k=0;k<200;k++)
    {
        for(j=0;j<200;j++)
        {
            for(i=0;i<200;i++)
            {
                a[i][j][k]=sin(1.5*pi*i)*cos(3.5*pi*j)*k;
                b[i][j][k]=cos(1.5*pi*i)*cos(2.5*pi*k)*j;
                c[i][j][k]=a[i][j][k]-b[i][j][k];
            }
        }
    }
}

Answer 1

这三个阵列需要大约190MB的空间，这几乎肯定超过了操作系统强加的堆栈大小限制。

尝试在堆上分配它们（使用new），而不是将它们放在堆栈上。

Answer 2

你将大量的double数组放入堆栈（假设你的架构是如何处理局部变量的）。几乎可以肯定，你的系统堆栈不能容纳那么多空间。

相反，使用向量代替在堆上分配：

std::vector<std::vector<std::vector<double> > > a(200, std::vector<std::vector<double> >(200, std::vector<double>(200)));

Answer 3

堆栈溢出 - ＆gt;分段错误

Answer 4

对于非对称的情况，我们必须使用它：

    double*** a=new double**[IE];
    for(int i=0;i<IE;i++){
        a[i]=new double *[JE];
        for(int j=0;j<JE;j++){
            a[i][j]=new double [KE];
        }
    }

使用上面的代码，我们可以构建巨大的矩阵，因为我们的计算机的ram允许。

Answer 5

此功能可以帮助您：

double ***alloc3d(int l, int m, int n) {
    double *data = new double [l*m*n];
    double ***array = new double **[l];
    for (int i=0; i<l; i++) {
        array[i] = new double *[m];
        for (int j=0; j<m; j++) {
            array[i][j] = &(data[(i*m+j)*n]);
        }
    }

    return array;
}

大矩阵计算

5 个答案: