Article View: sci.math.num-analysis
Article #6Pre-Conditioned CG NOT CONVERGING
From: Dennis Moon
Date: Tue, 27 Oct 1992 14:03
Date: Tue, 27 Oct 1992 14:03
21 lines
1224 bytes
1224 bytes
A friend without news access is having trouble with solving the matrices in his magneto-hydrodynamic code. The equation is the standard Ax=b, where the solution vector x is 4096 elements. The matrix is sparse, with 5 non-zero elements per row. He is using a CG solver with LU factorization. For most time steps it converges in about 10 iterations, but on some timesteps it will not converge no matter how many iterations. He has tried both a hand coded solver and routines from the IBM Risc6000 ESSL library with the same result. What other solution techniques should he be looking at? Is this an indication that in some sense the problem is not well posed? What matrix characteristics can be used to pre-diagnose this problem? Any comments would be welcome, his thesis depends on this Thanks Dennis Moon dam@ssesco.com
Technical Details
Message-ID:
<Bws9qM.F7o@ssesco.com>Path:
arc1-ca.pugleaf.net!sparky!uunet!sun-barr!cs.utexas.edu!qt.cs.utexas.edu!yale.edu!spool.mu.edu!umn.edu!uum1!salsa.ssesco.com!ssesco.com!damShow Full Headers
Complete Article Headers
Newsgroups: sci.math.num-analysis Path: arc1-ca.pugleaf.net!sparky!uunet!sun-barr!cs.utexas.edu!qt.cs.utexas.edu!yale.edu!spool.mu.edu!umn.edu!uum1!salsa.ssesco.com!ssesco.com!dam From: dam@ssesco.com (Dennis Moon) Subject: Pre-Conditioned CG NOT CONVERGING Message-ID: <Bws9qM.F7o@ssesco.com> Sender: news@ssesco.com Organization: SSESCO Date: Tue, 27 Oct 1992 14:03:57 GMT Lines: 21