Repo
kjana00@gmail.com
2011-01-07 01:21:45
今朝は……くもりっていうには晴れてるけど、日が陰っててなっていう感じ。朝ごはんはパンとサラダとベーコンエッグ。二日分半端にソーセージ買っても、ね。とりあえずちゃんと目は覚めたし、あんまりひどく寒くもないし、初日の頭としてはまずまず?
それから年末に片付けるつもりだったコードいじりを延々やりつつ脇で持って行くことになるかもしれない共用 PC に Cygwin を突っ込むだけ突っ込んでおくという作業。本当に長期出張になるんだか、どうなんだか。個人のアカウント作った方がいいのかなとかもちょっと考えるところ。
で、いじってたコードの方は何とも思わしくない感じで……変な値になってたのは一部負号を間違えていたからだったっていうのを直したら、とある引数での結果は正しくなった。んで、引数を変えて曲線のデータにするなんてやりだしたらどうもまともに動かない。おかしな引っ掛かり方をしているような。それで結果を double で fwrite() する脇で printf() してやったら動いてる様子は見えたけど、明らかに結果がおかしいという。解は単調減少する関数なんだけどな、これ。
しばらくああでもない、こうでもないといじくり回してみる。特異値分解の関数を変えてみたり、linear equation solver だの generic equation solver だの使ってみたり。でもやっぱり変なのだよね。そもそも収束しないし。延々といじってるとこころでちょっと怪しいところがあって、確認してみたら gsl_matrix_view_vector() ってベクタを列方向に並べるのを優先するんだねっていうのがわかったんで、解のベクタには列ベクタを並べる形で入ってるんですと直してみたらちょっとましになった、ぐらいで時間切れ。ほんと、ちょっとましになったぐらい。収束するけど引数に依存しない値に収束するもんな……
とりあえず今日のところはあきらめて帰って一休み。それで晩ごはんはおでん。辛味大根はやっぱり大根っていうより芋っぽいよねと両方入れて思ったりする。ちょっと辛味大根の方、細かく刻み過ぎたな。食べ難い。
後は WWW 見たりゲームしたりと。ちょっと熱が出たっぽかったけど、晩ごはん食べてしばらくしたら安定したっぽい。久しぶりに出っ放しで仕事だったから疲れた?