simData.h

Go to the documentation of this file.

#ifndef SIMDATA_H
#define SIMDATA_H

#include <vector>
#include <string>
#include "platformEnv.h"
#include "dataArgs.h"

// Macro for getting array index
#define ID(variable, idx, jdx, kdx) ((variable)*(d->Nx)*(d->Ny)*(d->Nz) + (idx)*(d->Ny)*(d->Nz) + (jdx)*(d->Nz) + (kdx))



class Data
{
  public:
    int
    nx, ny, nz;            

    double
    xmin, xmax,
    ymin, ymax,            
    zmin, zmax,
    endTime,               
    cfl;                   
    int Ng;                
    double
    gamma,                 
    sigma;                 
    int
    memSet,                
    bcsSet,                

    Ncons, Nprims, Naux;   

    double
    cp,                    

    mu1, mu2;              

    int
    frameSkip;             
    int
    reportItersPeriod;     
    bool
    functionalSigma;       
    double
    gam;                   
    double
    *cons, *prims, *aux,
    *f, *fnet,             
    *source,
    *sourceExtension,

    *x, *y, *z;            

    double
    alphaX, alphaY, alphaZ,

    t=-1,                  
    dt,                    

    dx, dy, dz;            

    int
    iters,                 

    Nx, Ny, Nz, Ntot;      

    std::vector<std::string>
    consLabels,
    primsLabels,           
    auxLabels;
    int
    dims,                  

    is, js, ks,
    ie, je, ke;            

    std::vector<double>
    optionalSimArgs;          
    std::vector<std::string>
    optionalSimArgNames;     
    int
    nOptionalSimArgs=0;      



    int id(int var, int i, int j, int k) {
      return var * this->Nx * this->Ny * this->Nz + i * this->Ny * this->Nz + j * this->Nz + k;
    }



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


     void initData(PlatformEnv *env, int nOptionalSimArgs=0, std::vector<double> optionalSimArgs=std::vector<double>(), std::vector<std::string> optionalSimArgNames=std::vector<std::string>());



    Data(int nx, int ny, int nz,
         double xmin, double xmax,
         double ymin, double ymax,
         double zmin, double zmax,
         double endTime,
         PlatformEnv *env,
         double cfl=0.5, int Ng=4,
         double gamma=5.0/3.0,
         double sigma=1e3,
         double cp=0.1,
         double mu1=-1.0e4, double mu2=1.0e4,
         int frameskip=10,
         int reportItersPeriod=1,
         bool funtionalSigma=false, double gam=12);


    Data(DataArgsBase args, PlatformEnv *env);

    ~Data() {};

};

#endif