Path: news.mathworks.com!not-for-mail
From: <HIDDEN>
Newsgroups: comp.soft-sys.matlab
Subject: condition number of FE matrix -> improving iterative solver
Date: Sun, 8 Nov 2009 15:50:02 +0000 (UTC)
Organization: The MathWorks, Inc.
Lines: 15
Message-ID: <hd6pba$84c$1@fred.mathworks.com>
Reply-To: <HIDDEN>
NNTP-Posting-Host: webapp-05-blr.mathworks.com
Content-Type: text/plain; charset="ISO-8859-1"
Content-Transfer-Encoding: 8bit
X-Trace: fred.mathworks.com 1257695402 8332 172.30.248.35 (8 Nov 2009 15:50:02 GMT)
X-Complaints-To: news@mathworks.com
NNTP-Posting-Date: Sun, 8 Nov 2009 15:50:02 +0000 (UTC)
X-Newsreader: MATLAB Central Newsreader 1337510
Xref: news.mathworks.com comp.soft-sys.matlab:583350


Hi,

I'm working on a FE software and right now I'm stuck with the solver.

I'm both trying cholmod2 (from Tim Davis' package) and standard pcg.
Right now, i've got a text matrix, which is just a beam with about 30 000dof and cholmod2 wins hands down. However, when using larger matrices, I can't use cholmod2 because there isn't sufficient memory.

However, I noted a few things which actually surpise me.
First of all, when using condest on my FE matrix of the beam, the condition number is in the order of 1e8. When checking the statistics of cholmod2, it returns the reciprocal of the condition number as about 1e-3. What is the cause for the big difference in the condition number? I have been told that the condition number for structural problem usually is somewhere arounc 1e4 to 1e5, so I think that the 1e8 is rather high. Especially since it's a simple problem with a regular mesh.

For the pcg I'm using a incomplete lu factorization with zero fill-in. Why is pcg so much slower than the direct solver, even if I use a preconditioner?

Are there any other tricks to improve the condition number of the my FE matrix, or ways to improve the solution speed of pcg? I need it for problems with more dof? Or is there a way to make cholmod2 use less memory (though it already applies reordering...)

Thanks for any help on my problem!!