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 15*d71ae5a4SJacob Faibussowitsch int main(int argc, char **argv) 16*d71ae5a4SJacob Faibussowitsch { 17c4762a1bSJed Brown TS ts; /* timestepping context */ 18c4762a1bSJed Brown PetscViewer viewer; 19c4762a1bSJed Brown 20327415f7SBarry Smith PetscFunctionBeginUser; 219566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc, &argv, NULL, help)); 229566063dSJacob Faibussowitsch PetscCall(PetscDLLibraryAppend(PETSC_COMM_WORLD, &PetscDLLibrariesLoaded, "advection-diffusion-reaction/ex1")); 239566063dSJacob Faibussowitsch PetscCall(TSCreate(PETSC_COMM_WORLD, &ts)); 249566063dSJacob Faibussowitsch PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD, "advection-diffusion-reaction/binaryoutput", FILE_MODE_READ, &viewer)); 259566063dSJacob Faibussowitsch PetscCall(TSLoad(ts, viewer)); 269566063dSJacob Faibussowitsch PetscCall(PetscViewerDestroy(&viewer)); 279566063dSJacob Faibussowitsch /* PetscCall(PetscFPTView(0)); */ 289566063dSJacob Faibussowitsch PetscCall(TSSetFromOptions(ts)); 299566063dSJacob Faibussowitsch PetscCall(TSSetUp(ts)); 309566063dSJacob Faibussowitsch PetscCall(TSView(ts, PETSC_VIEWER_STDOUT_WORLD)); 319566063dSJacob Faibussowitsch PetscCall(TSSolve(ts, NULL)); 329566063dSJacob Faibussowitsch PetscCall(TSDestroy(&ts)); 339566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 34b122ec5aSJacob Faibussowitsch return 0; 35c4762a1bSJed Brown } 36