コンパイルはできるが、実行ができない(エラーメッセージ無し)

Question

前提・実現したいこと

多腕バンディット問題をc言語で扱うプログラムをLinux環境でコンパイル・実行しようとしています。
エディタはemacsでコンパイラはgccを使っています。
エラーメッセージがないため解決方法が分からず質問させていただきました。
どなたかご回答お願いします。m(__)m

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

コンパイルまではできるのですが、実行するとエラーメッセージは何も出ず、そのまま止まってしまいます。
処理に時間がかかっているのかと思い待ちましたが全く進みません。

該当のソースコード

#include<stdio.h>
#include<time.h>
#define M (5)
#define N (500)
#define L (100)

int ID[L][N];
int Reward[L][N];
double ave_R[N];
double CDR[N];
double mu[L][M];
double R[L][M];
double T[L][M];

/* Mersenne Twister */
void init_genrand64(unsigned long long seed); // initialzes mt
long long genrand64_int63(void); // generates a random number on [0, 2^63-1]-interval
double genrand64_real1(void);

int mySlotMachine( int id, double *p ){
  if(p[id] > genrand64_real1())
    return 1;
  else return 0;
}

int main(void){
  int n, l, k, s, i, a;
  double r, e, E, max;

  unsigned long long seed;
  double p[M+1]={ -1, 0.3, 0.1, 0.9, 0.7, 0.5 }; // 確率は固定．p[0] は欠番
  seed = (unsigned long long)time(NULL);
  init_genrand64( seed );

  for(e = 1; e < 11; e++){
    a = e;
    if(e == 0)
      E = 0;
    else
      E = e / 10;

    for(l = 0; l < L; l++){
      
      for(n = 0; n < N; n++){
	r = genrand64_real1();
	max = 0.0;

	if(r < E)//探索フェーズ
	  {
	    k = genrand64_int63() % M;
	    ID[l][n] = k + 1;
	    T[l][k] += 1;
	    Reward[l][n] = mySlotMachine(k, p);//スロットの結果
	    R[l][k] += Reward[l][n];
	    mu[l][k] = R[l][k] / T[l][k];
	  }
	else //活用フェーズ
	  {
	    s = 0;
	    for(i = 0; i < M; i++){
	      if(max < mu[l][i]){
		s = i;
		max = mu[l][i];
	      }
	    }

	    ID[l][n] = s + 1;
	    T[l][s] += 1;
	    Reward[l][n] = mySlotMachine(s, p);//スロットの結果
	    R[l][s] += Reward[l][n];
	    mu[l][s] = R[l][s] / T[l][s];
	      }
      }
    }
    double sum[N];
    //平均報酬
    sum[0] = 0;
    for(l = 0; l < L; l++)
      sum[0] += Reward[l][0];

    for(n = 1; n < N; n++){
      sum[n] = sum[n-1];
      for(l = 0; l < L; l++){
	sum[n] += Reward[l][n];
      }
    }

    for(n = 0; n < N; n++)
      ave_R[n] = sum[n] / (L*(n+1));

    //平均正答率
    for(n = 0; n < N; i++){
      sum[n] = 0;
      for(l = 0; l < L; l++){
	if(ID[l][n] == 3)
	  sum[n] += 1;
      }
    }

    for(n = 0; n < N; n++)
      CDR[n] = sum[n] / L;

    printf("%f %f %f\n",E, ave_R[N-1], CDR[N-1]);
  }


  return 0;
}

試したこと

https://teratail.com/questions/296945
この質問と状況が似ていたため、ローカル変数をグローバル変数に変えました、が上手くいきませんでした。

補足情報

Ubuntu18.04 LTSでLinux環境を使い、エディタはemacsコンパイルにはgccコンパイラを使用しています。
コードの最初の方にある/* Mersenne Twister */の部分には配布された別のプログラムが必要で、それも一緒にコンパイルしています。別のプログラムと一緒にそのプログラムをコンパイルしたときには問題なく実行までできていたのでそのプログラムには問題がないと思われます。

Answer 1

実行できないというよりダンマリってことですね。

デバッグ環境をすぐに構築できないなら、要所要所にprintf文でも入れて、

• どのprintf文で止まっているか？
• 特定のprintf文が延々と出力されるか？

で、ある程度箇所は特定できると思います。
一部のループ処理はブレスがないので、そういうループ処理内にprintf入れる場合は全体をブレスで囲みましょう。
自分はループ処理が１センテンスでもブレスで囲むよう習慣付いていますが、、、

Comments

[kentohno]

printf文を入れたところ、関数mySlotMachine( int id, double *p )の中に入れた際に延々と同じ文字列が出力されました。
しかしmain文の先頭にprintf文を入れてもまたダンマリしてしまいました。
これはこれで理解ができません。。
main文の中で関数mySlotMachineよりも先にprintf文を入れてもそれを無視して関数mySlotMachineがfor文のような動作をしてしまうのはどこに原因があるんでしょう、、

[DreamTheater]

スクロールアウトしたんじゃないですか？
因みにループは最大で５５万回ループするようです。
単に時間が掛かっているだけでは？

[kentohno]

よく確認したらそもそもその関数をmain文の中で何万回と(意図的に)ループさせていたので当たり前の挙動でした。後半の方のfor文でループごとに値を増やす操作で間違っていたみたいです。
解決できました、ありがとうございました。

Answer 2

・統合開発環境で、デバッガを使う
・要所にprintf/fprintfを入れる
・頭で処理を追う
のどれかです。

Answer 3

C言語のコードを書くなら、デバッグできる環境を整えましょう。
Eclipseや、WindowsならVisualStudioなど。
コードの任意の場所で実行を止め、変数のナカミを見ることができます。そこから１行づつ実行して、コードの流れを見れるようになります
そうすれば、アテズッポでコードを書かなくて済むようになります。