Athena: integrators/integrate_1d_vl_sr.c File Reference

Integrate SRMHD equations using 1D version of MUSCL-Hancock (VL) integrator. More...

Go to the source code of this file.

Functions
static void	FixCell (GridS *pG, Int3Vect)
	Uses first order fluxes to fix negative d,P or superluminal v.
void	integrate_1d_vl (DomainS *pD)
	1D version of van Leer unsplit integrator for MHD.
void	integrate_init_1d (MeshS *pM)
	Allocate temporary integration arrays.
void	integrate_destruct_1d (void)
	Free temporary integration arrays.
Variables
static Prim1DS *	Wl_x1Face = NULL
static Prim1DS *	Wr_x1Face = NULL
static Cons1DS *	x1Flux = NULL
static Real *	Bxc = NULL
static Real *	Bxi = NULL
static Prim1DS *	W1d = NULL
static Prim1DS *	Wl = NULL
static Prim1DS *	Wr = NULL
static Cons1DS *	U1d = NULL
static Cons1DS *	Ul = NULL
static Cons1DS *	Ur = NULL
static Cons1DS *	x1FluxP = NULL
static PrimS *	W = NULL
static ConsS *	Uhalf = NULL
static PrimS *	Whalf = NULL

Detailed Description

Integrate SRMHD equations using 1D version of MUSCL-Hancock (VL) integrator.

PURPOSE: Integrate SRMHD equations using 1D version of MUSCL-Hancock (VL) integrator. Updates U.[d,M1,M2,M3,E,B2c,B3c,s] in Grid structure. Adds gravitational source terms, self-gravity.

CONTAINS PUBLIC FUNCTIONS:

Definition in file integrate_1d_vl_sr.c.

Function Documentation

static void FixCell	(	GridS *	pG,
		Int3Vect	indx
	)			`[static]`

Uses first order fluxes to fix negative d,P or superluminal v.

Definition at line 616 of file integrate_1d_vl_sr.c.

References ConsS::B2c, GridS::B2i, ConsS::B3c, GridS::B3i, Cons1DS::By, Cons1DS::Bz, ConsS::d, Cons1DS::d, GridS::dt, GridS::dx1, ConsS::E, Cons1DS::E, Int3Vect::i, GridS::ie, GridS::is, GridS::js, GridS::ks, ConsS::M1, ConsS::M2, ConsS::M3, Cons1DS::Mx, Cons1DS::My, Cons1DS::Mz, and GridS::U.

Referenced by integrate_1d_vl().

Here is the caller graph for this function:

void integrate_1d_vl ( DomainS * pD )

1D version of van Leer unsplit integrator for MHD.

The numbering of steps follows the numbering in the 3D version. NOT ALL STEPS ARE NEEDED IN 1D.

Definition at line 63 of file integrate_1d_vl_sr.c.

References PrimS::B1c, ConsS::B1c, GridS::B1i, PrimS::B2c, ConsS::B2c, GridS::B2i, PrimS::B3c, ConsS::B3c, GridS::B3i, Bxc, Bxi, Prim1DS::By, Cons1DS::By, Prim1DS::Bz, Cons1DS::Bz, cc_pos(), GridS::CGrid, check_Prim(), Cons1D_to_Prim1D(), Cons_to_Prim(), PrimS::d, Prim1DS::d, d(), ConsS::d, Cons1DS::d, dim, GridS::dt, GridS::dx1, ConsS::E, Cons1DS::E, fix_vsq(), FixCell(), fluxes(), four_pi_G, grav_mean_rho, DomainS::Grid, Int3Vect::i, GridS::ie, GridOvrlpS::ijke, GridOvrlpS::ijks, il, GridS::is, iu, Int3Vect::j, GridS::js, Int3Vect::k, GridS::ks, lr_states(), ConsS::M1, ConsS::M2, ConsS::M3, Cons1DS::Mx, Cons1DS::My, GridOvrlpS::myFlx, Cons1DS::Mz, n, GridS::NCGrid, GridS::NPGrid, PrimS::P, Prim1DS::P, GridS::PGrid, GridS::Phi, Prim1D_to_Cons1D(), Prim_to_Cons(), Cons1DS::s, ConsS::s, s, StaticGravPot, U, GridS::U, Ul, Ur, PrimS::V1, PrimS::V2, PrimS::V3, Prim1DS::Vx, Prim1DS::Vy, Vz(), Prim1DS::Vz, Wl, Wr, Wr_x1Face, and GridS::x1MassFlux.

