休養

スカイレイがまた放牧に出てしまった。
1月の京都2400M、睦月賞あたりを
目指しての短期放牧だといいんだけど。

GUI作りもそこそこに、最近は
コード書きを少し進めている。
圧力ベース解法（陰解法）。
これを使うと固体、気体、液体が
同時に解ける。これなら
使ってみたいという人も増えるかなー。

ライアン産駒好調

10月、11月と不調が続いていた
ライアン産駒が最近好調。

なかでも土曜の葉牡丹賞を
快勝したマイネルアナハイムには
来年のクラシックに是非出て欲しい。
次走はホープフルステークスが有力らしい。

見事勝利！！

今日の6R、東京2400Mを見事勝ち上がってくれました！！
今日は一日幸せな気分です。その後のレースは純粋に
楽しめました。

それにしても強い勝ち方でした。スタートも何事もなく
出てくれて中段やや後方。直線入り口では一瞬
前が壁になるようなシーンもありましたが、
残り400Mくらいから追い出し。見る見る先団に追いついて
最後は2馬身1/2突き放す完勝。上がりは34.5、
走破時計も2分26秒7となかなかです。

次走は12/23　阪神　1000万下　2400Ｍの江坂特別
少し間があくけど1/28　京都　1000万下　2400Ｍの睦月賞
あたりでしょうか。また東上って可能性もありますが。
それにしてもこれからも注目です。

出走

日曜東京6R500万下2400Mに出走決定。

メジロスカイレイ　栗坂重　53.7　39.7　26.4　13.5　馬なり

調子はよさそう。
騎手は…ルメール。
悪くはない、むしろいい。
でも武豊に戻るとおもっていたから。
ルメールは短期免許だから恐らく次は
乗れないし、武豊とは2戦2勝と相性も
よさそうだったから。まあ出遅れなくゲートを
出てくれて、無事に回ってきてくれて、
ついでに勝ってくれればそれが最高。

GUI作りはだいぶ形になってきた。
もうそろそろプログラムの配布の形式を
考えないと。最近はあまり使われていない
らしいけど静的共有ライブラリにしようかな。
データ出力部分のソースをユーザー側で
編集できるようにして、あわせて実行ファイルに
するようにすればいいかな。
まだ時間はかかりそうだけど。

順調

今週の東京に出走するのでしょうか。
順調に調教をこなしているようです。

メジロスカイレイ　栗坂路　57.3　41.2　27.5　14.0　馬なり

JCウィークなので武豊騎手も東京。
武騎手に戻るのか、岩田騎手のままなのか。
それとも別の騎手？もうあんな出遅れだけは
勘弁してくれ。

次走に向けて

スカイレイが既に帰厩していました。
さっそく時計も出しています。

メジロスカイレイ　DW良　80.9　65.7　52.1　38.7　－　12.3　一杯

だそうです。ロックスピリッツ（一杯）と併せて
0.2秒先着とのこと。なかなかいいじゃないでしょうか。
来週日曜、東京500万下2400M？
これを逃しても中京の特別レースに
2500M戦がいくつかありますが。
どうするんだろう。

昨日スカイレイとおなじブライト産駒の
マキハタサイボーグが1000万下の
特別レースを勝利！！

ジュビロも快勝！！
リーグ戦最後のヤマスタだし
行こうか迷っていたのです。
あー残念。12/9の天皇杯5回戦
行こうかなー。

ＧＵＩ作り、なかなか楽しいです。
工夫のしがいがあって。
購入した本を辞書変わりに
（といっても索引が少ないので
不便ではあるのですが）
使ってます。その分、流体計算コード
の方が書き進んでいません。
バランスをとりながら進められるといいんですが。

スカイレイ休養…

せっかっく約1年ぶりに出走したのに、
先月31日から休養に出たとのこと。
これで今週末の出走もなさそう。
長期休養にならないことを望む。

少しずつGUI作りも進行中。
数値計算とはまた違った楽しさが
あっていい。

トライ＆エラー

