xref: /petsc/include/petscts.h (revision 67d5d79c3450e77433783c7025980752e84ea2a2)
1*67d5d79cSLois Curfman McInnes /* $Id: ts.h,v 1.31 1998/07/23 22:50:59 bsmith Exp curfman $ */
2818ad0c1SBarry Smith /*
3f64a0f93SLois Curfman McInnes    User interface for the timestepping package. This is package
4f64a0f93SLois Curfman McInnes    is for use in solving time-dependent PDEs.
5818ad0c1SBarry Smith */
688d459dfSBarry Smith #if !defined(__TS_H)
788d459dfSBarry Smith #define __TS_H
86d84be18SBarry Smith #include "snes.h"
9818ad0c1SBarry Smith 
10f09e8eb9SSatish Balay typedef struct _p_TS* TS;
116d84be18SBarry Smith #define TS_COOKIE PETSC_COOKIE+18
12818ad0c1SBarry Smith 
1382bf6240SBarry Smith #define TS_EULER           "euler"
1482bf6240SBarry Smith #define TS_BEULER          "beuler"
1582bf6240SBarry Smith #define TS_PSEUDO          "pseudo"
16*67d5d79cSLois Curfman McInnes #define TS_CRANK_NICHOLSON "crank-nicholson"
1782bf6240SBarry Smith #define TS_PVODE           "pvode"
1882bf6240SBarry Smith 
1982bf6240SBarry Smith typedef char *TSType;
2082bf6240SBarry Smith 
2119bcc07fSBarry Smith typedef enum {TS_LINEAR, TS_NONLINEAR} TSProblemType;
22818ad0c1SBarry Smith 
236d84be18SBarry Smith extern int TSCreate(MPI_Comm,TSProblemType,TS*);
246d84be18SBarry Smith extern int TSSetType(TS,TSType);
25d4bb536fSBarry Smith extern int TSGetProblemType(TS,TSProblemType*);
266d84be18SBarry Smith extern int TSDestroy(TS);
27818ad0c1SBarry Smith 
2821c89e3eSBarry Smith extern int TSSetMonitor(TS,int(*)(TS,int,double,Vec,void*),void *);
295cd90555SBarry Smith extern int TSClearMonitor(TS);
3082bf6240SBarry Smith extern int TSGetType(TS,TSType*);
31818ad0c1SBarry Smith 
32d4bb536fSBarry Smith extern int TSSetOptionsPrefix(TS, char *);
33d4bb536fSBarry Smith extern int TSAppendOptionsPrefix(TS, char *);
34d4bb536fSBarry Smith extern int TSGetOptionsPrefix(TS, char **);
356d84be18SBarry Smith extern int TSSetFromOptions(TS);
366d84be18SBarry Smith extern int TSSetUp(TS);
37818ad0c1SBarry Smith 
386d84be18SBarry Smith extern int TSSetSolution(TS,Vec);
396d84be18SBarry Smith extern int TSGetSolution(TS,Vec*);
40818ad0c1SBarry Smith 
4121c89e3eSBarry Smith extern int TSSetDuration(TS,int,double);
426d84be18SBarry Smith extern int TSPrintHelp(TS);
43818ad0c1SBarry Smith 
4421c89e3eSBarry Smith extern int TSDefaultMonitor(TS,int,double,Vec,void*);
4521c89e3eSBarry Smith extern int TSStep(TS,int *,double*);
46818ad0c1SBarry Smith 
47c3e30b67SBarry Smith extern int TSSetInitialTimeStep(TS,double,double);
48c3e30b67SBarry Smith extern int TSGetTimeStep(TS,double*);
49d4bb536fSBarry Smith extern int TSGetTime(TS,double*);
5021c89e3eSBarry Smith extern int TSGetTimeStepNumber(TS,int*);
51c3e30b67SBarry Smith extern int TSSetTimeStep(TS,double);
52818ad0c1SBarry Smith 
5321c89e3eSBarry Smith extern int TSSetRHSFunction(TS,int (*)(TS,double,Vec,Vec,void*),void*);
54c3e30b67SBarry Smith extern int TSSetRHSMatrix(TS,Mat,Mat,int (*)(TS,double,Mat*,Mat*,MatStructure*,void*),void*);
55c3e30b67SBarry Smith extern int TSSetRHSJacobian(TS,Mat,Mat,int(*)(TS,double,Vec,Mat*,Mat*,MatStructure*,void*),void*);
56*67d5d79cSLois Curfman McInnes extern int TSSetRHSBoundaryConditions(TS,int (*)(TS,double,Vec,void*),void*);
57818ad0c1SBarry Smith 
58e34be4c2SBarry Smith extern int TSDefaultComputeJacobianWithColoring(TS,double,Vec,Mat*,Mat*,MatStructure*,void*);
59453b9e1bSBarry Smith extern int TSSetRHSJacobianDefault(TS,MatFDColoring,Mat,Mat);
60453b9e1bSBarry Smith extern int TSDefaultComputeJacobianSlow(TS,double,Vec,Mat*,Mat*,MatStructure*,void*);
61e34be4c2SBarry Smith 
62d4bb536fSBarry Smith extern int TSGetRHSMatrix(TS,Mat*,Mat*,void**);
63d4bb536fSBarry Smith extern int TSGetRHSJacobian(TS,Mat*,Mat*,void**);
64d4bb536fSBarry Smith 
65d8345c25SBarry Smith extern int TSPseudoSetTimeStep(TS,int(*)(TS,double*,void*),void*);
66d8345c25SBarry Smith extern int TSPseudoDefaultTimeStep(TS,double*,void* );
67d8345c25SBarry Smith extern int TSPseudoComputeTimeStep(TS,double *);
68d8345c25SBarry Smith 
69d8345c25SBarry Smith extern int TSPseudoSetVerifyTimeStep(TS,int(*)(TS,Vec,void*,double*,int*),void*);
70d8345c25SBarry Smith extern int TSPseudoDefaultVerifyTimeStep(TS,Vec,void*,double*,int*);
71d8345c25SBarry Smith extern int TSPseudoVerifyTimeStep(TS,Vec,double*,int*);
72bc223fc0SBarry Smith extern int TSPseudoSetTimeStepIncrement(TS,double);
730513a670SBarry Smith extern int TSPseudoIncrementDtFromInitialDt(TS);
7421c89e3eSBarry Smith 
75c3e30b67SBarry Smith extern int TSComputeRHSFunction(TS,double,Vec,Vec);
76818ad0c1SBarry Smith 
77488ecbafSBarry Smith extern FList TSList;
7882bf6240SBarry Smith extern int TSRegisterAll(char*);
79cf256101SBarry Smith extern int TSRegisterDestroy(void);
8084cb2905SBarry Smith extern int TSRegisterAllCalled;
81818ad0c1SBarry Smith 
826df38c32SLois Curfman McInnes extern int TSRegister_Private(char*,char*,char*,int(*)(TS));
836df38c32SLois Curfman McInnes #if defined(USE_DYNAMIC_LIBRARIES)
846df38c32SLois Curfman McInnes #define TSRegister(a,b,c,d) TSRegister_Private(a,b,c,0)
856df38c32SLois Curfman McInnes #else
866df38c32SLois Curfman McInnes #define TSRegister(a,b,c,d) TSRegister_Private(a,b,c,d)
876df38c32SLois Curfman McInnes #endif
886df38c32SLois Curfman McInnes 
896d84be18SBarry Smith extern int TSGetSNES(TS,SNES*);
906d84be18SBarry Smith extern int TSGetSLES(TS,SLES*);
91818ad0c1SBarry Smith 
92c3e30b67SBarry Smith extern int TSView(TS,Viewer);
9321c89e3eSBarry Smith 
9421c89e3eSBarry Smith extern int TSSetApplicationContext(TS,void *);
9521c89e3eSBarry Smith extern int TSGetApplicationContext(TS,void **);
9621c89e3eSBarry Smith 
973914022bSBarry Smith extern int TSLGMonitorCreate(char *,char *,int,int,int,int, DrawLG *);
983914022bSBarry Smith extern int TSLGMonitor(TS,int,double,Vec,void *);
993914022bSBarry Smith extern int TSLGMonitorDestroy(DrawLG);
1003914022bSBarry Smith 
10183e2fdc7SBarry Smith /*
10283e2fdc7SBarry Smith        PETSc interface to PVode
10383e2fdc7SBarry Smith */
104330e16abSSatish Balay #define PVODE_UNMODIFIED_GS PVODE_CLASSICAL_GS
105453b9e1bSBarry Smith typedef enum { PVODE_ADAMS, PVODE_BDF } TSPVodeType;
106330e16abSSatish Balay typedef enum { PVODE_MODIFIED_GS = 0, PVODE_CLASSICAL_GS = 1 } TSPVodeGramSchmidtType;
107453b9e1bSBarry Smith extern int TSPVodeSetType(TS,TSPVodeType);
10883e2fdc7SBarry Smith extern int TSPVodeGetPC(TS,PC*);
10990e6f1a8SBarry Smith extern int TSPVodeSetTolerance(TS,double,double);
110330e16abSSatish Balay extern int TSPVodeGetIterations(TS,int *,int *);
11176be9ce4SBarry Smith extern int TSPVodeSetGramSchmidtType(TS,TSPVodeGramSchmidtType);
11276be9ce4SBarry Smith extern int TSPVodeSetGMRESRestart(TS,int);
113df8cb225SBarry Smith extern int TSPVodeSetLinearTolerance(TS,double);
114e1311b90SBarry Smith extern int TSPVodeSetExactFinalTime(TS, PetscTruth);
11583e2fdc7SBarry Smith 
116818ad0c1SBarry Smith #endif
117818ad0c1SBarry Smith 
11821c89e3eSBarry Smith 
11921c89e3eSBarry Smith 
12021c89e3eSBarry Smith 
12121c89e3eSBarry Smith 
122