c++ - uBlasは非常に遅い

Question

重複の可能性:
BOOST uBLAS マトリックス製品が非常に遅い

uBlas ライブラリを (boost で) 使用しようとしています。

複数の 2 つの大きな行列に簡単なコードを書きました。結果が遅すぎる、MatlabとOctaveでパフォーマンスを比較したところ、極端に遅かった。問題を解決するのを手伝ってくれませんか

#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/io.hpp>
#include <boost/numeric/ublas/operation.hpp>
int main () {
   using namespace boost::numeric::ublas;
   const int D = 1000;
   matrix<double> a (D, D);
   matrix<double> b (D, D);
   matrix<double> c (D, D);
   for (unsigned i = 0; i < a.size1 (); ++ i)
      for (unsigned j = 0; j < a.size2 (); ++ j){
         a (i, j) = ((i-j)<150)?1:0;   
         b (i, j) = 3 * i + j;
      }
    std::cout <<"Hi\n";
    axpy_prod(a, b,c, false);  // C = A * B
    //noalias(c) = prod(a,b);
    return 0;
}