WideStudioで簡単なアプリケーションを
作り始めてみた。しかしそこに出てくる
文法がさっぱり分からない。
ネットで調べたら構造体やらなんやら。
->な記号にも初めて出あった。新鮮な感じ。

ということで2つの方向から進めることに。
一つは、当然ながらちゃんと文法を学ぶ。
もうひとつはトライ＆エラー。
WideStudioにのってるサンプルプログラムと
勉強した文法から、自分のやりたい事をやるには
こう書いたらいいのかな？と論理的に考え実証してみる。

間違えることで文法を正しく理解するし。
ということで学んだ事の防忘録。


......

FILE *fp, *fp2; \\fp, fp2はポインタ？
double para1, para2, para3;
fp = fopen("filename", "r"); \\r: 読み込み専用。
if(fp==NULL){
exit(-1); \\これじゃ強制終了
\\⇒新しいウィンドウ出して、警告を与えるようにしないと。
}
fscanf(fp, "%lf %lf", ¶1, &para2); \\%lfは倍精度実数。
\\ 一行に2つのパラメータ。その間は半角スペースの場合。
\\ポインタ指定で読み込み。
fclose(fp);

fp2 = fopen("filename2", "w"); \\w: 書き込み用
if(fp2==NULL){
exit(-1);
}

fprintf(fp2, "%lf %lf \n", para1, para2); \\今度はポインタ指定ではない。
\\読み込み時は改行(\n)必要なかったけど。。。
fprintf(fp2, "%f", para3) ;
fclose(fp2);
......

といった感じ。間違ってないよね？
よくやってしまうのが ; を忘れる間違い。
やっと慣れてきた、というより最初にこれを疑う。
大概これでエラーが出ている。
流体計算コードも書き進めないとな。

WideStudio

WideStudioを使って本格的にGUIを
作り始めようと思い、
「WideStudio 徹底ガイドブック」
「C言語によるUnix システムプログラミング」
という本をAmazonにて購入。

前者を読み始め、簡単なGUIを作り始めました。
まだまだやりたいことをするにはどうしたらよいか
さっぱりな状況ですが、頑張ってみようかと思ってます。

スカイレイ出走

約10ヶ月ぶりの出走。
直前の調教内容は

メジロスカイレイ　栗坂路　　　　　　　　　55.4　41.1　27.2　13.6　馬なり
メジロスカイレイ　DW良　83.9　68.9　53.5　40.0　　―　　12.0　強め
メジロスカイレイ　栗坂路　　　　　　　　　59.5 　53.7　28.8　14.3　馬なり
メジロスカイレイ　DW稍　84.8　68.6　54.1　40.3　　―　　12.0　仕掛け

そして今日、東京8R500万下芝2300Mに
岩田騎手で出走しました。

結果は…
１枠１番のためゲートに先入れ。
スタート時に立ち上がり、さらに騎手を振り落とすほど
ジャンプする感じでゲートを出る。大きな出遅れ…
正直あれは痛い。道中は最後方を追走し、
3コーナー過ぎから4コーナーあたりで
中段にとりつき、もしかしたら、と思わせましたが
3着まで。スタート直後に足を使ったと思われ、
それでも3着まで来たのだから弱くはないと思う。
それにしても残念で仕方ない。

これで6戦2勝。3着が2回。4歳でまだ6戦ですから
まだまだこれからなはず。次走は去年勝った
11月18日（土）、京都8R500万下2400Mかな？

競馬あれこれＩＩ

今週の出走はありませんでした。

メジロスカイレイ　助手　ＣＷ重　81.6　66.4　52.7　39.5　－　13.4　一杯

木曜の追い切りの動きがあまりよくなかったですし、
評判の良血とはいえ2歳馬に遅れてしまったので
まだ本調子ではないのでしょう。
先週のアレグレットのこともありますし、じっくり
行った方がいいでしょう。

そして今日の追い切り

メジロスカイレイ　助手　栗坂路　55.4　41.1　27.2　13.6　馬なり

京都2400Ｍは来月中旬までないけど、それまで
じっくりかなー。やっぱり気がはやる。

