[ublas] boost::ublas lu_substitute runtime error

Subject: [ublas] boost::ublas lu_substitute runtime error
From: mga-mga (mga-mga_at_[hidden])
Date: 2011-08-16 08:34:07

Hello,

I'm trying to use LU factorization to invert matrices of <double> items, as
describer here:
http://www.crystalclearsoftware.com/cgi-bin/boost_wiki/wiki.pl?LU_Matrix_Inversion.
That seems to be working OK with some matrices, but sometimes I'v got the
following runtime error on lu_substitute:

Check failed in file /usr/local/include/boost/numeric/ublas/lu.hpp at line
294:
detail::expression_type_check (prod
(triangular_adaptor&lt;const_matrix_type, upper&gt; (m), e), cm2)

I.e. I'v got such errors for the following 6x6 matrix:

Matrix:
[6,6]((15,29700,472042,7.8021e+06,1.32426e+08,2.29091e+09),(29700,1.32426e+08,2.29091e+09,4.01989e+10,7.13142e+11,1.27611e+13),(472042,2.29091e+09,4.01989e+10,7.13142e+11,1.27611e+13,2.29941e+14),(7.8021e+06,4.01989e+10,7.13142e+11,1.27611e+13,2.29941e+14,4.16694e+15),(1.32426e+08,7.13142e+11,1.27611e+13,2.29941e+14,4.16694e+15,7.58705e+16),(2.29091e+09,1.27611e+13,2.29941e+14,4.16694e+15,7.58705e+16,1.38694e+18))

What could be possibly wrong with types here? What may cause such kind of
errors?

Thank you in advance!

--
View this message in context: http://boost.2283326.n4.nabble.com/boost-ublas-lu-substitute-runtime-error-tp3747071p3747071.html
Sent from the Boost - uBLAS mailing list archive at Nabble.com.