Boostez le produit matriciel/vecteur uBLAS

Question

Quelqu'un peut-il donner un exemple de la façon d'utiliser le produit uBLAS pour multiplier les choses? Ou s'il y a une bibliothèque de matrice C++ plus agréable, vous pouvez me recommander aussi. Cela se transforme en un mal de tête majeur.

Voici mon code:

vector<double> myVec(scalar_vector<double>(3)); 
matrix<double> myMat(scalar_matrix<double>(3,3,1)); 
matrix<double> temp = prod(myVec, myMat); 

est ici l'erreur:

cannot convert from 'boost::numeric::ublas::matrix_vector_binary1<E1,E2,F>' to 'boost::numeric::ublas::matrix<T>' 

j'ai épuisé ma recherche. Stackoverflow a une question à propos de cette here. La documentation Boost a un exemple here. J'ai copié le code à partir d'un exemple, mais cela ne me sert à rien car le modèle de magie qui fonctionne pour stdout ne me sert à rien.

#include <boost/numeric/ublas/matrix.hpp> 
#include <boost/numeric/ublas/io.hpp> 

int main() { 
    using namespace boost::numeric::ublas; 
    matrix<double> m (3, 3); 
    vector<double> v (3); 
    for (unsigned i = 0; i < std::min (m.size1(), v.size()); ++ i) { 
     for (unsigned j = 0; j < m.size2(); ++ j) 
      m (i, j) = 3 * i + j; 
     v (i) = i; 
    } 

    std::cout << prod (m, v) << std::endl; 
    std::cout << prod (v, m) << std::endl; 
} 

Answer 1

Je ne l'ai pas regardé beaucoup à vous uBLAS Boost, mais Eigen est agréable, et a une bonne performance.

Answer 2

Le produit d'un vecteur et d'une matrice est un vecteur et non une matrice.