reconstruction/esystem_prim.c File Reference

Functions to evaluate the eigenvalues, and left- and right-eigenvectors for the linearized system in the PRIMITIVE variables. More...

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 rho_a2, 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];.

---

Detailed Description

Functions to evaluate the eigenvalues, and left- and right-eigenvectors for the linearized system in the PRIMITIVE variables.

PURPOSE: Functions to evaluate the eigenvalues, and left- and right-eigenvectors for the linearized system in the PRIMITIVE variables, i.e. W,t = AW,x, where W=(d,vx,vy,vz,[P],[By,Bz]). The eigenvalues are returned through the argument list as a vector of length NWAVE. The eigenvectors are returned as matrices of size (NWAVE)x(NWAVE), with right-eigenvectors stored as COLUMNS (R_i[*] = right_eigenmatrix[*][i]), and left-eigenvectors stored as ROWS (L_i[*] = left_eigenmatrix[i][*]).

To improve performance components of the eigenvectors which are zero are not set here (eigenmatrices must be initialized to zero in calling routine). However, for completeness statements which set these values are included, but are commented out.

REFERENCES:

• J. Stone, T. Gardiner, P. Teuben, J. Hawley, & J. Simon "Athena: A new code for astrophysical MHD", ApJS, (2008), Appendix A Equation numbers refer to this paper.

CONTAINS PUBLIC FUNCTIONS:

• esys_prim_iso_hyd() - isothermal hydrodynamics
• esys_prim_adb_hyd() - adiabatic hydrodynamics
• esys_prim_iso_mhd() - isothermal MHD
• esys_prim_adb_mhd() - adiabatic MHD ============================================================================

Definition in file esystem_prim.c.

---

Function Documentation

void esys_prim_adb_hyd	(	const Real	d,
		const Real	v1,
		const Real	rho_a2,
		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().

Here is the caller graph for this function:

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.

References a, and s.

Referenced by lr_states().

Here is the caller graph for this function:

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().

Here is the caller graph for this function:

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().

Here is the caller graph for this function: