コアダンプをなくしたい

Question

発生している問題・エラーメッセージ

Amideを用いて3次元画像を2次元画像に投影するプログラムを作成しています。コンパイルはできるのですが実行するとコアダンプが発生してしまいます。
初心者なのでなるべくかみ砕いて説明していただけるとありがたいです。

該当のソースコード

#include<stdio.h>
#include<stdlib.h>
#define nx 128
#define ny 128
#define nz 128

#define MASS_ATT_TISSUE 0.2264
#define DENSITY_TISSUE 1.06

#define BMASS_ATT_TISSUE 0.4242
#define BDENSITY_TISSUE 1.92
#define PIXEL 0.25

int main(void){
  
  int i,j,k;
  FILE *fp;

  float att_tis;
  float batt_tis;
  float img_att[nz][ny][nx];
  float proj_att[ny][nx];
  float att_tis_pix;
  float batt_tis_pix;
  
  att_tis=MASS_ATT_TISSUE*DENSITY_TISSUE;
  att_tis_pix=att_tis*PIXEL;

  batt_tis=BMASS_ATT_TISSUE*BDENSITY_TISSUE;
  batt_tis_pix=batt_tis*PIXEL;
  
  for(i=0;i<nx;i++){
    for(j=0;j<ny;j++){
      for(k=0;k<nz;k++){
	img_att[k][j][i]=0;
      }
    }
  }
  
  float x0,y0,a,b;
  float X,Y;
  
  //3D作成

  for(k=0;k<nz;k++){

    x0=64;y0=64;
    a=56;b=40;
    for(i=0;i<nx;i++){
      for(j=0;j<ny;j++){
	X=((float)i-x0)*((float)i-x0)/(a*a);
	Y=((float)j-y0)*((float)j-y0)/(b*b);
	if(X+Y <= 1){
	  img_att[k][j][i]=att_tis_pix;
	}
      }
    }

  }
 
  for(k=4;k<nz;k++){

    x0=40;y0=64;
    a=20;b=24;
    for(i=0;i<nx;i++){
      for(j=0;j<ny;j++){
	X=((float)i-x0)*((float)i-x0)/(a*a);
	Y=((float)j-y0)*((float)j-y0)/(b*b);
	if(X+Y <= 1){
	  img_att[k][j][i]=0;
	}
      }
    }
    x0=88;y0=64;
    a=20;b=24;
    for(i=0;i<nx;i++){
      for(j=0;j<ny;j++){
	X=((float)i-x0)*((float)i-x0)/(a*a);
	Y=((float)j-y0)*((float)j-y0)/(b*b);
	if(X+Y <= 1){
	  img_att[k][j][i]=0;
	}
      }
    }
  }
  


  //投影
  for(k=0;k<nz;k++){
    for(j=0;j<ny;j++){
      for(i=0;i<nx;i++){
	proj_att[k][i]+=img_att[k][j][i];
      }
    }
  }
  
  fp=fopen("proj_3DCT.raw","wb");
  if(fp == NULL)
    {
      printf("File open error\n");
      exit(1);
    }
    for(j=0;j<ny;j++)
      {
	fwrite(proj_att[j],sizeof(float),nx,fp);
      }
  fclose(fp);
  
  return 0;
}

Comments

[hoshi-takanori]

コアダンプが発生するソースコードの位置は分かりますか？ または、proj_3DCT.raw というファイルが公開されているなら、こちらで試してみますけど…。

Answer 1

ローカル変数に大きな容量の変数を定義してます。
ローカル変数はスタックに取られますんで、ここらへんの変数をグローバル変数としてみてはどうでしょう

Answer 2

たいがいのコンピュータ環境において、「スタック領域」はあまり大きくありません。通常のPC上のOSであれば2MB以下であることが多いです。
その範囲を超えるような操作をすると、プログラムが破綻して何が起こるかわかりませんが、C/C++ではそういう操作も「プログラマの責任」において禁止されず出来てしまいます。

で、ローカル変数はそのスタック領域上に確保されます。
float img_att[nz][ny][nx];
nz,ny,nxが各々128であれば、floatは4byteとしてざっと8MBの領域をこの配列が占めることになります。そこに書き込みを行うので、破綻してコアダンプという事態にあいなりました。

(どうすればよいか、は聞かれていないので、敢えて説明しません。自分で調べてください。'ローカル変数　サイズ'あたりをキーワードに検索すれば結構なことがわかると思います。もちろん「スタック」を知らなければそれもC言語をあつか必須の知識ですので調べておきましょう。)

Comments

[thkana]

> どうすればよいか、は聞かれていない
あ、本文にはなかったけどタイトルが「コアダンプをなくしたい」だったか...
まぁ、調べてみてください。

Answer 3

float img_att[nz][ny][nx]; で img_att のサイズは8メガバイトです。
Linux では、スタックサイズが 8メガバイトなので、すべての変数を
自動変数として割り付けることができません。
static float img_att[nz][ny][nx]; として、静的に領域を確保しましょう。

Answer 4

配列の添え字が範囲外になっているのではないかと思います。
配列を添え字でアクセスする直前で、添え字の値が範囲内かどうかチェックしましょう。