/* $Id: ts.h,v 1.1 1996/01/01 19:40:02 bsmith Exp bsmith $ */
/*
    User interface for the time-stepping package. This is package
  is for use in solving time dependent PDES.
*/
#if !defined(__TS_PACKAGE)
#define __TS_PACKAGE
#include "snes.h"

typedef struct _TS* TS;
#define TS_COOKIE PETSC_COOKIE+18

typedef enum { TS_EULER, TS_BEULER, TS_PSEUDO} TSType;
typedef enum { TS_LINEAR_CONSTANT_MATRIX, TS_LINEAR_VARIABLE_MATRIX,
               TS_LINEAR_NO_MATRIX, TS_NONLINEAR_JACOBIAN,
               TS_NONLINEAR_NO_JACOBIAN} TSProblemType;

extern int TSCreate(MPI_Comm,TSProblemType,TS*);
extern int TSSetType(TS,TSType);
extern int TSDestroy(TS);

extern int TSSetMonitor(TS,int(*)(TS,int,Scalar,Vec,void*),void *);
extern int TSGetType(TS,TSType*,char**);

extern int TSSetFromOptions(TS);
extern int TSSetUp(TS);

extern int TSSetSolution(TS,Vec);
extern int TSGetSolution(TS,Vec*);

extern int TSSetDuration(TS,int,Scalar);
extern int TSPrintHelp(TS);

extern int TSDefaultMonitor(TS,int,Scalar,Vec,void*);
extern int TSStep(TS,int *,Scalar*);

extern int TSSetInitialTimeStep(TS,Scalar,Scalar);
extern int TSGetTimeStep(TS,Scalar *);
extern int TSSetTimeStep(TS,Scalar);

extern int TSSetRHSFunction(TS,int (*)(TS,Scalar,Vec,Vec,void*),void*);
extern int TSSetRHSMatrix(TS,Mat,Mat,int (*)(TS,Scalar,Mat*,Mat*,MatStructure*,void*),void*);
extern int TSSetRHSJacobian(TS,Mat,Mat,int(*)(TS,Scalar,Vec,Mat*,Mat*,MatStructure*,void*),void*);

extern int TSComputeRHSFunction(TS,Scalar,Vec,Vec);

extern int TSRegisterAll();
extern int TSRegister(int,char*,int (*)(TS));

extern int TSGetSNES(TS,SNES*);
extern int TSGetSLES(TS,SLES*);

#endif