Eigen :: MatrixXd typedefの置換

Question

Eigen::MatrixXdとEigen::VectorXdをベクトルと行列で置き換える最も簡単な方法は、long doubleの要素ですか？

コード内のすべての基本浮動小数点変数は、タイプlong doubleです。また、私は行列やベクトルを使うたびに、次のtypedefを使います。

typedef Eigen::VectorXd Vec; 
typedef Eigen::MatrixXd Mat; 

これらのtypedefをどのように切り替えるのが最適ですか？もし私がそれらをそのまま残すとどうなりますか？

Answer 1

Eigenのglobal matrix typedefsに基づいて、独自のtypedefを定義するだけです。

あなたがEigen::MatrixXdを使用して、タイプlong doubleの要素とそれを埋める場合は、これらの値は、精度の損失や、最悪の場合には、オーバーフローエラーになり行列のdouble要素に収まるように狭められます。しかし、多くのアーキテクチャでは、倍精度浮動小数点演算は80ビットの拡張精度で行われるため、結果は同じになる可能性があります。あなたは確かにこれに頼るべきではありません！詳細は、例えばlong double vs doubleをご覧ください。

#include <Eigen/Core> 

typedef Eigen::Matrix< long double, Eigen::Dynamic, 1    > Vec; 
typedef Eigen::Matrix< long double, Eigen::Dynamic, Eigen::Dynamic > Mat; 

int main() 
{ 
    long double ld = 2; 

    Mat m(1,1); 
    m(0,0) = ld; 
}