Poisson eq and conjugate gradient method
Hi,
I got some problems trying to simulate a Poisson equation in 2D with periodic boundary condition (for testing purpose). I use finite element discretization with bilinear functions on a regular mesh and I solve the equation with a conjugate gradient iterative solver. The solver doesn't seem to converge quickly. Is it expected and well known? I am also looking for a preconditioner. Is it possible to find one which can be use on a element level without building the global matrix (and which is efficient of course)? Thanks for help Gary 
Re: Poisson eq and conjugate gradient method
Oups I forgot to say that I impose a nodal value inside the domain.

Re: Poisson eq and conjugate gradient method
If the matrix (Ax=b) you finally have to solve is symmetric and positive definite, it will converge at any intial values of x0 (initially assumed approximate solution vector). And for a preconditioner, I think, Jacobi or SSOR can be used. And ICCG, ILUCG, MILUCG have been reported to give more fast convergence. But the convergence rates depend on the matrix type you have to solve. So, at first, check the matrix you have to solve.

All times are GMT 4. The time now is 00:25. 