#ifndef __timestepping_h
#define __timestepping_h

/*
CLASS
   TimeInfo
   */
class TimeInfo
{
double t,k,endzeit;
int it,lit,ps,f,outit;
public:
  TimeInfo (double time=0.,double timestep=.01,int lastit=100,
	    double endtime=1., int printstep=1)
  {
    t=time;
    k=timestep;
    endzeit=endtime;
    it=0;
    lit=lastit;
    ps = printstep;
    f=0;
    outit=0;
  }
  set (double time=0.,double timestep=.01,int lastit=100,
	    double endtime=1., int printstep=1)
  {
    t=time;
    k=timestep;
    endzeit=endtime;
    it=0;
    lit=lastit;
    ps = printstep;
    f=0;
    outit=0;
  }
  double& time() {return t;}
  double& timestep() {return k;}
  double& endtime() {return endzeit;}
  int& lastit() {return lit;}
  int& iteration() {return it;}
  int& printstep() {return ps;}
  int& flag() {return f;}
  int& out() {return outit;}
  int reached()
  {
    if ((it<lit)&&(t<endzeit))
	return 0;
    return 1;
  }
};

//////////
template<class PROBLEM, class VECTOR>
void timestepping(PROBLEM& P,VECTOR& u,VECTOR& b,TimeInfo& info)
{
  P.time_initial(u);
  for(;!info.reached();)
  {
    info.iteration()++;
    info.time() += info.timestep();
    P.time_right(b);
    P.time_solve(u,b);
    P.time_post(u,info);
  }
}

#endif