srmhd.h

Go to the documentation of this file.

#ifndef SRMHD_H
#define SRMHD_H

#include "model.h"


/*
This is the human readable description of this models variables.

  SRMHD has nine conserved variables:
    D, Sx, Sy, Sz, tau, Bx, By, Bz, phi
  Eight primitive variables:
    rho, vx, vy, vz, p, Bx, By, Bz
  Thirteen auxiliary variables:
    h, W, e, c, b0, bx, by, bz, bsq, vsq, BS, Bsq, Ssq
*/




class SRMHD : public Model
{

  public:

    int smartGuesses;     

    double * solution;    


    SRMHD();     


    SRMHD(Data * data);

    virtual ~SRMHD();     



    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) { };

};




int residual(void *p, int n, const double *x, double *fvec, int iflag);




typedef struct
{
  double
  D,    
  g,    
  Bsq,  
  Ssq,  
  BS,   
  tau;  
  int i;
} Args;



typedef struct
{
  // Store coordinates of the failed cell
  int
  x, y, z;  

} Failed;

#endif