28 #if MHD_FORMULATION == DIV_CLEANING
31 MX2 = 255, BX2 = 255, MX3 = 255, BX3 = 255,
35 RHO, MX1, MX2, BX1, BX2,
39 #if MHD_FORMULATION == DIV_CLEANING
46 RHO, MX1, MX2, MX3, BX1, BX2, BX3,
50 #if MHD_FORMULATION == DIV_CLEANING
56 VX1 = MX1, VX2 = MX2, VX3 = MX3,
64 VX = VX1, VY = VX2, VZ = VX3,
65 MX = MX1, MY = MX2, MZ = MX3,
66 BX = BX1, BY = BX2, BZ = BX3
84 #define NFLX (2*COMPONENTS + (EOS == IDEAL ? 2:1) + (MHD_FORMULATION == DIV_CLEANING))
104 #if MHD_FORMULATION != DIV_CLEANING
115 #if MHD_FORMULATION == DIV_CLEANING
116 , KPSI_GLMM, KPSI_GLMP
126 #define AX1 (NVAR + 1)
127 #define AX2 (NVAR + 2)
128 #define AX3 (NVAR + 3)
140 #if GEOMETRY == CYLINDRICAL
156 #if GEOMETRY == POLAR
172 #if GEOMETRY == SPHERICAL
194 #if BACKGROUND_FIELD == YES
195 double **GetBackgroundField (
int,
int,
int,
Grid *);
198 int ConsToPrim (
double **, real **,
int ,
int,
unsigned char *);
199 void Eigenvalues (
double **,
double *,
double **,
int,
int);
201 void PrimEigenvectors (
double *,
double,
double,
double *,
double **,
double **);
203 double **,
double **,
double *);
205 void Enthalpy (
double **,
double *,
int,
int );
206 void Entropy (
double **,
double *,
int,
int );
208 void Flux (
double **,
double **,
double *,
double **,
double **,
210 void HLL_Speed (
double **,
double **,
double *,
double *,
double **,
211 double *,
double *,
int,
int);
212 void MaxSignalSpeed (
double **,
double *,
double *,
double *,
double **,
int,
int);
213 void PrimToCons (
double **,
double **,
int,
int);
214 void PrimRHS (
double *,
double *,
double,
double,
double *);
216 double *,
double *,
double **,
Grid *);
218 #if MHD_FORMULATION == EIGHT_WAVES
220 void ROE_DIVB_SOURCE (
const State_1D *,
int,
int,
Grid *);
221 void HLL_DIVB_SOURCE (
const State_1D *,
double **,
int,
int,
Grid *);
223 #elif MHD_FORMULATION == DIV_CLEANING
227 #elif MHD_FORMULATION == CONSTRAINED_TRANSPORT
233 Riemann_Solver HLL_Solver, HLLC_Solver, HLLD_Solver;
234 Riemann_Solver LF_Solver, Roe_Solver;
236 #if RESISTIVE_MHD != NO