Help with Biconjugate Gradient to solve a system

 Register Blogs Members List Search Today's Posts Mark Forums Read

 December 5, 2003, 06:38 Help with Biconjugate Gradient to solve a system #1 Isa Guest   Posts: n/a Hello to everybody: I have the following doubt: I want to solve the following system: Ax=f I use an iterative solver, at first A is non-symmetric and positive. I use the Biconjugate gradient with A as preconditioner matrix. Several iterations later, A gets to be negative in some of its elements and it is then when the program aborts. If I used a preconditioner matrix that was always symmetric and positive (for example a diagonal one like a lumped-mass matrix), would this problem disapear? What is the right solution? Thank-you in advance. Best regards Isa

 December 5, 2003, 14:51 Re: Help with Biconjugate Gradient to solve a syst #2 shobhit agarwal Guest   Posts: n/a hi! for using a preconditioned conjugate gradient matrix there are three requirements. (1). The main(center)diagonal elements should never be zero. (2).The absolute value of the main diagonal element should be greater than all the values of the elements in that row.i.e. example abs(a22)>a2j (3).The signs of all the diagonal elements should be same and also matrix should be symmetric in terms of signs in your case since the sign of some elements of main diagonal change after a few iterations then the rule-3 gets violated hence the problem begins....so it seems that there is some formulation problem check that,because even if you use the preconditioned gradient method the requiremnt no.3 is there for more refference see MSIMSL library and particularly help files jcgrc. shobhit

 December 20, 2003, 08:19 Re: Help with Biconjugate Gradient to solve a syst #3 kenn Guest   Posts: n/a you should use GMRES or GPBiCG(m,l) linears solvers. these are two best general-purposed linear solvers.

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post jussi OpenFOAM Installation 0 April 24, 2008 14:25 Mahesh Main CFD Forum 0 March 26, 2005 16:32 Paul Hsieh Main CFD Forum 3 November 11, 2003 05:52 west_wing FLUENT 0 August 25, 2003 10:00 Atit Koonsrisuk Main CFD Forum 2 January 10, 2002 11:52

All times are GMT -4. The time now is 17:56.