numerical solution of the 3D Poisson equation using the incomplete Cholesky factorization conjugate gradient method