|[eigen] Large matrix problem|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
- To: <eigen@xxxxxxxxxxxxxxxxxxx>
- Subject: [eigen] Large matrix problem
- From: "Jian Yang" <Jian.Yang@xxxxxxxxxxx>
- Date: Thu, 1 Jul 2010 16:19:05 +1000
- Thread-index: AcsVWhnoxl3rJH5QTLu6JAqRqg4cVwDipiug
- Thread-topic: Large matrix problem
I'm trying to create a 14000 x 14000*22 matrix using the following code
MatrixXf a=MatrixXf::Zero(14000, 14000*22);
The program stopped with the error: "Segmentation fault".
I have 64G RAM in my machine. I think I have enough memory to create
Anyone knows what's going on?
From: Listengine [mailto:listengine@xxxxxxxxxxxxxxxxx] On Behalf Of
Sent: Sunday, 27 June 2010 4:05 AM
Subject: Re: [eigen] [patch] LDLt decomposition with rank-deficient
Thanks for the patches and sorry that I don't have time to look at
them today: we are currently in an incredible rush to get beta1 out
the door, and that means a lot of documentation work. I think that
your stuff is OK to apply post-beta1 anyway. Just wanted to let you
know why we might have long response times these days.
2010/6/26 Ben Goodrich <bgokgm@xxxxxxxxxxxxxx>:
> On Fri, Jun 25, 2010 at 1:50 AM, Ben Goodrich <bgokgm@xxxxxxxxxxxxxx>
>>>> 5) I copied-and-pasted a block inside /test/cholesky.cpp and
>>>> the pivot=false option. It seems to work when you do ./check.sh
>>>> cholesky. I did some other tests locally with singular matrices,
>>>> /test/cholesky.cpp does not seem to have any tests with singular
>>>> matrices, so maybe some should be added?
>>> why not.
>> I have not added the singular tests yet, but I can do that soon.
> This patch does so. It seems to work okay with both the Pivoting and
> NoPivoting options from the previous patch.
>>> You should also make the solve function skips the transpositions
>>> no pivoting has been computed.
>> I have not done this yet either. What did you decide about making a
>> new class versus putting a flag in the class definition?