TCUP2次元コードのテスト

密度,圧力の変わりに圧力,温度を基本変数とした

TCUP法の2次元コードが完成しテスト中。

キャビティ内の自然対流の計算を行っている。

壁境界を扱うのは初めてなので新鮮。

でも左右の壁の温度差が2[K]弱の計算なので

なかなか対流が成長しない。この結果は合っているのだろうか。。。

物理時間2秒の計算をするのに、200万time-stepもかかってるし。

瑛太

瑛太見たさにアンフェア（ドラマ版）をレンタル。

黒幕のことは知っててみたけど、結構楽しめた。

余命一ヶ月の花嫁も見たいけど、さすがに男一人では行きづらいんで

レンタル始まったら借りて見よう。

今月見た映画

5/1：GOEMON
映画の日で安かったんで。キャシャーンでがっかりしてたので
保険をかけて、映画の日にしておいた。
キャシャーンよりは大分面白くなっている。
それでも、少し長いとか説教くさくなるとか不満は多少あり。
あと映像きれいだけど、目がチカチカする。

5/16：天使と悪魔
予告で非常に面白そうな雰囲気だったんで、初日に鑑賞。
この日の夜にダビンチコードがテレビ放送だったし。
CMで見て気になっていた、ミカエルが剣を鞘に収める像があるサンタンジェロ城。
映画を見て実際に行きたくなった。。。
現在ダビンチコードの原作を読んでいる。天使と悪魔の方が出版されたのは先なんだ。

5/23：スラムドッグ＄ミリオネア
ミリオネアのセットも音楽も、どこかの会社がセット販売してるんだ（wiki情報）。
アカデミー賞とっただけあって、確かによかった。
主人公を演じた子供、実際に住むとこないらしい。

浮動小数点演算の一貫性

数値計算をしていて奇妙な現象に遭遇。

実行ファイル、インプットファイルが

完全に一致している（diffで確認）のに、

実行ディレクトリが違うだけで結果が一致しない。

write文をいたるところに挿入して確認したところ、

最初は倍精度変数の一番下の桁だけ違う

（当然どんどん誤差は拡大していくが）。

Z = a*b + c*dの計算で

a,b,c,dは一致しているのにZが違う。

ここまで確認して、もうこれはどうしようもないかと諦めかけた。

今日色々と調べたところ、Intel Fortran Compilerであれば

コンパイルオプションに-fp-model preciseを加えることで

回避できることが判明。