for文の中が途中から実行されない

Question

以下のコードを実行すると、
#t=0
のみ書かれたファイルが作られてしまいます。
エラーは表示されていないのですが、どうしたら配列hの値を計算、出力されるようになるのでしょうか。

該当のソースコード

#include<iostream>
#include<string>
#include<cstdlib>
#include<ctime>
#include<fstream>
using namespace std;

string name = "sand.data";
ofstream datafile;
double h[100][100];
double D = 0.1;
double Q = 0.1;
int L, L0 = 1, b = 1;
int i = 0, j = 0;

//boundary conditions
int i_bc(int i) {
	if (i > 100) return i - 100;
	else if (i < 0) return i + 100;
	else return i;
}

int j_bc(int j) {
	if (j > 100) return j - 100;
	else if (j < 0) return j + 100;
	else return j;
}


int main()
{
	datafile.open(name);

	srand(time(NULL));
	//yuragi
	for (i = 0; i < 100; i++) {
		for (j = 0; j < 100; j++) {
			h[i][j] = rand();
		}
	}

	for (int t = 0; t < 100; t++) {

		datafile << "#t=" << t << endl;

		//creep
		for (i = 0; i < 100; i++) {
			for (j = 0; j < 100; j++) {
				if (h[i][j] > h[i_bc(i + 1)][j_bc(j)] && h[i][j] > h[i_bc(i)][j_bc(j + 1)] && h[i][j] > h[i_bc(i)][j_bc(j - 1)] && h[i][j] > h[i_bc(i - 1)][j_bc(j)] && h[i][j] > h[i_bc(i + 1)][j_bc(j + 1)] && h[i][j] > h[i_bc(i - 1)][j_bc(j + 1)] && h[i][j] > h[i_bc(i - 1)][j_bc(j - 1)] && h[i][j] > h[i_bc(i + 1)][j_bc(j - 1)]) {
					h[i][j] -= D * h[i][j];
					h[i_bc(i + 1)][j_bc(j)] += D * h[i][j] / 6;
					h[i_bc(i)][j_bc(j + 1)] += D * h[i][j] / 6;
					h[i_bc(i - 1)][j_bc(j)] += D * h[i][j] / 6;
					h[i_bc(i)][j_bc(j - 1)] += D * h[i][j] / 6;
					h[i_bc(i + 1)][j_bc(j + 1)] += D * h[i][j] / 12;
					h[i_bc(i - 1)][j_bc(j + 1)] += D * h[i][j] / 12;
					h[i_bc(i - 1)][j_bc(j - 1)] += D * h[i][j] / 12;
					h[i_bc(i + 1)][j_bc(j - 1)] += D * h[i][j] / 12;
				}
				else continue;
			}
		}

		//saltation
		for (int k = 0; k < 10; k++) {
			int i_ran = rand() % 100;
			int j_ran = rand() % 100;
			if (h[i_ran][j_ran] > 0) {
				L = L0 + b * h[i_ran][j_ran];
				h[i_bc(i_ran)][j_bc(j_ran)] -= Q;
				h[i_bc(i_ran + L)][j_bc(j_ran)] += Q;
			}
			else continue;
		}

		for (i = 0; i < 100; i++) {
			for (j = 0; j < 100; j++) {
				datafile << i << " " << j << " " << h[i][j] << endl;
			}
			datafile << endl << endl;
		}

	}

	datafile.close();
	cout << "Output data in " << name << "." << endl;
	return 0;

}

Answer 1

Visual C++ でデバッグ実行してみました。
71行目で配列の境界を超えてアクセスしている為に落ちているようです。
Lが4209となっているのは想定された値でしょうか？

Comments

[24_k]

回答ありがとうございます。
初歩的なことで申し訳ないのですが、71行目で i_bc(i_ran+L) が100以上の値となってしまうため落ちてしまうということでしょうか。
Lの値は想定していたものではなく、Lの式を変えてみたり、65~72行目を削除して実行してみたりしたのですが、改善されていない状況です。

[cx20]

> 71行目で i_bc(i_ran+L) が100以上の値となってしまうため落ちてしまう
そこもそうですが、それよりも前に配列の境界を越えてアクセスしている箇所があるようです。

例えば、52行目の h[i_bc(i)][j_bc(j + 1)] は、
i が 0 、j が 99 の場合、i_bc(i) が 0 を、j_bc(j + 1) が 100 を返します。
配列は h[100][100] と宣言されていた場合はアクセスできるのは h[0～99][0～99] までです。
h[0][100] をアクセスした場合は、配列の境界を越えてアクセスしていることになります。

そのような場合は、予期しないメモリ領域を壊してしまい、アプリケーションが不安定になったり落ちたりします。

[24_k]

ありがとうございます。理解しました。
17~27行目を剰余の式に変えて、正常に値が出力されるようになりました。