Re: [ublas] ublas::symmetric_marix and potrs for equation solving for positive definite matrices

From: Karl Meerbergen (Karl.Meerbergen_at_[hidden])
Date: 2007-03-01 07:52:07

potrf and potrs require a dense matrix.
symmetric_matrix has a packed format. The same issue was discussed 
earlier this week.

Best,

Karl

Preben Hagh Strunge Holm wrote:

>Hi
>
>I get this error:
>-------------------------------------
>/home/preben/working/RobWork -I /home/preben/working/RobWork/src -c 
>-DNDEBUG test.cpp
>BD Software STL Message Decryptor v2.47a for gcc
>/home/preben/working/RobWork/ext/boostbindings/boost/numeric/bindings/lapack/po
>     sv.hpp: In function
>    `int boost::numeric::bindings::lapack::potrs(
>         const boost::numeric::ublas::symmetric_matrix<
>             double, boost::numeric::ublas::lower
>           , boost::numeric::ublas::column_major
>           , boost::numeric::ublas::unbounded_array<double>
>         > &
>       , boost::numeric::ublas::matrix<
>             double, boost::numeric::ublas::column_major
>           , boost::numeric::ublas::unbounded_array<double>
>         > &
>     )':
>test.cpp:20:   instantiated from here
>/home/preben/working/RobWork/ext/boostbindings/boost/numeric/bindings/lapack/posv.hpp:337:
>     error: incomplete type
>    `boost::STATIC_ASSERTION_FAILURE< false>' used in nested name specifier
>/home/preben/working/RobWork/ext/boostbindings/boost/numeric/bindings/lapack/posv.hpp:337:
>     error: size of array has non-integral type `<type error>'
>/home/preben/working/RobWork/ext/boostbindings/boost/numeric/bindings/traits/ma
>     trix_traits.hpp: In function
>    `int boost::numeric::bindings::traits::leading_dimension(
>         const boost::numeric::ublas::symmetric_matrix<
>             double, boost::numeric::ublas::lower
>           , boost::numeric::ublas::column_major
>           , boost::numeric::ublas::unbounded_array<double>
>         > &
>     )':
>/home/preben/working/RobWork/ext/boostbindings/boost/numeric/bindings/lapack/posv.hpp:294:
>     instantiated from
>    `int boost::numeric::bindings::lapack::detail::potrs(
>         char
>       , const boost::numeric::ublas::symmetric_matrix<
>             double, boost::numeric::ublas::lower
>           , boost::numeric::ublas::column_major
>           , boost::numeric::ublas::unbounded_array<double>
>         > &
>       , boost::numeric::ublas::matrix<
>             double, boost::numeric::ublas::column_major
>           , boost::numeric::ublas::unbounded_array<double>
>         > &
>     )'
>/home/preben/working/RobWork/ext/boostbindings/boost/numeric/bindings/lapack/posv.hpp:347:
>     instantiated from
>    `int boost::numeric::bindings::lapack::potrs(
>         const boost::numeric::ublas::symmetric_matrix<
>             double, boost::numeric::ublas::lower
>           , boost::numeric::ublas::column_major
>           , boost::numeric::ublas::unbounded_array<double>
>         > &
>       , boost::numeric::ublas::matrix<
>             double, boost::numeric::ublas::column_major
>           , boost::numeric::ublas::unbounded_array<double>
>         > &
>     )'
>test.cpp:20:   instantiated from here
>/home/preben/working/RobWork/ext/boostbindings/boost/numeric/bindings/traits/matrix_traits.hpp:131:
>     error: `leading_dimension' is not a member of
>    `boost::numeric::bindings::traits::matrix_traits<
>         const boost::numeric::ublas::symmetric_matrix<
>             double, boost::numeric::ublas::lower
>           , boost::numeric::ublas::column_major
>           , boost::numeric::ublas::unbounded_array<double>
>         >
>     >'
>-------------------------------------
>
>
>When I try to compile this simple example:
>-------------------------------------
>#include <boost/numeric/ublas/matrix.hpp>
>#include <boost/numeric/ublas/symmetric.hpp>
>#include <boost/numeric/ublas/io.hpp>
>
>#include <boost/numeric/bindings/traits/ublas_matrix.hpp>
>#include <boost/numeric/bindings/traits/ublas_symmetric.hpp>
>
>#include <boost/numeric/bindings/lapack/posv.hpp>
>
>#include <vector>
>
>namespace ublas = boost::numeric::ublas;
>namespace lapack = boost::numeric::bindings::lapack;
>
>
>int main() {
>   ublas::symmetric_matrix<double, ublas::lower, ublas::column_major> 
>A(3,3);
>   ublas::matrix<double, ublas::column_major> B(3,1);
>
>   lapack::potrs(A,B);
>
>   return 0;
>}
>-------------------------------------
>
>
>What have I done wrong? What's the missing piece?
>
>
>Thanks
>
>/ Preben
>_______________________________________________
>ublas mailing list
>ublas_at_[hidden]
>http://lists.boost.org/mailman/listinfo.cgi/ublas
>  
>

Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm