//program file name <mrproc.cpp>
#include <iostream>
#include <fstream>
#include <cstring>
#include <cstdlib>

using namespace std;

//大きな配列はグローバルとして宣言する。
	float	v1[1202][1202];		//ラスターデータ１
	float	v2[1202][1202];		//ラスターデータ２

int main(int argc, char *argv[]) {
	int	i, j;
	char	rn[128];		//地域名
	char	vn[128];		//変数名
	char	un[128];		//単位名
	int	dt;			//データのタイプ（未使用）
	int	xas, yas;		//要素数
	float	xs, ys;			//モデルの大きさ
	char	str[128];		//ヘッダーとデータを区切る文字列

	//２つのラスターデータは、すべてゼロで初期化しておく。
	for (j=0; j<=1201; j++) {
		for (i=0; i<=1201; i++) {
			v1[i][j] = 0;
			v2[i][j] = 0;
		}
	}

	//コマンドライン引数が１つでない場合は、使用法を表示して終了する。
	if (argc != 2) {
		cerr << "Usage: mrproc <filename>\n";
		return 1;
	}
	
	//ファイルが開けなかった場合は、その旨表示して終了する。
	ifstream fin(argv[1]);
	if (!fin) {
		cerr << "Cannot open file.\n";
		return 1;
	}

	//データの読み込み。
	fin.getline(rn, 80);
	fin.getline(vn, 80);
	fin.getline(un, 80);
	fin >> dt >> xas >> yas;
	fin >> xs >> ys >> str;
	for (j=1; j<=yas; j++) {
		for (i=1; i<=xas; i++) {
			fin >> v1[i][j];
		}
	}
	fin.close();

	//データの外挿
	for (j=1; j<=yas; j++) {
		v1[0][j] = v1[1][j] * 2 - v1[2][j];
		v1[xas + 1][j] = v1[xas][j] * 2 - v1[xas - 1][j];
	}

	for (i=1; i<=xas; i++) {
		v1[i][0] = v1[i][1] * 2 - v1[i][2];
		v1[i][yas + 1] = v1[i][yas] * 2 - v1[i][yas - 1];
	}

	v1[0][0] = v1[1][1] * 2 - v1[2][2];
	v1[xas + 1][0] = v1[xas][1] * 2 - v1[xas - 1][2];
	v1[0][yas + 1] = v1[1][yas] * 2 - v1[2][yas - 1];
	v1[xas + 1][yas + 1] = v1[xas][yas] * 2 - v1[xas - 1][yas - 1];

	//データの処理
	for (j=1; j<=yas; j++) {
		for (i=1; i<=xas; i++) {
			//ここに実行したい処理を記述する。
			//以下の例はラプラシアン演算子の適用。
			v2[i][j] = v1[i][j] * 4 - v1[i][j-1] - v1[i-1][j] - v1[i+1][j] - v1[i][j+1];
		}
	}

	//変数名に処理内容を適宜追加する。
	strcat(vn, "-Laplacian");

	//結果の出力。
	//標準出力（画面）に出力するので、ファイルに保存する場合は、リダイレクトする。
	cout << rn << "\n" << vn << "\n" << un << "\n" ;
	cout << dt << "\n" << xas << "\n" << yas << "\n";
	cout << xs << "\n" << ys << "\n" << str << "\n";
	for (j=1; j<=yas; j++) {
		for (i=1; i<=xas; i++) {
			cout << v2[i][j] << "\n";
		}
	}
	return 0;
}	//プログラムの終わり。