Here is the call graph for this function:

void integrate_destruct_1d ( void )

Free temporary integration arrays.

Definition at line 584 of file integrate_1d_vl_sr.c.

References Bxc, Bxi, dhalf, free_1d_array(), geom_src, phalf, Ul, Ur, Ur_x1Face, Wl, Wr, Wr_x1Face, and x1Flux.

Here is the call graph for this function:

void integrate_init_1d ( MeshS * pM )

Allocate temporary integration arrays.

Definition at line 536 of file integrate_1d_vl_sr.c.

References ath_error(), Bxc, Bxi, calloc_1d_array(), CoolingFunc, dhalf, MeshS::Domain, MeshS::DomainsPerLevel, geom_src, DomainS::Grid, integrate_destruct(), MeshS::NLevels, GridS::Nx, phalf, StaticGravPot, Ul, Ur, Ur_x1Face, Wl, Wr, Wr_x1Face, and x1Flux.

Here is the call graph for this function:

Variable Documentation

Real* Bxc = NULL [static]

Definition at line 36 of file integrate_1d_vl_sr.c.

Referenced by integrate_1d_vl(), integrate_destruct_1d(), and integrate_init_1d().

Real * Bxi = NULL [static]

Definition at line 36 of file integrate_1d_vl_sr.c.

Referenced by integrate_1d_vl(), integrate_destruct_1d(), and integrate_init_1d().

Cons1DS* U1d = NULL [static]

Definition at line 38 of file integrate_1d_vl_sr.c.

ConsS* Uhalf = NULL [static]

Definition at line 47 of file integrate_1d_vl_sr.c.

Cons1DS * Ul = NULL [static]

Definition at line 38 of file integrate_1d_vl_sr.c.

Referenced by integrate_1d_vl(), integrate_destruct_1d(), and integrate_init_1d().

Cons1DS * Ur = NULL [static]

Definition at line 38 of file integrate_1d_vl_sr.c.

Referenced by integrate_1d_vl(), integrate_destruct_1d(), and integrate_init_1d().

PrimS* W = NULL [static]

Definition at line 44 of file integrate_1d_vl_sr.c.

Prim1DS* W1d = NULL [static]

Definition at line 37 of file integrate_1d_vl_sr.c.

PrimS* Whalf = NULL [static]

Definition at line 48 of file integrate_1d_vl_sr.c.

Referenced by integrate_2d_vl(), and integrate_3d_vl().

Prim1DS * Wl = NULL [static]

Definition at line 37 of file integrate_1d_vl_sr.c.

Referenced by integrate_1d_vl(), integrate_destruct_1d(), and integrate_init_1d().

Prim1DS* Wl_x1Face = NULL [static]

Definition at line 32 of file integrate_1d_vl_sr.c.

Prim1DS * Wr = NULL [static]

Definition at line 37 of file integrate_1d_vl_sr.c.

Referenced by integrate_1d_vl(), integrate_destruct_1d(), and integrate_init_1d().

Prim1DS * Wr_x1Face = NULL [static]

Definition at line 32 of file integrate_1d_vl_sr.c.

Referenced by integrate_1d_vl(), integrate_destruct_1d(), and integrate_init_1d().

Cons1DS* x1Flux = NULL [static]

Definition at line 33 of file integrate_1d_vl_sr.c.

Cons1DS* x1FluxP = NULL [static]

Definition at line 40 of file integrate_1d_vl_sr.c.