[ublas] Sparse blas

From: Frank Astier (fastier_at_[hidden])
Date: 2007-03-22 12:34:09

May I ask, any good library for sparse matrix operations??

Frank
fastier_at_[hidden]

On 平成 19/03/22, at 9:29, Nico Galoppo wrote:

> 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]
>
> _______________________________________________
> ublas mailing list
> ublas_at_[hidden]
> http://lists.boost.org/mailman/listinfo.cgi/ublas