REGIME.h

Go to the documentation of this file.

#ifndef REGIME_H
#define REGIME_H

#include "modelExtension.h"
#include "flux.h"

// Macros for accessing matricies
// dwdsb
#define IDWS(ldx, mdx, idx, jdx, kdx)  ((ldx)*(3)*(d->Nx)*(d->Ny)*(d->Nz) + (mdx)*(d->Nx)*(d->Ny)*(d->Nz) + (idx)*(d->Ny)*(d->Nz) + (jdx)*(d->Nz) + (kdx))
// dfxdw, dfydw, dfzdw
#define IDFW(ldx, mdx, idx, jdx, kdx)  ((ldx)*(12)*(d->Nx)*(d->Ny)*(d->Nz) + (mdx)*(d->Nx)*(d->Ny)*(d->Nz) + (idx)*(d->Ny)*(d->Nz) + (jdx)*(d->Nz) + (kdx))
// Mx, My, and Mz matrix
#define IDM(ldx, mdx, idx, jdx, kdx)  ((ldx)*(3)*(d->Nx)*(d->Ny)*(d->Nz) + (mdx)*(d->Nx)*(d->Ny)*(d->Nz) + (idx)*(d->Ny)*(d->Nz) + (jdx)*(d->Nz) + (kdx))



class REGIME : public ModelExtension
{
  public:

    double
    //{
    *dfxdw, *dfydw, *dfzdw,      
    //}
    *dwdsb,                      
    *E,                          
    *q,                          
    *K,                          
    //{
    *Mx, *My, *Mz,               
    //}
    //{
    *fbx, *fby, *fbz,               
    //}
    //{
    *diffuX, *diffuY, *diffuZ,   
    //}
    *alpha;                      

    FluxMethod * fluxMethod;     

    REGIME();

    REGIME(Data * data, FluxMethod * fluxMethod);

    virtual ~REGIME();


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

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

    //{
    void set_Dx(double * cons, double * prims, double * aux);
    void set_Dy(double * cons, double * prims, double * aux);
    void set_Dz(double * cons, double * prims, double * aux);
    //}


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

    //{
    void set_dfxdw(double * cons, double * prims, double * aux);
    void set_dfydw(double * cons, double * prims, double * aux);
    void set_dfzdw(double * cons, double * prims, double * aux);
    //}

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

};

#endif