<div class="gmail_quote">On Thu, Nov 17, 2011 at 05:33, Konstantinos Kontzialis <span dir="ltr"><<a href="mailto:ckontzialis@lycos.com">ckontzialis@lycos.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":1bp">[0]PETSC ERROR: VecGhostUpdateBegin() line 170 in /home/kontzialis/petsc-3.2-p5/<u></u>src/vec/vec/impls/mpi/<u></u>commonmpvec.c<br>
[0]PETSC ERROR: residual() line 43 in "unknowndirectory/"../src/<u></u>residual.c<br>
[0]PETSC ERROR: base_residual_implicit() line 28 in "unknowndirectory/"../src/<u></u>base_residual_implicit.c<br>
[0]PETSC ERROR: TSComputeIFunction() line 339 in /home/kontzialis/petsc-3.2-p5/<u></u>src/ts/interface/ts.c<br>
[0]PETSC ERROR: TSFunction_Sundials() line 99 in /home/kontzialis/petsc-3.2-p5/<u></u>src/ts/impls/implicit/<u></u>sundials/sundials.c</div></blockquote></div><br><div>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.</div>
<div><br></div><div>Also, if your IFunction does not have the form</div><div><br></div><div>G(x,xdot) = xdot + F(x)</div><div><br></div><div>(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.</div>