プログラムの方は一応熱伝導も入れてみた。
でも熱伝導は研究で扱ったことがないので、
よいテスト問題を知らない。元指導教官に
聞いてみよう。

競馬あれこれ

メジロスカイレイの未勝利戦は
2400Mじゃなくて2200Mだった。
来週の土曜の京都500万下2400Mに
出るかな？

メジロスカイレイ　助手　CW重　69.7　54.2　39.7　　－　　12.6　馬なり
メジロスカイレイ　助手　栗坂路　　　　　57.3　41.3　27.2　13.8　馬なり

順調に調教をこなしているようです。
後者の調教は噂の良血馬、
ザサンデーフサイチと併せ馬。
馬なりで併入だそうです。

今日はメジロアレグレットが競争中止で予後不良…
今日はグリーンチャンネルの栄光の名馬で
メジロドーベルが取り上げられる日。その日にまさか
娘が亡くなるとは。信じられん。

悪いことばかり言ってると暗くなるんでいい話も。
メジロブライトと3/4の同血、メジロベリンダが入厩。

あと知らなかったのだが、メジロブライトの全妹
メジロプレストが繁殖にあがっているんだね。
ブライト死んじゃったし、貴重な血。
若いからまだ種付けはしないのかな。

久々更新

更新が滞ってしまった。
勤務地変更があり、引越しやらなんやら。
車通勤にも大分慣れました。
でも慣れすぎには注意です。

計算コード作りですが、
１次元のMHDコードが最近完成。
近似リーマン解法の１種であるRoe法に
空間精度を高次精度化するTVDの１種の
MUSCL内挿を組み合わせました。
時間方向に関しては、２次のルンゲクッタ法を
使用。全体で空間２次時間２次の精度です。
次は熱伝導を加えようかと思ってます。
さらにはソルバーをCIP法かRoe法かを
選択出来るようにして、オープンにしようかと
思ってます。いつになるかは分かりませんが…。

競馬では明日の夜、凱旋門賞です。
嫌いなサンデー産駒ですが、
ディープは応援したくなります。
特にこの凱旋門賞は勝って欲しい。
日本競馬のレベルの高さを示して欲しい。

自分の好きな馬、メジロスカイレイが久々に帰厩。
さっそく時計を出してますよ。

メジロスカイレイ　助手　栗坂路 　　　　　　　58.6 42.9 14.1 馬なり
メジロスカイレイ　助手　DW良　85.7 69.0 53.7 40.3 13.0 馬なり

だそうです。500万下からの再スタートです。
２戦２勝の京都2400Mのレースですが、
10/14（土）　8R　500万下　2400M
11/18（土）　8R　500万下　2400M
　　　　　　　10R　八瀬特別（1000万下）　2400M
ですかね。とにかく無事に行って欲しい。

細々と修正

Freeの某Fortranコンパイラーをインストール。
ＧＵＩソフトを作るプログラムをインストール。

これで準備完了と思いきや、
Vine2.5ではＧＵＩソフトが上手く出来ない
（ビルドオールしようとすると、make.exeがstat出来ない、
makefileがねーよ、ってエラーが出る）。
検索した結果、OSをVine2.6にして、
GUIソフト作成プログラムのヴァージョンを上げろ、
とのことで、Vine2.5→Vine2.6に。
上手く行くことを確認したんで今日はここまで。

中田引退。あー。ひとつの時代の終了って感じ。
やっぱりフランス大会の代表が1番好き。
名波、中田、山口（素）の中盤を忘れない。

VMware Player

前回書いたGUIで動かせるCFDソフトの製作を
前進させるため、手始めに開発環境を整えるという事で、
VMware PlayerなるVMwareのフリー版（？）を
インストールすることにした。

http://blog.yasaka.com/archives/2005/10/vmware_playervi.html
http://d.hatena.ne.jp/ISISI/20060211/1139665963
http://matanet.ath.cx/vmwareplayer/
http://animemo.seesaa.net/category/1349221.html

