On Thu, Nov 17, 2011 at 05:33, Konstantinos Kontzialis < ckontzialis at lycos.com> wrote:
> [0]PETSC ERROR: VecGhostUpdateBegin() line 170 in > /home/kontzialis/petsc-3.2-p5/**src/vec/vec/impls/mpi/**commonmpvec.c > [0]PETSC ERROR: residual() line 43 in "unknowndirectory/"../src/** > residual.c > [0]PETSC ERROR: base_residual_implicit() line 28 in > "unknowndirectory/"../src/**base_residual_implicit.c > [0]PETSC ERROR: TSComputeIFunction() line 339 in > /home/kontzialis/petsc-3.2-p5/**src/ts/interface/ts.c > [0]PETSC ERROR: TSFunction_Sundials() line 99 in > /home/kontzialis/petsc-3.2-p5/**src/ts/impls/implicit/** > sundials/sundials.c > We can't pass ghosted Vecs through the Sundials interface. You'll have to copy into a ghosted work vector if you use VecGhost in your function. Also, if your IFunction does not have the form G(x,xdot) = xdot + F(x) (that is, if dG/dxdot is not the identitiy), then this method will not produce correct answers. Unfortunately, this is a limitation of the method/implementation. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111117/908514b3/attachment-0001.htm>
