От Boost :: ublas до ATLAS, - PullRequest
Новая
       46

От Boost :: ublas до ATLAS,

1 голос
/ 25 августа 2011

Я написал программу, использующую 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;
}
...