ここらへんのページを参考にさせてもらいました。
ただし1番上のページにある"Browser Appliance"を使った
方法をそのまま試しても、希望したOS（私の場合はVine
です）をインストール出来ず、勝手にデフォルトのOSが
インストールされてしまったので、それ以外のページある
方法でインストールしました。

何とかVine3.2をインストールしたんですが、
今度はGUI作成ソフトであるWideStudio（以下、WS）の
インストールがうまく行かず…。
WSのHPからrpmファイルをdownloadして、
rpm -ivh ***.rpmとしても依存関係の欠如（だったっけ？）で
インストール出来ず。C++のライブラリが足りないみたい
なんです。色々試しても上手く行かず、rpmファイルが
Vine2.5用らしいのでVine2.5をインストールし直して
しまいました。

上手くいったんですが、今度はインストール終了後の
rebootでF2（であってる？）を押し忘れ。
bios設定のboot優先順位１位をCD-ROMにしたまま
2回も起動してしまう。インストール作業が始まって
しまって、どうやってrebootしていいか分からず…。
VMware Playerでどうやって仮想マシンの電源を落とすか、
どなたか教えて下さーい。
どこかにそんなボタンありましたっけ？

ともかくも何度もインストール作業を繰り返し
（1回1時間超!）、やっとVine2.5をインストール
出来ました。
そしてWSのインストール。
WideStudio-3.20-1vl0.i386.rpm
ではなぜかインストール出来なかったが、
WideStudio-3.20-1vl2.i386.rpm
では出来ました。

しかししかし、今度は起動時にlibws.so.3.20が
見つからないというエラー。
また色々と検索してみて、
/etc/ld.so.configに
/usr/local/lib
と一行追加し、
/sbin/ldconfig コマンドを実行
することで解決することが判明。

まだまだ環境を整える段階で先は長そう。

さて昨日は大阪で「ダビンチコード」を見ました。
正直、キリスト教関係の用語が多すぎて置いてけぼり。
でも謎解き部分は結構楽しめました。もう一度じっくり見て
頭の中を整理したいっす。
あとはバーゲンセールということで、服を2着。
しめて2万6千円。予算をオーバー。
そんなに買うつもりはなかったんだけどな。

Fortran90 メモ（疑問）

仕事で市販のＣＦＤソフトを使ってるんですが、
かなり凄いです、当然ですが
（ライセンス料が年間1千万くらい？）。

それに触発されて、ＧＵＩで全て操作出来る
（可視化までは難しいかなー）流体計算プログラムを
作ってみようかと思ってます。かなり長期的な目標ですが。
メッシュ生成まではやらないかな。

それでFortranメモ（疑問）をつらつら書いていこうかと
思ってます。

double precision, allocatable :: w(1:100)
は良くても
double precision, dimension :: w(1:100)
real*8,dimension :: w(1:100)
は駄目？コンパイルでエラーが出ちゃいます。

Fortran90で、Fortran77みたいに配列を仮引数として
受け渡すには、interface文を利用。

例:
受け渡し側で
interface
subroutine test(w,n)
double precision, dimension(:), intent(in) :: w
integer :: n
end subroutine test
end interface
.....
call test(de,100) !interface文では適当は配列wを使って、
call test(pr,100) !subroutineの機能を定義。

受け取り側で
subroutine test(ww,nn)
double precision,dimension(1:nn), intent(in) :: ww
integer :: nn

といった感じ。
間違いがあったら是非指摘してください。

皐月賞馬はやっぱり強い…

◎アドマイヤムーン　 ⇒7着
〇メイショウサムソン ⇒1着
▲サクラメガワンダー ⇒10着
△アドマイヤメイン　 ⇒2着
×マルカシェンク　　 ⇒4着

結果はこんな感じ｡
武豊めずらしく内に突っ込んで前が塞がる。
ただそれを考慮しても離され過ぎ。
共同通信杯まではダービーの1番候補だったけど…。
まあ今年は、これという応援馬がいなかったから
どうでもいいや。

好みだけの日本ダービー予想

発走まであと1分！
◎アドマイヤムーン
〇メイショウサムソン
▲サクラメガワンダー
△アドマイヤメイン
×マルカシェンク

またまた

