16 #ifndef RESISTIVE_RMHD
17 #define RESISTIVE_RMHD NO
28 RHO, MX1, BX1, ENG, PRS = ENG,
29 #if MHD_FORMULATION == DIV_CLEANING
32 MX2 = 255, BX2 = 255, MX3 = 255, BX3 = 255,
36 RHO, MX1, MX2, BX1, BX2, ENG, PRS = ENG,
37 #if MHD_FORMULATION == DIV_CLEANING
44 RHO, MX1, MX2, MX3, BX1, BX2, BX3, ENG, PRS = ENG,
45 #if MHD_FORMULATION == DIV_CLEANING
51 VX1 = MX1, VX2 = MX2, VX3 = MX3,
55 DN = RHO, PR = PRS, EN = ENG,
56 VX = VX1, VY = VX2, VZ = VX3,
57 MX = MX1, MY = MX2, MZ = MX3,
58 BX = BX1, BY = BX2, BZ = BX3
75 #define NFLX (2 + 2*COMPONENTS + (MHD_FORMULATION==DIV_CLEANING))
89 KFASTM, KFASTP, KENTRP
91 #if MHD_FORMULATION != DIV_CLEANING
102 #if MHD_FORMULATION == DIV_CLEANING
103 , KPSI_GLMM, KPSI_GLMP
108 #define SUBTRACT_DENSITY YES
123 typedef struct MAP_PARAM{
145 #define AX1 (NVAR + 1)
146 #define AX2 (NVAR + 2)
147 #define AX3 (NVAR + 3)
158 #if GEOMETRY == CYLINDRICAL
174 #if GEOMETRY == POLAR
190 #if GEOMETRY == SPHERICAL
208 int ConsToPrim (
double **,
double **,
int,
int,
unsigned char *);
209 void PRIM_EIGENVECTORS (
double *,
double,
double,
double *,
double **,
double **);
211 void Enthalpy (
double **,
double *,
int,
int);
212 void Entropy (
double **,
double *,
int,
int);
225 void Flux (
double **,
double **,
double *,
double **,
double *,
int,
int);
226 void HLL_Speed (
double **,
double **,
double *,
double *,
double *,
double *,
227 double *,
double *,
int,
int);
228 int MaxSignalSpeed (
double **,
double *,
double *,
double *,
double *,
int,
int);
230 void PrimToCons (
double **,
double **,
int,
int);
231 void VelocityLimiter (
double *,
double *,
double *);
233 int Magnetosonic (
double *vp,
double cs2,
double h,
double *lambda);
234 int QuarticSolve (
double,
double,
double,
double,
double *);
235 int CubicSolve (
double,
double,
double,
double *);
237 Riemann_Solver LF_Solver, HLL_Solver, HLLC_Solver, HLLD_Solver;
239 #if MHD_FORMULATION == EIGHT_WAVES
241 void POWELL_DIVB_SOURCE(
const State_1D *,
int,
int,
Grid *);
242 void HLL_DIVB_SOURCE (
const State_1D *,
double **,
int,
int,
Grid *);
244 #elif MHD_FORMULATION == DIV_CLEANING
248 #elif MHD_FORMULATION == CONSTRAINED_TRANSPORT