Re: [ublas] Help with Matrix Inversion

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
>
> http://savingyoutime.wordpress.com/2009/09/21/c-matrix-inversion-boostublas/
>
> 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.
>

Hi,

Don't know the cause of your problem, but if you submit a full working
c++ test file I can give you a check.

Cheers,

-- Marco