なんだかublasではまった日誌と化してきたが. でも大したことでなくてもここに記すことでだれかの手間は省けるかもしれない.Assertion failed in file /usr/include/boost/numeric/ublas/matrix_sparse.hpp at line 2679: !preservecompressed_matrixをresiz…

uBLASを使って複素数の線形方程式を解きたいとき, 通常のmatrixなら問題ないのにcompressed_matrixでlu_factorizeしようとするとエラーが出てコンパイルできない./usr/include/boost/numeric/ublas/lu.hpp:146: error: no match for 'operator/' in ...原因…

ublasでLU分解を行う際に使用するpermutation_matrixについ てはまったのでメモしておく. permutation_matrixはlu_factorize, lu_substituteと同様にlu.hppに含まれる. 問題はlu_factorize(m, pm)を行った際に, 行列mがsingularでないのにアサーション(BOOST…