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;
}