|
|
Ublas : |
From: Thomas Lemaire (thomas.lemaire_at_[hidden])
Date: 2006-09-12 09:24:47
I am successfully using the lapack bindings which offer this functionality,
maybe this could meet your needs...
thomas
On Tuesday 12 September 2006 15:09, Marco wrote:
> Atlas bindings available on
>
> http://news.tiker.net/filebrowser/software/boost-bindings
>
> and described in
>
> http://www.crystalclearsoftware.com/cgi-bin/boost_wiki/wiki.pl?Linear_Algeb
>ra_With_UBLAS
>
> are rising a lot of problems on my debian box.
>
> - they are searching for ublas.h and clapack.h in the wrong place
> - if I correct their includes, I still get
>
> /tmp/ccW4mZZR.o: In function
> `boost::numeric::bindings::atlas::detail::getrf(CBLAS_ORDER, int, int,
> double*, int, int*)':
> a.cpp:(.gnu.linkonce.t._ZN5boost7numeric8bindings5atlas6detail5getrfE11CBLA
>S_ORDERiiPdiPi[boost::numeric::bindings::atlas::detail::getrf(CBLAS_ORDER,
> int, int, double*, int, int*)]+0x30): undefined reference to
> `clapack_dgetrf'
> /tmp/ccW4mZZR.o: In function
> `boost::numeric::bindings::atlas::detail::getri(CBLAS_ORDER, int,
> double*, int, int const*)':
> a.cpp:(.gnu.linkonce.t._ZN5boost7numeric8bindings5atlas6detail5getriE11CBLA
>S_ORDERiPdiPKi[boost::numeric::bindings::atlas::detail::getri(CBLAS_ORDER,
> int, double*, int, int const*)]+0x29): undefined reference to
> `clapack_dgetri
>
> with the example source code
>
> #include <boost/numeric/bindings/atlas/clapack.hpp>
> #include <boost/numeric/bindings/traits/ublas_matrix.hpp>
> #include <boost/numeric/bindings/traits/std_vector.hpp>
>
> namespace ublas = boost::numeric::ublas;
> namespace atlas = boost::numeric::bindings::atlas;
>
> int main() {
> std::size_t n = 5;
> ublas::matrix<double> A (n, n);
> // fill matrix A
>
> std::vector<int> ipiv (n); // pivot vector
> atlas::lu_factor (A, ipiv); // alias for getrf()
> atlas::lu_invert (A, ipiv); // alias for getri()
> }
>
> How could I use ATLAS (I need SVD) with uBLAS?
>
> Thanks,
> Marco
> _______________________________________________
> ublas mailing list
> ublas_at_[hidden]
> http://lists.boost.org/mailman/listinfo.cgi/ublas
-- thomas http://www.laas.fr/~tlemaire