timeInt.h

Go to the documentation of this file.

#ifndef TIMEINT_H
#define TIMEINT_H

#include "simData.h"
#include "model.h"
#include "boundaryConds.h"
#include "flux.h"
#include "modelExtension.h"
#include "hybrid.h"


class TimeIntegratorBase
{
  public:

    Data * data;                

    Model * model;              

    Bcs * bcs;                  

    FluxMethod * fluxMethod;    

    ModelExtension * modelExtension;      


    TimeIntegratorBase(Data * data, Model * model, Bcs * bcs, FluxMethod * fluxMethod, ModelExtension * modelExtension = NULL) :
                   data(data), model(model), bcs(bcs), fluxMethod(fluxMethod), modelExtension(modelExtension) { }


    virtual ~TimeIntegratorBase() { }     


    virtual void step(double * cons, double * prims, double * aux, double dt=0) = 0;



    virtual void finalise(double * cons, double * prims, double * aux) = 0;

};



class TimeIntegrator : public TimeIntegratorBase
{
  public:


    TimeIntegrator(Data * data, Model * model, Bcs * bcs, FluxMethod * fluxMethod, ModelExtension * modelExtension = NULL) :
                    TimeIntegratorBase(data, model, bcs, fluxMethod, modelExtension) { }


    virtual ~TimeIntegrator() { }     


    void finalise(double * cons, double * prims, double * aux)
    {
      // Perfrom C2P
      try {
        this->model->getPrimitiveVars(cons, prims, aux);
      }
      catch (const std::exception& e) {
        printf("Time integration raises exception with following message:\n%s\n", e.what());
        throw e;
      }

      // Finalise via model
      model->finalise(cons, prims, aux);
      // Apply BCs
      this->bcs->apply(cons, prims, aux);
    }

};

#endif