Re: [ublas] Help with Matrix Inversion

Subject: Re: [ublas] Help with Matrix Inversion
From: Marco Guazzone (marco.guazzone_at_[hidden])
Date: 2010-08-31 10:02:09

On Tue, Aug 31, 2010 at 2:52 AM, Romin Puri <rpuri6_at_[hidden]> wrote:
> 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
>

Hello!

I think you forget to attach files.
If they are very big in size you can use some of the available free
file-hosting sites (e.g., see
http://en.wikipedia.org/wiki/Category:File_sharing)

> 1. Each csv file is a hermitian matrix and is a subset of the next larger
> one.
> 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
>
> Romin
>
> -----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
>>
>> http://savingyoutime.wordpress.com/2009/09/21/c-matrix-inversion-boost
>> ublas/
>>
>> 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.
>

Best,

-- Marco