可能重复:
BOOST uBLAS 矩阵乘积极慢
我正在尝试使用 uBlas 库(在 boost 中)。
我为多个两个大型矩阵编写了一个简单的代码。结果太慢了,当我将性能与 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;
}