久々の更新となってしまいました。
引越しの関係でネットにつなげなかったのですが、
先週やっとつながりました。これからも時間を
見つけて更新するつもりですのでよろしくです。

久々に更新

ここのところ、引越し準備などで
更新は停滞気味。

順を追って振り返ると…
20日：追いコン＠聖護院御殿荘

21日：荷物だけ引越し。朝から細々としたものを荷造り。
　　　　WBCの決勝を見つつ、引越し業者を待つも、
　　　　フリー便ということでいっこうに来ない。結局来た
　　　　のは6時過ぎ。この日から友人宅にお世話になった。
　　　　ちなみにジュビロは今季初勝利。

23日：学位授与式。博士になりました。
　　　　だからと言って、何も変わりませんが。

24日：家電引き取り。午前に来るという話だったのに、
　　　　まったく来ない。電話したらしんどくなって休んで
　　　　いたとのこと。何ー！？さ、最悪だー。

25日：夜にいつも夕飯を食っていたメンバーで鍋パーティー。
　　　　最初は焼肉したり、普通の醤油ベースな鍋だった
　　　　のに、チーズを入れた辺りからおかしくなり、ヤクルト
　　　　やら蒟蒻畑（グレープ味）が入り、ちょっと変な
　　　　テンション。食べられなくはないんだけど…。
　　　　臭いは食べ物のそれではない（牛小屋の臭いに
　　　　近い。チーズの臭いが効いているんでしょう）。
　　　　デザートうどん（グレープ味）にはさすがに手が
　　　　出なかった。でもホット蒟蒻畑は食べたよー。

26日：京都ともお別れ。部屋の片付けを済まして、
　　　　大学を出たのが4時過ぎ。重い荷物を持って
　　　　京都駅へ。指定席は取れず、並んで列の先頭を
　　　　確保。18時発のひかりで、何とか座れた。
　　　　すぐに就寝。家では荷物が多すぎて自分の部屋には
　　　　寝れず、居間で寝る生活。今週は実家（かな？）

ジュビロvsサンガ

昨日は久々にジュビロの試合を見れた。
京都テレビに感謝。うちでは映りが非常に
悪いのだが何とか誰かは判別出来た。

見れたのはいいのだが、非常にストレスが
溜まる試合だった。その原因は監督の
采配にある。
1. 4バックの両サイドが守備的過ぎる。
　 右に鈴木、左に服部。意外と（と言ったら失礼
　 かもしれないが）鈴木は攻撃参加をしていたのだが、
　 左の服部は少ない。特にサイドを駆け上がらない、
　 １対１で勝負して抜いていくことがない。
　 １点目の時は村井と絡んでいい仕事をしたけど。
　 リードしての終盤になったときの選択肢としてはありだと
　 思うが。右に菊地（太田）、左に村井という攻撃的かつ
　 若い人選を望む。
2. 名波をボランチへ下げるべき。
　 やはりプレスのきついトップ下ではきつい。よさが
　 生きていない気がする。福西が怪我気味だし、
　 名波、ファブリシオ、福西、菊地、船谷あたりで
　 調子いい選手をボランチに。
3. 成岡を使え。個人的希望。

全体的に運動量が足りない感じがした。動けてない。
だからこそ若手の起用が必要だと思うんだけど。
まずは監督を変えて。これまた個人的希望。
なんか好かん。もっと攻撃的にいける監督希望。
結果も出せていないし。去年から何の進歩もないじゃん。
あーむしゃくしゃする。

話は変わって、金、土と二日連続で
「女王の教室」のスペシャルがあった。
続編希望には反対だったのだが、
実際は過去の話だった。意外とよかったと思う。

さらに、TUTAYAの無料レンタル券で
リベリオンというマトリックスっぽい映画を
借りた。どこかのサイトで評判がよかったから。
期待しすぎたせいか、それほどでもないと
感じたが、まあ楽しめることは楽しめる。
新鮮な（というか、ちょっと笑える）のが、
拳銃で殴ったりするガン=カタという戦法。
はてなダイアリーにある解説（もとはパンフらしい）
によると

