[ublas] More inputs on internal_logic() exception

From: Manoj Rajagopalan (rmanoj_at_[hidden])
Date: 2006-04-04 00:03:54

I must add: This compile-time exception occurs when i switch to larger 
sizes of A and x. For a small size I get no error. What am I missing?

thanks in advance,
Manoj

-------- Original Message --------
Subject: [ublas] internal_logic() exception
Date: Mon, 03 Apr 2006 23:55:15 -0400
From: Manoj Rajagopalan <rmanoj_at_[hidden]>
Reply-To: ublas mailing list <ublas_at_[hidden]>
To: ublas_at_[hidden]

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