Subject: Re: [ublas] Help with Matrix Inversion
From: Marco Guazzone (marco.guazzone_at_[hidden])
Date: 2010-08-27 04:28:47
On Fri, Aug 27, 2010 at 4:32 AM, Romin Puri <rpuri6_at_[hidden]> wrote:
> I am using the code as outlined in this link
> and trying out Hermition matricies of different sizes to get timings/memory
> requirements. The requirement is to solve a NxN where N = 1000.
> While Matrix sizes upto N =30 work just fine, at N=50 I get this error. The
> data set is the same (truncated to ensure that Hermitian Matrix is of the
> desired size).
> Check failed in file c:\program
> files\boost\boost_1_43_0\boost\numeric\ublas\lu.hpp at line 294:
> detail::expression_type_check (prod (triangular_adaptor<const_matrix_type,
> upper> (m), e), cm2)
> internal logic "
> Any pointers would be most helpful.
Don't know the cause of your problem, but if you submit a full working
c++ test file I can give you a check.