結局、ublasのlu_factorize関数周りの疑問はイマイチすっきりしないのでcholesky分解アルゴリズムを手で書いちゃうことに。あんまり難しいこと考えなければ、結構簡単なんだな。(トートロジーです) /*! g++ main.cpp */ #include <boost/numeric/ublas/matrix.hpp> #include <boost/numeric/ublas/triangular.hpp> #include <boost/numeric/ublas/banded.hpp> #incl</boost/numeric/ublas/banded.hpp></boost/numeric/ublas/triangular.hpp></boost/numeric/ublas/matrix.hpp>…