23 void Flux (
double **ucons,
double **wprim,
double *a2,
double **bck,
24 double **fx,
double *p,
int beg,
int end)
44 for (i = beg; i <= end; i++) {
53 #elif EOS == BAROTROPIC
54 ptot += BAROTROPIC_PR(w[
RHO]);
55 #elif EOS == ISOTHERMAL
58 print (
"! Flux(): not defined for this EoS\n");
62 vB = EXPAND(w[
VX1]*w[BX1] , + w[
VX2]*w[BX2], + w[
VX3]*w[BX3]);
64 #if BACKGROUND_FIELD == YES
65 ptot += EXPAND(bck[i][BX1]*w[BX1], + bck[i][BX2]*w[BX2], + bck[i][BX3]*w[BX3]);
67 EXPAND(bt1 = w[
BXn] + bck[i][
BXn]; ,
72 EXPAND(fx[i][
MX1] = w[
VXn]*u[
MX1] - bt1*w[BX1] - w[
BXn]*bck[i][BX1]; ,
76 EXPAND(fx[i][
BXn] = 0.0; ,
80 fx[i][ENG] = (u[ENG] + ptot)*w[
VXn] - bt1*vB;
88 EXPAND(fx[i][
BXn] = 0.0; ,
92 fx[i][ENG] = (u[ENG] + ptot)*w[
VXn] - w[
BXn]*vB;
void Flux(double **u, double **w, double *a2, double **fx, double *p, int beg, int end)
void print(const char *fmt,...)
double glm_ch
The propagation speed of divergence error.
#define QUIT_PLUTO(e_code)