#ifndef _ADAPT_DT_H_
#define _ADAPT_DT_H_

#ifndef NDIM
#define NDIM 3
#endif

typedef struct {
  Vec Y;
} TSAdapt_MyCFLAdapt;

extern double cfl;
extern double gam;

typedef struct {
    PetscReal r;
    PetscReal ru[NDIM];
    PetscReal rE;
} EulerNode;

PetscErrorCode MyTSAdaptRegister(TS, TSAdapt, double, double);
PetscErrorCode TSAdaptChoose_MyCFLAdapt(TSAdapt, TS, PetscReal,
                                        PetscInt*, PetscReal*,
                                        PetscBool*, PetscReal*,
                                        PetscReal*, PetscReal*);
PetscErrorCode TSAdaptReset_MyCFLAdapt(TSAdapt);
PetscErrorCode TSAdaptDestroy_MyCFLAdapt(TSAdapt);
PetscErrorCode TSAdaptCreate_MyCFLAdapt(TSAdapt);

#endif