hybrid.h

Go to the documentation of this file.

#ifndef HYBRID_H
#define HYBRID_H

#include "model.h"
#include "srrmhd.h"
#include "srmhd.h"
#include "REGIME.h"
#include "flux.h"


class Hybrid : public Model
{

  public:

    double
    *icons, *iprims, *iaux,                     // Ideal cons, prims and aux. Size is \f$N_{var}*N_x*N_y*N_z\f$.
    *sicons, *siprims, *siaux,                  // Ideal cons, prims and aux. Size is \f$N_{var}\f$.
    *iflux, *rflux,                             // Flux vectors for ideal, resistive and REGIME. Size is \f$N_{cons}*N_x*N_y*N_z\f$.
    *isource, *rsource, *regimeSource,          // Source vectors for ideal, resistive and REGIME. Size is \f$N_{cons}*N_x*N_y*N_z\f$.
    sigmaCrossOver,                             // Centre conductivity of penalty function
    sigmaSpan;                                  // Span of conductivity of penalty function

    bool useREGIME;                             // Should we use REGIME? (Default to true)



    SRRMHD * resistiveModel;                    // Pointer to SRRMHD model

    SRMHD * idealModel;                         // Pointer to SRMHD model

    REGIME * subgridModel = NULL;               // Pointer to REGIME model

    int *mask;                                  // Flag: can we set REGIME source?

    Hybrid();                                   



    Hybrid(Data * data, double sigmaCrossOver=150, double sigmaSpan=50, bool useREGIME=true);


    virtual ~Hybrid();  


    void setupREGIME(FluxMethod * fluxMethod);

  private:

    double idealWeight(double * cons, double * prims, double * aux);


    double idealWeightID(double * cons, double * prims, double * aux, int i, int j, int k);


    bool useResistive(double * cons, double * prims, double * aux);


    void setIdealCPAs(double *cons, double * prims, double * aux);


    void setIdealCPAsAll(double *cons, double * prims, double * aux);


    void setMasks(double * cons, double * prims, double * aux);

  public:

    void sourceTermSingleCell(double *cons, double *prims, double *aux, double *source, int i=-1, int j=-1, int k=-1);


    void sourceTerm(double *cons, double *prims, double *aux, double *source);


    void getPrimitiveVarsSingleCell(double *cons, double *prims, double *aux, int i=-1, int j=-1, int k=-1);


    void getPrimitiveVars(double *cons, double *prims, double *aux);


    void primsToAll(double *cons, double *prims, double *aux);


    void fluxVector(double *cons, double *prims, double *aux, double *f, const int dir);


    void finalise(double *cons, double *prims, double *aux);
};

#endif