done boostバージョンの平滑化プログラムをtrunkにコミット svnのhelpでignoreが見つけられないんだけどこれってtortoise svnの独自機能だったのカナ？ smoothingの挙動を確認。大丈夫そう。 todo 内部の行列演算を帯行列用に最適化。 ちょっと行列の数値計…

帯行列の逆行列は一般には密行列なので、帯行列の逆行列の計算は。しかし、帯行列のLU分解はで出来るので、帯行列とベクトルに対してはで計算できる。これを実装。 /*! g++ main.cpp -g -DNDEBUG */ #include <boost/numeric/ublas/banded.hpp> #include <boost/numeric/ublas/io.hpp> #include <boost/numeric/ublas/operation.hpp> #include <boost/numeric/ublas/lu.hpp> #include </boost/numeric/ublas/lu.hpp></boost/numeric/ublas/operation.hpp></boost/numeric/ublas/io.hpp></boost/numeric/ublas/banded.hpp>