「銃撃戦の統計学的分析の元、敵の死角に身を置き、
瞬時に相手を倒すそれは、剣道、空手を始めとする
東洋武術と銃の技術を融合させた無敵の「型」である。

銃を構える多くの敵の前に平然と姿を晒し、相手の弾丸を
ギリギリでかわしつつ、たちどころに敵の生命を絶つ。」

だそうだ。はっきりいって主人公強すぎ。さらにかっこいい。
何も考えずただスタイリッシュなアクションを楽しみたい人に
お勧め。

レッズvsジュビロ

土曜日の試合、関東ではテレビ中継が
あったのに、こちらでは阪神のオープン戦…。
なんで…。

前々から希望していた4バックだったから
期待していたが、機能しなかったみたい。
サイドバックが鈴木、服部じゃあね。
スポーツニュースでは三都主に
ぶっちぎられるシーンを何度も見た。
やはりベテランはベンチスタートに
して、若手を起用してほしい。
同じ負けるならね。

左SBに村井、右SBには菊地、もしくは
太田なんてどう？攻撃的過ぎるかな。
CBは田中（もしくは鈴木）、金でいいかな。
ボランチには福西（ファブリシオ）、名波（菊地、船谷）、
攻撃的MFに左から成岡、前田、西（カレン）、
１トップにはカレン（西野）、これでどうだ。
前田が復帰するまではトップ下に成岡、
左にカレンとかかな。
その前に山本監督解任が必要かな。

今週の土曜日は京都戦。京都には
悪いがジュビロ勝て。今の京都
（川崎Fに7-2で敗戦）に勝てないようじゃ
どこに勝つんだ。さらにホーム（ヤマハスタジアム）。
負けたら暴動が起きるか？

京都テレビが録画中継してくれるんで、
久々にテレビで見れる。京都のホーム
だったら西京極に行ったんだけどな。

写真は本文にはまったく関係ありません。
平等院の近くの橋から撮った宇治川の
写真です。

免許取得

昨日は6時起きして、免許試験場へ行ってきました。
京阪で中書島へ、バスで試験場へ。かなりの
人数がバス停の近くで待っていました。
さらに試験場につくと凄い人数が並んでいた。
毎日のようにあんな人数が来るんでしょうか。

試験は思ったより難しく、これは落ちたかなと
心配しました。しかしなんとか合格できました。
合格発表時には歓声が…。そういうものなのかな？

合格発表から免許受け取りまで2時間以上待たされた。
やはりあの人数じゃあしかたないのか。
iPodを持っていくのを忘れたんで、何もすることなく、
コンビニに行って立ち読みとかで時間をつぶす。

帰りはバスの時間まで1時間近くあったんで、
簡略な地図を頼りに中書島駅まで歩いてみた。
所要時間は45分くらいだったかな？
迷うことなく辿り着けた。迷ったら迷ったで
それなりに楽しくなったのだろうが。

夜は追いコン。花束をもらってしまいました。
本当に有難うございました。しかし1次会で
食べ過ぎたのか、かなり気持ち悪く早起きの
せいで眠いし。ということで2次会はおとなしく
してました。3次会ではとあるビデオを
見ながらのゼミがあったようだけど、限界を
迎えていたので帰って寝てしまった。

宇宙戦争【ネタバレあり】

教習所の待合所で何度か見かけていて、
結局最後まで見ていなかったので、
教習所卒業記念（？）ということでDVDを
レンタルして見てみた。

あまり気にしていなかったが前半1時間弱くらいは
教習所で見てしまっていた。ということで未見部分は
後半1時間くらい。

個人的感想としては、1番面白い部分を教習所で
見てしまっていたという感じ。宇宙人の操縦する
マシーン登場から大量殺戮、それからの逃亡シーンは
非常によかったと思う。だからこそ続きが気になって
レンタルしたわけだが。

結局この映画、宇宙人出現を物語の始まりとしているが、
描いているのは家族とか人間の本性とかそういう部分が
メインな気がしている。ありがちだが、本当に恐ろしいのは
人間、と感じさせる描写もあったし。実際宇宙人は我々の
預かり知らない所で勝手に滅びてしまうし。
未見部分だけの満足度は低いかな。

