Comment symétriser une matrice clairsemée dans Eigen C++?

Question

J'ai une matrice clairsemée A dans Eigen C++. Maintenant, je veux symétriser à une autre matrice clairsemée Asym:

J'espérais que ce serait aussi simple que:

Eigen::SparseMatrix<FLOATDATA> A; 
... 
Eigen::SparseMatrix<FLOATDATA> Asym = 0.5*(A+A.transpose()); // error here 

Mais pour des raisons évidentes, il donne l'erreur d'échec assert suivant:

error: static assertion failed: THE_STORAGE_ORDER_OF_BOTH_SIDES_MUST_MATCH 

Ma question est comment faire proprement l'opération ci-dessus dans Eigen C++?

Answer 1

La meilleure façon de rendre votre compilation de code est d'évaluer la matrice transposée dans un temporaire de l'ordre de stockage correct:

Eigen::SparseMatrix<FLOATDATA> Asym = 0.5*(A+Eigen::SparseMatrix<FLOATDATA>(A.transpose())));