[ublas] internal_logic() exception

From: Manoj Rajagopalan (rmanoj_at_[hidden])
Date: 2006-04-03 23:55:15

hi

   I am trying to solve a dense system using the LU decomposition 
procedure. Why am I getting the following compile time error?

Assertion failed in file /usr/include/boost/numeric/ublas/lu.hpp at line 
271:
detail::expression_type_check (prod 
(triangular_adaptor<const_matrix_type, upper> (m), e), cv2)
terminate called after throwing an instance of 
'boost::numeric::ublas::internal_logic'
   what():  internal logic

My code is as follows:

matrix<double> A;
vector<double> x;
// initialize both
permutation_matrix<> piv(A.size1());
lu_factorize(A, piv); // with partial pivoting, lu.hpp
lu_substitute(A, piv, x); // solve the system

thanks in advance,
Manoj