卒業検定

去年の9月に通い始め、D論を書きつつ、
ちまちま通った教習も今日が大詰め。
午前中に卒業検定を受けてきました。
日曜日には送迎バスの無い所なので
自転車で30分程かけて教習所へ。

一般コースと自主経路、さらに
所内検定。目安としてそれぞれ
2.5km、2.0kmと書いてあったんで、
そんなに走るのか、と少しびびった。
でも自主経路は簡単な道で教習所に
帰ればいいとのことで、結局1.0kmくらい
になった。

順番は2番目。少し緊張していたが、
最初の人の運転を色々観察しつつ、
自分ならこうするな、とか考えていたら
落ち着いた。一般コースが簡単なコース
だったのもラッキーだった。心配していた
駐車（練習では寄せすぎてタイヤをすっていた）も
うまくいったし、所内での方向転換も問題なく
うまく行った。ということで無事合格。
あとは学科試験。あと少しだー。

Fortran90コードが

物理的に怪しい挙動を示す。
ということで、またFortran77コードに戻った。
やはりFortran77で書き始めたコードを
途中でFortran90に書き換えたからどっかに
バグがあるんでしょう。恐らくcommonの
関係（配列の共有）だと思う。

先日発見したバグを修正したら
以前から気になっていた点も改善されたので、
暫くは（ずーっとでもいいけど）Fortran77コードで
行く事にしよう。

今週は

火曜日は大阪の水族館（海遊館）＋アウトレットへ、

水曜日には中学時代の親友が京都へ遊びに
来たので、七条油ノ小路（伊東甲子太郎が絶命した場所）、
三十三間堂へ、

木曜日はやはり親友と宇治（平等院近辺）、西本願寺、
東本願寺（両本願寺とも、平成の大改修をしていて、
かなり残念。特に西本願寺は新撰組の屯所に使われていた
から楽しみだったんだが）、京都タワーへ、

金曜日は、親友と一緒に観光行きたかったが、午後に
予定が入り断念。しかし時間が決まっておらず、結局
夕方になったので観光に行けた…。
しかしもうそろそろ健康診断を受けておかなければ、
会社へ提出する期限が迫っているのでヤバイと思い
下宿近くの病院へ。すごく近いのに初めて行った。

こんな感じの1週間だった。そのうち写真をアップしたいと
思っています。いつになるかは分からないけど…。

高速教習

昨日は、路上教習も終盤、高速教習へと
行ってきました。風邪の影響か、前日は腹痛が
酷く、高速教習へ行っても大丈夫かと心配
しましたが、なんとか当日は腹痛もましになった
ので行ってきました。幾分熱っぽかったのですが。
だってキャンセルしたら3時限分＋高速代700円を
もう一度払わなければいけないので、もしそう
なったら痛い出費です。

朝9時10分集合で、3人1組。高速までは自主経路。
私は3人中3番目の順番。なかなか順番が回ってこず、
少しドキドキしながら順番を待つことに。教習所で
知り合った人の話だと、2番目がいいということを
聞いていたのですが…。

2番目の人は、高速道路へ行く途中で変わり、
そのまま高速を乗って下りるところまでずーっと。
それで今日は終わりということらしい。つまり
一度も乗り変わりがないので確かに楽。
次は俺かと思ったら、最初の人が高速道路を運転。
まだかよーと思いつつ、緊張からトイレへ行きたく
なったり。やっと自分の順番。思っていたよりは
スピード出しているという感覚もなく、ETCを体験
したり、初めてセカンドギアを使ったりして終了。
一度も中断なく終了したんでまあ楽だった。
最初の人は、最後教習所へ帰る道も運転だったんで、
最初から最後まで気が休まることがなかっただろうな。

あとは山岳、自主経路、教習効果確認が終われば
卒業検定。頑張って朝早く行って予約を取らないと、
引越しまでに間に合わなくなってしまう。

応急救護処置

4月からの就職に備え、この歳（27）になってようやく
教習所に通いはじめ、今現在、2段階の終盤。
今日は学科1時限＋応急救護処置（1時限学科、
2時限実習）。気道確保、人工呼吸、心臓マッサージ
の3つが応急救護のABCだそうです。
いざとなったら出来なさそう。免許取った多くの
人がそれ以降復習することが無さそうだし。

All OK

コードを見直し、全ての原因が判明。
凄く単純な理由だった。
Fortran77→Fortran90の書き替えの動機となった
怪しい振る舞いも改善。これですっきりした。

今後何をしようかな。Cを勉強するか、さらに
コードを発展させるか。思い切って非構造格子に
よる数値計算の勉強でも始めるか。

原因判明

MPI部分が原因ではなく、それ以外のところが原因だった。
今自分が解いている問題に依存した事であることが判明。
具体的には、円柱座標系で解いているのだが、軸（z軸）上の
値を周りの点の平均として計算している。それにより
軸上でも軸方向成分以外のベクトル成分が0でない値を
持てるようにしている。その部分が悪さをしている模様
（今日はそこまでしか分からなかったが、大きな進歩）。
明日にでもゆっくり考えよう。

codeのコアな部分（並列化部分を含む）のバグじゃなかった
のには一安心。しかし逆にそういう細かい所のバグで、
ずーっと悩んでいたのは悔しい。Fortran77→Fortran90の
書き替えを決心した不可解な振る舞いもこれに関係ある
のか？まあどちらにしろFortran90に少し詳しくなったし、
よしとしよう。

並列化

MPI（Message Passing Interface）を使って、codeを
並列化してるんですが、どうもうまくいかない。
並列化とは１CPUでは出来ない計算を、多くのCPUを
使って計算しようということなんで、当然幾つCPUを
使っても結果が一致しないといけない。
でもそれが一致しない…。

以前のFortran90コードは一致したんだが、
今回のは一般座標系（直交曲線座標系）コードにしたのと、
時間精度をあげるための工夫をしたため、頻繁に
MPIによるデータ転送が発生。
ここのところ一生懸命原因究明を試みるものの、
まったくもって分かりません。もう嫌になってきました。
他の事したい。

話題は変わって。昨日は夜中帰って、NHKで
スピードスケートを見た。金を期待された加藤選手も
メダルには届かず。今回はまったく駄目？あと
期待できるのは女子フィギュアスケートと
女子スピードスケート、ジャンプくらい？
まったくメダルが取れなくても不思議ではない。

今週のサッカーマガジンはジュビロ特集。
名波、成岡のインタビューが読めるので
買った。ちなみにサッカーダイジェストでは
成岡のインタビューが載っていたが、立ち読みで
済ます。あと１月程でJリーグも開幕です。

Fortran77→Fortran90

最近は自作したMHD（MagnetoHydroDynamic）
codeをFortran77からFortran90に書き換える作業に
取り掛かってます。

何故かというと、Fortran77で出来た配列を仮引数として
受け渡す事がFortran90では出来ないと思っていたのです。
しかしInterface文というのを使うことによって、Fortran90でも
出来る事を知ったので。やはりコンピュータ言語を極める
のは難しい。分かったつもりで書いていたのですが、
全然分かっていなかったということ。学部の頃の教科書に
しっかりと書いてあったとは…。

話題は変わって。サッカー。
アメリカ戦、私の住んでる所では
放送がなくて見れませんでしたが、
ボロ負け？アメリカが強いのか、日本が弱いのか。
両方かな？日本はシーズンオフだし。
アメリカのメジャーリーグはどうなんだろう。
野球と同じでオフ？
もう最近はジュビロ＞代表なんでどうでも
いいけど、Jリーグのレベルを甘く見られるのは
嫌だし、やはり勝っては欲しい。なんだかんだで
W杯が近づけば応援するんだろうな。

blog移行

今までYahooでblogを書いてましたが、
あまりの重さに辟易。
ということで引っ越してきました。
とりあえず初投稿。
色々日々思ったことを書き留めて行く
予定です。