Re: [ublas] prod() performance question

From: Nico Galoppo (nico_at_[hidden])
Date: 2007-03-22 12:29:10

Here:

http://tinyurl.com/37h28v

#include <boost/numeric/ublas/matrix.hpp>
namespace ublas = boost::numeric::ublas;

#include <boost/numeric/bindings/atlas/cblas3.hpp>
#include <boost/numeric/bindings/traits/ublas_matrix.hpp>
namespace atlas = boost::numeric::bindings::atlas;

int main()
{
     size_t n = 100;
     ublas::matrix<double, ublas::column_major> A(n,n),B(n,n),C(n,n);

     // fill A and B

     atlas::gemm(A,B,C);
}

Stephen Gross wrote:
>> BTW, it really pays off to use the ATLAS bindings for this 
>> kind of dense 
>> product. It is orders of magnitude faster than the default 
>> ublas implementation.
> 
> Can you show me a simple example of how that would work?
> 
> _______________________________________________
> ublas mailing list
> ublas_at_[hidden]
> http://lists.boost.org/mailman/listinfo.cgi/ublas

-- 
Nico Galoppo        UNC-CH PhD. student        http://www.ngaloppo.org
                       [phone]