/**
 * ここはフレームワークから、
 * ここは1秒間に数十回、繰り返しコールされます。
 */
 void onUpdate(int renderer){

	// UPDATE_PER_FRAME回だけ繰り返し力学演算（時間発展）する
	for(int cycle=0; cycle<UPDATE_PER_FRAME; cycle++){

		timeCount++;
		if(timeCount == END_TIME_COUNT){
			alert("終了時刻です。より長い時間のシミュレーションを行うには、END_TIME_COUNT の値を大きく設定してください。");
		} else if(timeCount > END_TIME_COUNT){
			// 終了時刻以降は何しない。
			return;
		}

		// X / Y方向の線張力（※TENSIONは単位幅あたりの張力）
		float tensionX = TENSION * Y_LENGTH / Y_N;
		float tensionY = TENSION * X_LENGTH / X_N;

		// 波のタイプがサイン波の場合は、境界をsin関数で振動させる
		if(waveType == WAVE_TYPE_SIN){
			for(int i=0; i<Y_N; i++){
				vertex[i][0][Z] = SIN_HEIGHT * sin(timeCount * DT * SIN_SPEED);
			}
		}

		// 質点にはたらく力を計算
		if(FREE_END){
			updateForceFreeEnd(tensionX, tensionY);
		}else{
			updateForceFixedEnd(tensionX, tensionY);
		}

		// 位置と速度の計算
		for(int i=0; i<Y_N; i++){
			for(int j=0; j<X_N; j++){

				// 頂点の速度変化の計算
				vertexVZ[i][j] += vertexFZ[i][j] / mass[i][j] * DT;

				// 頂点の座標変化の計算
				vertex[i][j][Z] += vertexVZ[i][j] * DT;
			}
		}

		// モデルの頂点配列を更新
		setModelVertex(waveModel, vertex, QUADRANGLE_GRID);

	}
}