Boost logo

Ublas :

From: Song Chang (chang8217_at_[hidden])
Date: 2008-01-16 10:15:54

I am using LU decomposition to solve linear system with ublas,
and I just cannot figure the error apeared like this:

Check failed in file d:\lib\boost_1_34_1\boost\numeric\ublas\lu.hpp at
line 272:
detail::expression_type_check (prod
(triangular_adaptor<const_matrix_type, upper
> (m), e), cv2)

lu_factorize(A,pm) seems working well, but lu_substitute(A, pm, x)
sometimes throws
the exception above.

my code:
template<class T>
int lusolver(ublas::matrix<T>& A, ublas::vector<T>& x, ublas::vector<T>& b){
x = b;
ublas::permutation_matrix<std::size_t> pm(A.size1());
int res = ublas::lu_factorize(A,pm);
if( res != 0 ) return 1;
ublas::lu_substitute(A, pm, x);
return 0;}

anybody can help me? thanks a lot.

Easy + Joy + Powerful = Yahoo! Bookmarks x Toolbar