Repo
kjana00@gmail.com
2010-12-09 00:44:05
今朝はくもり。朝ごはんはパンとサラダとソーセージ。卵が無かった……何かエネルギーが足らない感じ。単に昨日の晩から冷え込んできたから、そろそろ毛布でも出そうかっていうだけの話という説もあるけど。
係数をかけて差分を小さくするのは難しくても、単純に大き過ぎる項を制限するだけだったら簡単だなという話はあるんでちょっと試してみた。……あっさり上手く行った。制限がかかった時に動く方向は解に近付きはしないかもしれないけど、真っ当な計算が出来る範囲内に動くから後のステップで収束の軌道に乗れるっていう感じみたい。まあ少なくとも場外ホームランな感じに動いて無限大とか、行ったり来たり振動して結局何も起こらないとかよりはずっとましな動きになるっていうこと。
とりあえず昨日までで動いてなかった部分が動いたから続き。直列な列は計算出来たから、これを並列に並べましょう、と。それでまあ、悩む、悩む。概形は出来ても計算が変っていうことで。まあその計算が変っていうのもヤコビアン行列作り間違えたとか、方程式の中間値ベクトルの要素を並べる時に添字を間違えてたとか、挙句直列列のヤコビアン行列を並べる形にしてみてるけど、解ベクトルの要素の順番に縛られるんだから勝手にそんなわかり易い並べ方をしてちゃだめだったとかだったりするわけだけど。焦ったとか手を抜きたかったとか、まあもっと落ち着きなさいって。
最終的には概ね真っ当に計算出来るようになったっぽい仕上がり。u_n = u_n-1 - J(u)^{-1} . f(u) っていうのを素直にそのまま計算させてるから、ヤコビアン行列の要素にオーダー違いなところがあるとちょっと破綻気味っぽい辺りをどうにかしないといけないか。一応、逆行列を使わない方法も見つけてはあるけど、見つけた方法っていうのが概要だから中身を詰めなきゃいけない。
まあそれでも動いたことは動いたっていうことで昨日よりもいい気分で帰る。それで晩ごはんはカレー。サラダぐらい付けるべきだ、と思いつつもカレーのみ。目を離し過ぎて危うく焦がしかけたけど、まあ十分にいい出来でしょう……っていうか、カレーは普通悪い出来にならない。
後は WWW 見てゲームして。火力が足りないかと思いかけたけどそうでもない、という感じで迷宮のミノタウルスを倒して帰還のロッドを入手、即吸収して巻物を売り飛ばして荷物の空きを増やして……と。売り物のアーティファクトが結構良かったのも助かった。対混乱が付いたから、だいぶ安定感は増した感じ。