Help with Biconjugate Gradient to solve a system
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 nonsymmetric 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 lumpedmass matrix), would this problem disapear? What is the right solution? Thankyou in advance. Best regards Isa 
Re: Help with Biconjugate Gradient to solve a syst
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 rule3 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

Re: Help with Biconjugate Gradient to solve a syst
you should use GMRES or GPBiCG(m,l) linears solvers. these are two best generalpurposed linear solvers.

All times are GMT 4. The time now is 09:02. 