Question

我正在编写一个代码，目的是操纵来自3D数组的一维数组。显示代码将解释所有内容：

首先，我正在使用结构数组：

typedef struct range_in_memory {
      double E, R;
} RANGE;

我在简单的调试实现中调用了这个函数，就像魅力一样：

    RANGE *dq_Eloss_load_range_file (double Aion, double Zion, double Atar, double Ztar, int *n){

      char *filename;
      char errormsg[80], dummy;
      int N = 2;
      int i = 0;
      FILE *fp;
      RANGE *memrange;

      filename = (char*)malloc (50*sizeof (char));
      dq_Eloss_set_filename (filename, Aion, Zion, Atar, Ztar);  
      fp = dq_myfopen (filename, "r", errormsg);

      while ( !feof (fp) ){
        fscanf (fp, "%c", &dummy);

        if (dummy == '\n') {
          N++;

        }
      }

      rewind (fp);

      memrange = (RANGE*) calloc (N, sizeof (RANGE));

      while ( !feof (fp) ){
        fscanf (fp, "%lf\t%lf\n", &memrange[i].E, &memrange[i].R);
        i++;
      }
      *n = N;

    //   for (i=0; i<N; i++){
    //         printf ("\n%lf  %lf", memrange[i].E, memrange[i].R);
    //   }

      fclose (fp);

      return (memrange);
}

评论周期至关重要，因此请记住。现在，遇到我的问题，我需要使用3D数组：

RANGE ***memrangeTAR;

函数*dq_Eloss_load_range_file在RANGE类型的数组中加载文件的内容。我想对几个文件执行此操作，这些文件根据两个参数Aion和Zion命名。在我看来，这应该是这样的：

(double) memrangeTAR [Aion][Zion][i].E
(double) memrangeTAR [Aion][Zion][i].R

这将是与Aion和Zion相关的文件中第i行的值。

这是代码：

void dq_load_range_files (RANGE ***memrangeTAR){
  int NTAR;
  double *ZMAX;
  int Aion, Zion, AMAX = 250.;
  int i;



  printf ("\n   Allocating Memory for range data..."); fflush (stdout);
  ZMAX = (double*) malloc ((AMAX+1)*sizeof (double));
  memrangeTAR   = (RANGE***) malloc ((AMAX+1)*sizeof (RANGE**));


  for (Aion=1; Aion<=AMAX; Aion++){
    ZMAX[Aion] = ceil (dq_range_table_get_max_Zion (Aion)); 
    memrangeTAR  [Aion] = (RANGE**) malloc ((ZMAX[Aion]+1)*sizeof(RANGE*)); 
  }

  printf ("   ...ALLOCATED!                                  \n\n");


  printf ("\n   Loading range data...\n"); fflush (stdout);

  for (Aion=1; Aion <=AMAX; Aion++){

        for (Zion=1; Zion<=ZMAX[Aion]; Zion++){

      memrangeTAR  [Aion][Zion] = dq_Eloss_load_range_file ((double)Aion, (double)Zion, cor_sosau16.mass.Mtar,  cor_sosau16.charge.Ztar,  &NTAR); 

//   for (i=0; i<N; i++){
//         printf ("\n%lf  %lf", memrangeTAR[Aion][Zion][i].E, memrangeTAR[Aion][Zion][i].R);
//   }   

    }
  }

}

它编译（使用gnu99标准在unix机器上的gcc）。除非我尝试访问3D阵列，否则它会运行。如果我在第一个函数中取消注释打印周期，我可以得到我想要的东西：它在终端上打印当前文件的内容。如果我在第二个函数中取消注释for循环，它会给我分段错误。

我做错了什么？

Answer 1

我设法解决但不明白为什么它是一个解决方案。通过定义

RANGE *tmp;

变量，将通常的函数调用为

tmp = dq_Eloss_load_range_file ((double)Aion, (double)Zion, cor_sosau16.mass.Mtar, cor_sosau16.charge.Ztar, &NTAR);

然后分配想要的指针：

memrangeTAR[Aion][Zion] = tmp;

我现在可以以任何方式使用变量的内容，例如：

   for (i=0; i<N; i++){
        printf ("\n%lf  %lf", memrangeTAR[Aion][Zion][i].E, 
                              memrangeTAR[Aion][Zion][i].R);
   }

通过这种方式没有遇到分段错误错误，让我觉得问题是我将指针传递给函数的方式，但仍然不明白......

使用3D动态阵列的功能中的分段故障错误

1 个答案: