SoloCodigo

Hola,

Estoy trabajando con arreglos de 4 dimensiones y estoy definiendo los arreglos con la funcion malloc. Sin embargo, cuando incremento el tamano de alguno de los elementos, obtengo un segmentation fault pero no logro conseguir el por que. A continuacion pondre como defini mi funcion a ver si alguien puede ayudarme:

De esta manera defino lasa matrices:
    double ****Fd, ****Fdnew;
    double ****gxFd, ****gxFdnew;
    double ****gyFd, ****gyFdnew;
    double ****gvxFd, ****gvxFdnew;
    double ****gvyFd, ****gvyFdnew;

    double ****FIn;
    double ****gxFIn;
    double ****gyFIn;
    double ****gvxFIn;
    double ****gvyFIn;

    double ****Fo;
    double ****Fa, ****gxFa, ****gvxFa;
    double ****Fo1;  
    double ff, ff2, gg, gg2, vv, vv2;
   
            Fd = (double****)malloc(sizeof(double***)*ymax+2);
            Fdnew = (double****)malloc(sizeof(double***)*ymax+2);
            FIn = (double****)malloc(sizeof(double***)*ymax+2);
            Fa = (double****)malloc(sizeof(double***)*ymax+2);
            Fo = (double****)malloc(sizeof(double***)*ymax+2);
            Fo1 = (double****)malloc(sizeof(double***)*ymax+2);
           
            gxFd = (double****)malloc(sizeof(double***)*ymax+2);
            gxFIn = (double****)malloc(sizeof(double***)*ymax+2);
            gxFa = (double****)malloc(sizeof(double***)*ymax+2);
            gxFdnew = (double****)malloc(sizeof(double***)*ymax+2);
           
            gyFd = (double****)malloc(sizeof(double***)*ymax+2);
            gyFIn = (double****)malloc(sizeof(double***)*ymax+2);
            gyFdnew = (double****)malloc(sizeof(double***)*ymax+2);
           
            gvxFd = (double****)malloc(sizeof(double***)*ymax+2);
            gvxFIn = (double****)malloc(sizeof(double***)*ymax+2);
            gvxFdnew = (double****)malloc(sizeof(double***)*ymax+2);
           
            gvyFd = (double****)malloc(sizeof(double***)*ymax+2);
            gvyFIn = (double****)malloc(sizeof(double***)*ymax+2);
            gvyFdnew = (double****)malloc(sizeof(double***)*ymax+2);
           
            for(i=0;i<ymax+2;i++){
            Fd = (double***)malloc(sizeof(double**)*xmax+2);
            Fdnew = (double***)malloc(sizeof(double**)*xmax+2);
            Fo = (double***)malloc(sizeof(double**)*xmax+2);
            Fo1 = (double***)malloc(sizeof(double**)*xmax+2);
            FIn = (double***)malloc(sizeof(double**)*xmax+2);
            Fa = (double***)malloc(sizeof(double**)*xmax+2);
           
            gxFd = (double***)malloc(sizeof(double**)*xmax+2);
            gxFIn = (double***)malloc(sizeof(double**)*xmax+2);
            gxFa = (double***)malloc(sizeof(double**)*xmax+2);
            gxFdnew = (double***)malloc(sizeof(double**)*xmax+2);
           
            gyFd = (double***)malloc(sizeof(double**)*xmax+2);
            gyFIn = (double***)malloc(sizeof(double**)*xmax+2);
            gyFdnew = (double***)malloc(sizeof(double**)*xmax+2);
           
            gvxFd = (double***)malloc(sizeof(double**)*xmax+2);
            gvxFIn = (double***)malloc(sizeof(double**)*xmax+2);
            gvxFdnew = (double***)malloc(sizeof(double**)*xmax+2);
           
            gvyFd = (double***)malloc(sizeof(double**)*xmax+2);
            gvyFIn = (double***)malloc(sizeof(double**)*xmax+2);
            gvyFdnew = (double***)malloc(sizeof(double**)*xmax+2);
             
             
            for(j=0;j<xmax+2;j++){
            Fd[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            Fdnew[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            Fo[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            Fo1[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            FIn[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            Fa[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
           
            gxFd[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            gxFdnew[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            gxFIn[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            gxFa[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
           
            gyFd[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            gyFdnew[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            gyFIn[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
           
            gvxFd[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            gvxFdnew[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            gvxFIn[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
           
            gvyFd[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            gvyFdnew[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
            gvyFIn[j] = (double**)malloc(sizeof(double*)*(2*Nmax+2));
           
            for(p=0;p<2*Nmax+2;p++){
            Fd[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            Fdnew[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            Fo[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            Fo1[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            FIn[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            Fa[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
           
            gxFd[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            gxFdnew[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            gxFIn[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            gxFa[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
           
            gyFd[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            gyFdnew[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            gyFIn[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
           
            gvxFd[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            gvxFdnew[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            gvxFIn[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
           
            gvyFd[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            gvyFdnew[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
            gvyFIn[j][p] = (double*)malloc(sizeof(double)*(2*Mmax+2));
           
           
            for(q=0;q<2*Mmax+2;q++){
            Fd[j][p][q] = 0.0;
            Fdnew[j][p][q] = 0.0;
            Fo[j][p][q] = 0.0;
            Fo1[j][p][q] = 0.0;
            FIn[j][p][q] = 0.0;
            Fa[j][p][q] = 0.0;
           
            gxFd[j][p][q] = 0.0;
            gxFdnew[j][p][q] = 0.0;
            gxFIn[j][p][q] = 0.0;
            gxFa[j][p][q] = 0.0;
           
            gyFd[j][p][q] = 0.0;
            gyFdnew[j][p][q] = 0.0;
            gyFIn[j][p][q] = 0.0;
           
            gvxFd[j][p][q] = 0.0;
            gvxFdnew[j][p][q] = 0.0;
            gvxFIn[j][p][q] = 0.0;
           
            gvyFd[j][p][q] = 0.0;
            gvyFdnew[j][p][q] = 0.0;
            gvyFIn[j][p][q] = 0.0;
           
             }                    
             }
             }
             }

Y esto las uso en una funcion como la siguiente:
void Init(double**** Fd, double**** gxFd, double****gyFd, double**** gvxFd, double**** gvyFd, double**** FIn,
double**** gxFIn, double**** gyFIn, double**** gvxFIn, double**** gvyFIn, double**** Fo, double****Fo1,
double* Vx, double* Vy, double nx)
{
   int i, j, p, q;
   
   for(i=0; i<=xmax+1; i++){ //Inicialmente lleno TODO x
        for(j=0; j<=ymax+1; j++){
             for(p=0; p<=2*Nmax+1; p++){
                  for(q=0; q<=2*Mmax+1; q++){
     
      Fo[j][p][q]= nx*(me/(2.0*Pi*k*Tg))*exp(-(me*((Vx[p]*Vx[p])+(Vy[q]*Vy[q]))/(2.0*k*Tg)));
      Fo1[j][p][q]= nx*(me/(2.0*Pi*k*Tg))*exp(-(me*(((Vx[p])*(Vx[p]))+(Vy[q]*Vy[q]))/(2.0*k*Tg))); //Fo desplazado pero tg
     
      Fd[j][p][q]= nx*(me/(2.0*Pi*k*Tei))*exp(-(me*(((Vx[p]-Vx[33])*(Vx[p]-Vx[33]))+(Vy[q]*Vy[q]))/(2.0*k*Tei)));
      gxFd[j][p][q]= 0.0;
      gyFd[j][p][q]= 0.0;
      gvxFd[j][p][q]= -((2.0*me*((Vx[p]-Vx[33])))/(k*Tei))*nx*sqrt(me/(2.0*Pi*k*Tei))*exp(-(me*(((Vx[p]-Vx[33])*(Vx[p]-Vx[33]))+(Vy[q]*Vy[q]))/(2.0*k*Tei)));
      gvyFd[j][p][q]= -((2.0*me*(Vy[q]))/(k*Tei))*nx*sqrt(me/(2.0*Pi*k*Tei))*exp(-(me*(((Vx[p]-Vx[33])*(Vx[p]-Vx[33]))+(Vy[q]*Vy[q]))/(2.0*k*Tei)));
     
      //Injection Function at the Inlet of the channel
      FIn[j][p][q]= nx*(me/(2.0*Pi*k*Tei))*exp(-(me*(((Vx[p]-Vx[33])*(Vx[p]-Vx[33]))+(Vy[q]*Vy[q]))/(2.0*k*Tei)));
      gxFIn[j][p][q]= 0.0;
      gyFIn[j][p][q]= 0.0;
      gvxFIn[j][p][q]= -((2.0*me*((Vx[p]-Vx[33])))/(k*Tei))*nx*sqrt(me/(2.0*Pi*k*Tei))*exp(-(me*(((Vx[p]-Vx[33])*(Vx[p]-Vx[33]))+(Vy[q]*Vy[q]))/(2.0*k*Tei)));
      gvyFIn[j][p][q]= -((2.0*me*(Vy[q]))/(k*Tei))*nx*sqrt(me/(2.0*Pi*k*Tei))*exp(-(me*(((Vx[p]-Vx[33])*(Vx[p]-Vx[33]))+(Vy[q]*Vy[q]))/(2.0*k*Tei)));//*/
     
                    }
                  }
                }      
              }
     }

Donde:
#define xmax 70
#define ymax 50
#define Nmax 30
#define Mmax 30

Se que es un poco tedioso. El programa funciona bien si xmax vale 50 o incluso con 60 funciona. El error se mantiene aun usando una computadora de mayor capacidad asi que no creo que ese sea el problema.