1c4762a1bSJed Brown static char help[] = "Loads a previously saved TS.";
2c4762a1bSJed Brown
3c4762a1bSJed Brown /*
4c4762a1bSJed Brown It loads a TS saved with TSView()
5c4762a1bSJed Brown
6c4762a1bSJed Brown */
7c4762a1bSJed Brown /*
8c4762a1bSJed Brown Include "petscts.h" to use the PETSc timestepping routines. Note that
9c4762a1bSJed Brown this file automatically includes "petscsys.h" and other lower-level
10c4762a1bSJed Brown PETSc include files.
11c4762a1bSJed Brown */
12c4762a1bSJed Brown #include <petscts.h>
13c4762a1bSJed Brown
main(int argc,char ** argv)14*d71ae5a4SJacob Faibussowitsch int main(int argc, char **argv)
15*d71ae5a4SJacob Faibussowitsch {
16c4762a1bSJed Brown TS ts; /* timestepping context */
17c4762a1bSJed Brown PetscViewer viewer;
18c4762a1bSJed Brown
19327415f7SBarry Smith PetscFunctionBeginUser;
209566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc, &argv, NULL, help));
219566063dSJacob Faibussowitsch PetscCall(PetscDLLibraryAppend(PETSC_COMM_WORLD, &PetscDLLibrariesLoaded, "advection-diffusion-reaction/ex1"));
229566063dSJacob Faibussowitsch PetscCall(TSCreate(PETSC_COMM_WORLD, &ts));
239566063dSJacob Faibussowitsch PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD, "advection-diffusion-reaction/binaryoutput", FILE_MODE_READ, &viewer));
249566063dSJacob Faibussowitsch PetscCall(TSLoad(ts, viewer));
259566063dSJacob Faibussowitsch PetscCall(PetscViewerDestroy(&viewer));
269566063dSJacob Faibussowitsch /* PetscCall(PetscFPTView(0)); */
279566063dSJacob Faibussowitsch PetscCall(TSSetFromOptions(ts));
289566063dSJacob Faibussowitsch PetscCall(TSSetUp(ts));
299566063dSJacob Faibussowitsch PetscCall(TSView(ts, PETSC_VIEWER_STDOUT_WORLD));
309566063dSJacob Faibussowitsch PetscCall(TSSolve(ts, NULL));
319566063dSJacob Faibussowitsch PetscCall(TSDestroy(&ts));
329566063dSJacob Faibussowitsch PetscCall(PetscFinalize());
33b122ec5aSJacob Faibussowitsch return 0;
34c4762a1bSJed Brown }
35