|
|
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
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.
Cheers,
-- Marco
_______________________________________________
ublas mailing list
ublas_at_[hidden]
http://lists.boost.org/mailman/listinfo.cgi/ublas
Sent to: rpuri6_at_[hidden]