Go to the source code of this file.
Functions | |
void | esys_prim_iso_hyd (const Real d, const Real v1, Real eigenvalues[], Real right_eigenmatrix[][4], Real left_eigenmatrix[][4]) |
ISOTHERMAL HYDRO Input: d, v1 = primitive variables Output: eigenvalues[4], right_eigenmatrix[4,4], left_eigenmatrix[4,4];. | |
void | esys_prim_adb_hyd (const Real d, const Real v1, const Real p, Real eigenvalues[], Real right_eigenmatrix[][5], Real left_eigenmatrix[][5]) |
ADIABATIC HYDRO Input: d, v1, p = primitive variables Output: eigenvalues[5], right_eigenmatrix[5,5], left_eigenmatrix[5,5];. | |
void | esys_prim_iso_mhd (const Real d, const Real v1, const Real b1, const Real b2, const Real b3, Real eigenvalues[], Real right_eigenmatrix[][6], Real left_eigenmatrix[][6]) |
ISOTHERMAL MHD Input: d, v1, b1, b2, b3 = density, velocities, and B field Output: eigenvalues[6], right_eigenmatrix[6,6], left_eigenmatrix[6,6];. | |
void | esys_prim_adb_mhd (const Real d, const Real v1, const Real rho_a2, const Real b1, const Real b2, const Real b3, Real eigenvalues[], Real right_eigenmatrix[][7], Real left_eigenmatrix[][7]) |
ADIABATIC MHD Input: d, v1, p, b1, b2, b3 = density, velocities, pressure, and B field Output: eigenvalues[7], right_eigenmatrix[7,7], left_eigenmatrix[7,7];. | |
void | lr_states_destruct (void) |
NoOp for first order, but included for compatibility with integrator (needed for 2nd and 3rd order). | |
void | lr_states_init (MeshS *pM) |
NoOp for first order, but included for compatibility with integrator (needed for 2nd and 3rd order). | |
void | lr_states (const GridS *pG, const Prim1DS W[], const Real Bxc[], const Real dt, const Real dx, const int is, const int ie, Prim1DS Wl[], Prim1DS Wr[], const int dir) |
Computes L/R states Input Arguments: W = PRIMITIVE variables at cell centers along 1-D slice Bxc = B in direction of slice at cell center dtodx = dt/dx il,iu = lower and upper indices of zone centers in slice W and Bxc must be initialized over [il-2:iu+2]. |
void esys_prim_adb_hyd | ( | const Real | d, | |
const Real | v1, | |||
const Real | p, | |||
Real | eigenvalues[], | |||
Real | right_eigenmatrix[][5], | |||
Real | left_eigenmatrix[][5] | |||
) |
ADIABATIC HYDRO Input: d, v1, p = primitive variables Output: eigenvalues[5], right_eigenmatrix[5,5], left_eigenmatrix[5,5];.
Definition at line 120 of file esystem_prim.c.
References a.
Referenced by lr_states().
void esys_prim_adb_mhd | ( | const Real | d, | |
const Real | v1, | |||
const Real | rho_a2, | |||
const Real | b1, | |||
const Real | b2, | |||
const Real | b3, | |||
Real | eigenvalues[], | |||
Real | right_eigenmatrix[][7], | |||
Real | left_eigenmatrix[][7] | |||
) |
ADIABATIC MHD Input: d, v1, p, b1, b2, b3 = density, velocities, pressure, and B field Output: eigenvalues[7], right_eigenmatrix[7,7], left_eigenmatrix[7,7];.
Definition at line 388 of file esystem_prim.c.
Referenced by lr_states().
void esys_prim_iso_hyd | ( | const Real | d, | |
const Real | v1, | |||
Real | eigenvalues[], | |||
Real | right_eigenmatrix[][4], | |||
Real | left_eigenmatrix[][4] | |||
) |
ISOTHERMAL HYDRO Input: d, v1 = primitive variables Output: eigenvalues[4], right_eigenmatrix[4,4], left_eigenmatrix[4,4];.
Definition at line 51 of file esystem_prim.c.
References Iso_csound.
Referenced by lr_states().
void esys_prim_iso_mhd | ( | const Real | d, | |
const Real | v1, | |||
const Real | b1, | |||
const Real | b2, | |||
const Real | b3, | |||
Real | eigenvalues[], | |||
Real | right_eigenmatrix[][6], | |||
Real | left_eigenmatrix[][6] | |||
) |
ISOTHERMAL MHD Input: d, v1, b1, b2, b3 = density, velocities, and B field Output: eigenvalues[6], right_eigenmatrix[6,6], left_eigenmatrix[6,6];.
Definition at line 212 of file esystem_prim.c.
References Iso_csound, Iso_csound2, and s.
Referenced by lr_states().
void lr_states | ( | const GridS * | pG, | |
const Prim1DS | W[], | |||
const Real | Bxc[], | |||
const Real | dt, | |||
const Real | dx, | |||
const int | il, | |||
const int | iu, | |||
Prim1DS | Wl[], | |||
Prim1DS | Wr[], | |||
const int | dir | |||
) |
Computes L/R states Input Arguments: W = PRIMITIVE variables at cell centers along 1-D slice Bxc = B in direction of slice at cell center dtodx = dt/dx il,iu = lower and upper indices of zone centers in slice W and Bxc must be initialized over [il-2:iu+2].
Computes L/R states Input Arguments:
Computes L/R states.
Output Arguments: Wl,Wr = L/R-states of PRIMITIVE variables at interfaces over [il:iu+1]
Input Arguments:
Output Arguments:
Output Arguments:
Definition at line 62 of file lr_states_plm.c.
References beta, Bz(), C, d(), dWm, esys_prim_adb_hyd(), esys_prim_adb_mhd(), esys_prim_iso_hyd(), esys_prim_iso_mhd(), Gamma, n, omg, pW, GridS::r, r, GridS::ri, ri, vel, Prim1DS::Vx, Whalf, and Wim1h.
void lr_states_destruct | ( | void | ) |
NoOp for first order, but included for compatibility with integrator (needed for 2nd and 3rd order).
Free memory used by work arrays.
Definition at line 73 of file lr_states_dc.c.
References dWm, free_2d_array(), pW, vel, Whalf, and Wim1h.
Referenced by lr_states_init(), and main().
void lr_states_init | ( | MeshS * | pM | ) |
NoOp for first order, but included for compatibility with integrator (needed for 2nd and 3rd order).
Allocate enough memory for work arrays.
Definition at line 63 of file lr_states_dc.c.
References ath_error(), calloc_2d_array(), MeshS::Domain, MeshS::DomainsPerLevel, dWm, DomainS::Grid, lr_states_destruct(), MeshS::NLevels, GridS::Nx, pW, vel, Whalf, and Wim1h.
Referenced by main().