[petsc-users] sundials and ts

[Jed Brown]

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>