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
