1c4762a1bSJed Brown 2c4762a1bSJed Brown static char help[] ="Loads a previously saved TS."; 3c4762a1bSJed Brown 4c4762a1bSJed Brown /* 5c4762a1bSJed Brown It loads a TS saved with TSView() 6c4762a1bSJed Brown 7c4762a1bSJed Brown */ 8c4762a1bSJed Brown /* 9c4762a1bSJed Brown Include "petscts.h" to use the PETSc timestepping routines. Note that 10c4762a1bSJed Brown this file automatically includes "petscsys.h" and other lower-level 11c4762a1bSJed Brown PETSc include files. 12c4762a1bSJed Brown */ 13c4762a1bSJed Brown #include <petscts.h> 14c4762a1bSJed Brown 15c4762a1bSJed Brown int main(int argc,char **argv) 16c4762a1bSJed Brown { 17c4762a1bSJed Brown TS ts; /* timestepping context */ 18c4762a1bSJed Brown PetscErrorCode ierr; 19c4762a1bSJed Brown PetscViewer viewer; 20c4762a1bSJed Brown 21c4762a1bSJed Brown ierr = PetscInitialize(&argc,&argv,NULL,help);if (ierr) return ierr; 22*5f80ce2aSJacob Faibussowitsch CHKERRQ(PetscDLLibraryAppend(PETSC_COMM_WORLD,&PetscDLLibrariesLoaded,"advection-diffusion-reaction/ex1")); 23*5f80ce2aSJacob Faibussowitsch CHKERRQ(TSCreate(PETSC_COMM_WORLD,&ts)); 24*5f80ce2aSJacob Faibussowitsch CHKERRQ(PetscViewerBinaryOpen(PETSC_COMM_WORLD,"advection-diffusion-reaction/binaryoutput",FILE_MODE_READ,&viewer)); 25*5f80ce2aSJacob Faibussowitsch CHKERRQ(TSLoad(ts,viewer)); 26*5f80ce2aSJacob Faibussowitsch CHKERRQ(PetscViewerDestroy(&viewer)); 27*5f80ce2aSJacob Faibussowitsch /* CHKERRQ(PetscFPTView(0)); */ 28*5f80ce2aSJacob Faibussowitsch CHKERRQ(TSSetFromOptions(ts)); 29*5f80ce2aSJacob Faibussowitsch CHKERRQ(TSSetUp(ts)); 30*5f80ce2aSJacob Faibussowitsch CHKERRQ(TSView(ts,PETSC_VIEWER_STDOUT_WORLD)); 31*5f80ce2aSJacob Faibussowitsch CHKERRQ(TSSolve(ts,NULL)); 32*5f80ce2aSJacob Faibussowitsch CHKERRQ(TSDestroy(&ts)); 33c4762a1bSJed Brown ierr = PetscFinalize(); 34c4762a1bSJed Brown return ierr; 35c4762a1bSJed Brown } 36