ここでdxとdyはそれぞれx方向とy方向の格子間隔である。(T(i+1,j)+T(i-1,j)+T(i,j+1)+T(i,j-1))/4+heat_source*h^2/4.が得られる。ヤコビ法ではこの式が成り立つまで温度場Tを繰り返し計算する。heat_sourceが一定かつ一様で4という値であれば、右辺第二項は4*h^2/4となる。この項は一定なので、ヤコビ法の繰り返し計算で毎回計算するのは馬鹿らしい。doのようにheatという定数を足すだけで済むようにしている。コメントアウトされていない行では0.25を掛けている。これは割り算よりも掛け算の方が計算が速いからである。†thermal_diffusion.f90をコンパイルし、計算機「scalar」で実行せよ。今回は2次元データの可視化を行う。前回と同様gnuplotを利用する。正方形上(x,y平面上)に分布する温度をすべて書き出すものである。†thermal_diffusion.f90をmpif90コマンドでコンパイルし、ジョブスクリプトthermal_diffusion.jsで並列計算せよ。連番つきファイルtemp.2d.????が出力されるはずである。その中身は以下のようになっているはず。†ソースコードの変更。†thermal_diffusion02.f90で求めたチェッカーボード型熱源分布の平衡解(最終状態解)を可視化せよ。温度の最大値と最小値も記せ。今回の演習内容はどうでしたか?(どれか一つ、一度だけ押してください。「お名前」欄は空欄で可。最新の10件を表示しています。感想、質問など気楽に書き込んでください。名前欄を空白のままでOK。00:13:47これはテストです。