Boost logo

Ublas :

Subject: Re: [ublas] Help with Matrix Inversion
From: Romin Puri (rpuri6_at_[hidden])
Date: 2010-08-30 20:52:12

My apologies to all subscribers for sending such a large message.

-----Original Message-----
From: Romin Puri [mailto:rpuri6_at_[hidden]]
Sent: Monday, August 30, 2010 5:37 PM
To: 'ublas mailing list'
Subject: RE: [ublas] Help with Matrix Inversion

Hi Marcos,

I think I have a serious memory corruption issue. Unfortunately I do not
have access to Purify or like otherwise I would have tracked it down.

I am attaching test code for your preview along with data files

1. Each csv file is a hermitian matrix and is a subset of the next larger
2. I have included a 963x963 data file as the problem becomes more
pronounced after this inversion.

Please let me know if you know any additional information. Your help would
be greatly appreciated. I am on boost release 1.43 and am still using Visual
Studio 2005 on Windows 7 32-bit ( We are a very poor startup).

While it is very possible I am making a subtle mistake with my use of boost,
your inputs would be greatly appreciated.

Thank you very much for your help

Best regards


-----Original Message-----
From: ublas-bounces_at_[hidden] [mailto:ublas-bounces_at_[hidden]]
On Behalf Of Marco Guazzone
Sent: Friday, August 27, 2010 1:29 AM
To: ublas mailing list
Subject: Re: [ublas] Help with Matrix Inversion

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
> ublas/
> and trying out Hermition matricies of different sizes to get
> timings/memory requirements. The requirement is to solve a NxN where N =
> 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.


-- Marco
ublas mailing list
Sent to: rpuri6_at_[hidden]