[ublas] Help with Matrix Inversion

Subject: [ublas] Help with Matrix Inversion
From: Romin Puri (rpuri6_at_[hidden])
Date: 2010-08-26 22:32:57

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.

Thanks

Romin