boost - 从 Boost::ublas 到 ATLAS，

Question

我使用 Boost::ublas 编写了一个程序，该程序使用了广泛的稀疏矩阵向量乘法。我对它的速度一点也不满意，我想试试 ATLAS。是否有明确的程序来转换代码？为了讨论，我在下面附上了一个非常简单的代码。

#include <boost/numeric/ublas/matrix_sparse.hpp>
#include <boost/numeric/ublas/operation.hpp>

namespace ublas = boost::numeric::ublas;
int main()
{
ublas::compressed_matrix<double> m1(700,700,0);
ublas::compressed_matrix<double> m2(700,700,0);
ublas::compressed_matrix<double> r(700,700,0);
ublas::axpy_prod(m1,m2,r,true);
return 0;
}