[petsc-users] sundials and ts

[Konstantinos Kontzialis]

Dear all,

I want to use sundials with -ts_sundials_type adams, but I get the 
following error:

Timestep   0: dt = 0.01, T = 0, Res[rho] = 7.80568e-18, Res[rhou] = 
0.0258457, Res[rhov] = 1.25e-12, Res[E] = 4.12368e-09, CFL = 1999.99
[0]PETSC ERROR: --------------------- Error Message 
------------------------------------
[0]PETSC ERROR: Invalid argument!
[0]PETSC ERROR: Vector is not ghosted!
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.2.0, Patch 5, Sat Oct 29 
13:45:54 CDT 2011
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR: 
------------------------------------------------------------------------
[0]PETSC ERROR: ./hoac on a linux-gnu named PlusSodaL by kontzialis Thu 
Nov 17 13:32:48 2011
[0]PETSC ERROR: Libraries linked from 
/home/kontzialis/petsc-3.2-p5/linux-gnu-c-debug/lib
[0]PETSC ERROR: Configure run at Sat Nov  5 20:58:12 2011
[0]PETSC ERROR: Configure options --with-debugging=1 
---with-mpi-dir=/usr/lib64/mpich2/bin --with-shared-libraries 
--with-shared-libraries --with-large-file-io=1 --with-precision=double 
--with-blacs=1 --download-blacs=yes --download-f-blas-lapack=yes 
--with-plapack=1 --download-plapack=yes --with-scalapack=1 
--download-scalapack=yes --with-superlu=1 --download-superlu=yes 
--with-superlu_dist=1 --download-superlu_dist=yes --with-ml=1 
--download-ml=yes --with-umfpack=1 --download-umfpack=yes --with-mpi=1 
--download-mpich=1 --with-sundials=1 --download-sundials=1 
--with-parmetis=1 --download-parmetis=1 --with-hypre=1 --download-hypre=1
[0]PETSC ERROR: 
------------------------------------------------------------------------
[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
application called MPI_Abort(comm=0x84000000, 62) - process 0
[cli_0]: aborting job:
application called MPI_Abort(comm=0x84000000, 62) - process 0
==15195==
==15195== HEAP SUMMARY:
==15195==     in use at exit: 3,110,164 bytes in 18,559 blocks
==15195==   total heap usage: 36,568 allocs, 18,009 frees, 300,373,551 
bytes allocated
==15195==
==15195== LEAK SUMMARY:
==15195==    definitely lost: 39,298 bytes in 24 blocks
==15195==    indirectly lost: 24 bytes in 3 blocks
==15195==      possibly lost: 0 bytes in 0 blocks
==15195==    still reachable: 3,070,842 bytes in 18,532 blocks
==15195==         suppressed: 0 bytes in 0 blocks
==15195== Rerun with --leak-check=full to see details of leaked memory
==15195==
==15195== For counts of detected and suppressed errors, rerun with: -v
==15195== ERROR SUMMARY: 0 errors from 0 contexts (suppressed: 6 from 6)

Any suggestions?

Thank you,

Kostas