2011-02-01 todo & done 作業記録 lu_factorizeが出来るようになったと思ったら、lu_substituteでこけた。原因不明。気になるけど、放置。帯行列は諦めて、密行列として扱おう。どうせそんなにサイズでかくないんだし。時間掛かるようならまた考えるということで。 早過ぎる最適化は諸悪の根源、肝に銘じよ。 よく考えると、帯行列をpivoting有りでLU分解すると、出てくるのは帯行列でなくなるんではなかろうか。従って、計算時間がサイズに対して線形時間でなくなるんではなかろうか。要確認。