|
|
Ublas : |
Subject: Re: [ublas] boost::ublas lu_substitute runtime error
From: Ungermann, Jörn (j.ungermann_at_[hidden])
Date: 2011-08-16 09:45:49
Hello!
The LU substitution works fine for our testcases. What may not be obvious is that the routine you link to requires the inverse matrix to be of the correct size.
You might want to add some asserts to this effect.
If this doesn't help, I'd be willing to debug some simple testcase that causes this problem on your end, if you provide it.
Regards,
Joern
> -----Original Message-----
> From: ublas-bounces_at_[hidden] [mailto:ublas-
> bounces_at_[hidden]] On Behalf Of mga-mga
> Sent: Dienstag, 16. August 2011 14:34
> To: ublas_at_[hidden]
> Subject: [ublas] boost::ublas lu_substitute runtime error
>
> 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<const_matrix_type, upper> (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.324
> 26e+08,2.29091e+09,4.01989e+10,7.13142e+11,1.27611e+13),(472042,2.29091
> e+09,4.01989e+10,7.13142e+11,1.27611e+13,2.29941e+14),(7.8021e+06,4.019
> 89e+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.
> _______________________________________________
> ublas mailing list
> ublas_at_[hidden]
> http://lists.boost.org/mailman/listinfo.cgi/ublas
> Sent to: j.ungermann_at_[hidden]
------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------
Forschungszentrum Juelich GmbH
52425 Juelich
Sitz der Gesellschaft: Juelich
Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498
Vorsitzender des Aufsichtsrats: MinDirig Dr. Karl Eugen Huthmacher
Geschaeftsfuehrung: Prof. Dr. Achim Bachem (Vorsitzender),
Dr. Ulrich Krafft (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt,
Prof. Dr. Sebastian M. Schmidt
------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------