On Wed, Mar 30, 2011 at 18:43, I-Kamenetskiy, Dmitry S < dmitry.s.kamenetskiy at boeing.com> wrote:
> I have run in to the following issue with the PETSc SNES usage. The logic > of the SNESSolve() routine in my code is that the solution vector x should > include ghost points at the boundary of subdomains in the distributed case. Not true, the solvers take "global" vectors, each degree of freedom appears exactly once on exactly one process. Ghost points are typically needed inside of residual evaluation and Jacobian assembly, but they should not be present in the vector you give to the solver. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110331/cdf0a68c/attachment.htm>
