Substitutionen für Eigen :: MatrixXd typedefs

Question

Was ist der einfachste Weg, alle Eigen::MatrixXd s und Eigen::VectorXd s mit Vektoren und Matrizen zu ersetzen, die long double Elemente haben?

Jede grundlegende Gleitkommavariable in meinem Code ist vom Typ long double. Außerdem verwende ich jedes Mal, wenn ich eine Matrix oder einen Vektor verwende, die folgenden typedefs.

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

Was ist das Beste, um diese Typdefs zu wechseln? Was passiert, wenn ich sie so belasse wie sie sind?

Answer 1

Definieren Sie einfach Ihre eigenen Typedefs basierend auf Eigen global matrix typedefs.

Wenn Sie Eigen::MatrixXd verwenden und füllen Sie es mit Elementen des Typs long double, werden diese Werte in die verengte double Elemente der Matrix passen, was zu einem Verlust an Präzision führt oder, im schlimmsten Fall, Überlauffehler. Bei vielen Architekturen wird Gleitkomma-Arithmetik mit doppelter Genauigkeit in 80-Bit-Extended-Genauigkeit ausgeführt, so dass die Ergebnisse identisch sein können. Darauf sollten Sie sich sicher nicht verlassen! Für mehr siehe z. B. 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; 
}