15 #ifndef ResRMHD_RESISTIVITY 16 #define ResRMHD_RESISTIVITY IMEX 41 #if DIVB_CONTROL != CONSTRAINED_TRANSPORT 46 #if DIVB_CONTROL == DIV_CLEANING 47 #define PSI_GLM (CRG + 1) 48 #define PHI_GLM (CRG + 2) 59 #define NFLX (10 + DIV_COMP + HAVE_ENERGY + 1) 61 #define NFLX (10 + DIV_COMP + HAVE_ENERGY) 76 KFASTM, KFASTP, KENTRP
78 #if DIVB_CONTROL != DIV_CLEANING 82 , KSLOWM, KSLOWP, KALFVM, KALFVP
84 #if DIVB_CONTROL == DIV_CLEANING 85 , KPSI_GLMM, KPSI_GLMP
101 typedef struct Map_Param_{
124 #define AX1 (NVAR + 1) 125 #define AX2 (NVAR + 2) 126 #define AX3 (NVAR + 3) 137 void ComputeCurrent(
double *,
double *);
138 void ComputeElectricField(
double *,
double *);
140 int ConsToPrim (
double **,
double **,
int,
int, uint16_t *);
141 void ConvertTo4vel (
double **,
int,
int);
142 void ConvertTo3vel (
double **,
int,
int);
148 double *,
double *SR,
int beg,
int end);
149 void HLL_Speed_RMHD (
const State *,
const State *,
double *,
double *,
int,
int);
153 void MaxSignalSpeed_RHD (
const State *,
double *,
double *,
int,
int);
157 int PressureFix (
double *,
double *,
double *);
158 void PrimToCons (
double **,
double **,
int,
int);
161 Riemann_Solver FluxSplitting_LF_Solver, LF_Solver, HLL_Solver;
162 Riemann_Solver GMUSTA1_Solver, MHLL_Solver, MHLLC_Solver;
184 #ifndef CHARGE_SCHEME 185 #define CHARGE_SCHEME 0 188 #ifndef RMHD_REDUCED_ENERGY 189 #define RMHD_REDUCED_ENERGY NO 194 #ifndef USE_NEWTON_BROYDEN // YES /NO 195 #define USE_NEWTON_BROYDEN YES 198 #if DIVB_CONTROL == EIGHT_WAVES 199 void POWELL_DIVB_SOURCE(
const Sweep *,
int,
int,
Grid *);
201 #elif DIVB_CONTROL == DIV_CLEANING 203 #elif DIVB_CONTROL == CONSTRAINED_TRANSPORT int RMHD_EnergySolve(Map_param *)
Definition: rmhd_energy_solve.c:28
double W
Definition: mod_defs.h:111
double m2
Definition: mod_defs.h:105
void HLL_Speed_RHD(const State *stateL, const State *stateR, double *SL, double *SR, int beg, int end)
Definition: hll_speed.c:68
void VelocityLimiter(double *, double *, double *)
Definition: vel_limiter.c:16
double S2
Definition: mod_defs.h:107
void MaxSignalSpeed(const State *, double *, double *, int, int)
Definition: eigenv.c:34
int ConsToPrim(double **, double **, int, int, uint16_t *)
Definition: mappers.c:88
void HLL_DIVB_SOURCE(const Sweep *, double **, int, int, Grid *)
Definition: source.c:94
int RHD_PressureFix(double *, double *)
Definition: rhd_pressure_fix.c:28
void IMEXImplicitUpdate(Data *data, double dt1, Data_Arr S, Data_Arr Jf, Grid *grid)
Definition: imex_source.c:27
void PrimToCons(double **, double **, int, int)
Definition: mappers.c:26
double * u
Definition: mod_defs.h:114
Definition: structs.h:124
double lor
Definition: mod_defs.h:112
void Flux(const State *, int, int)
Definition: fluxes.c:23
void HLL_Speed(const State *, const State *, double *, double *, int, int)
Definition: hll_speed.c:24
Header file for GLM Divergence Cleaning.
int RHD_EntropySolve(double *, double *)
Definition: rhd_entropy_solve.c:21
double rho
Definition: mod_defs.h:110
Definition: mod_defs.h:101
Header file for Constrained-Transport (CT) module.
int Magnetosonic(double *vp, double cs2, double h, double *lambda)
Definition: eigenv.c:68
double sigma_c
Definition: mod_defs.h:103
Definition: structs.h:514
double D
Definition: mod_defs.h:102
double B2
Definition: mod_defs.h:108
Definition: structs.h:261
int RHD_EnergySolve(double *, double *)
Definition: rhd_energy_solve.c:25
double prs
Definition: mod_defs.h:113
double S
Definition: mod_defs.h:106
void Resistive_eta(double *v, double x1, double x2, double x3, double *J, double *eta)
Definition: res_eta.c:17
double E
Definition: mod_defs.h:104
Definition: structs.h:289