読書メモ - シンメトリーの地図帳

book

シンメトリーの地図帳 (新潮クレスト・ブックス)作者: マーカス・デュ・ソートイ,冨永星出版社/メーカー: 新潮社発売日: 2010/02/26メディア: 単行本購入: 1人クリック: 17回この商品を含むブログ (24件) を見る群論の本。群論は、昔習ったし、その後も何度…

2011-01-30

ublasのlu_factorizeはbanded_matrixには使えないのか?

boost ublas c++

http://d.hatena.ne.jp/hotoku/20110121/1295538662なんて書いて喜んだのも束の間、ublasで帯行列とLU分解の組み合わせは一筋縄では行かぬかもしれない。 LU分解関数 lu_factorizeの中のpivotingの部分で、matrix_rowクラスのswap関数を呼ぶんだけど、こいつ…

2011-01-25

todo & done

作業記録

ini modeで#をコメントにすべし。引き続きlu分解デバッグ中。

2011-01-24

todo & done

作業記録

gimpをインストール中上手く行くかどうか。。。 db46のインストールに失敗 java のヘッダが無いと怒られている。遅いので明日。逆行列計算関数のデバッグ boost ublasのlu_factorizeでこけている。 BOOST_UBLAS_CHECK (singular != 0 || detail::expressio…

2011-01-23

todo & done

作業記録

lu_factorizeに入れると元の行列が壊れちゃうので、Q, Rは直接入れられない。ので、自動変数にコピーすることにする。サイズはそんなにでかくならないので、問題無かろう。後、Rをsymmetric matrixクラスからbanded matrixクラスに変更。 lu_factorizeでこ…

2011-01-21

今日やったこと & やるべき事

作業記録

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

2011-01-21

ublasのbanded_matrixとlu_factorizeで、帯行列の逆行列とベクトルの積を線形時間で

boost ublas

帯行列の逆行列は一般には密行列なので、帯行列の逆行列の計算は。しかし、帯行列の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>

2011-01-20

banded_matrixの積はちゃんと線形時間だった件

boost ublas

http://d.hatena.ne.jp/hotoku/20101201/1291219911の疑問が解決した。簡単なことだ。NDEBUGマクロを定義すれば良かったのだ。ublasは、NDEBUGが無い場合はデバッグビルドと看做してサイズチェック等のテストコードを埋め込む。なので、最大の速度を求めるな…