从Boost:：ublas到ATLAS，

我已经使用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;
}