BDS Namespace Reference

Functions
void	ComputeEdgeState (amrex::Box const &bx, int ncomp, amrex::Array4< amrex::Real const > const &q, AMREX_D_DECL(amrex::Array4< amrex::Real > const &xedge, amrex::Array4< amrex::Real > const &yedge, amrex::Array4< amrex::Real > const &zedge), AMREX_D_DECL(amrex::Array4< amrex::Real const > const &umac, amrex::Array4< amrex::Real const > const &vmac, amrex::Array4< amrex::Real const > const &wmac), amrex::Array4< amrex::Real const > const &divu, amrex::Array4< amrex::Real const > const &fq, amrex::Geometry geom, amrex::Real l_dt, amrex::Vector< amrex::BCRec > const &h_bcrec, amrex::BCRec const *pbc, int const *iconserv, bool is_velocity)
void	ComputeSlopes (amrex::Box const &bx, const amrex::Geometry &geom, int icomp, amrex::Array4< amrex::Real const > const &s, amrex::Array4< amrex::Real > const &slopes, amrex::Vector< amrex::BCRec > const &h_bcrec)
void	ComputeConc (amrex::Box const &bx, const amrex::Geometry &geom, int icomp, amrex::Array4< amrex::Real const > const &s, AMREX_D_DECL(amrex::Array4< amrex::Real > const &sedgex, amrex::Array4< amrex::Real > const &sedgey, amrex::Array4< amrex::Real > const &sedgez), amrex::Array4< amrex::Real const > const &slopes, AMREX_D_DECL(amrex::Array4< amrex::Real const > const &umac, amrex::Array4< amrex::Real const > const &vmac, amrex::Array4< amrex::Real const > const &wmac), amrex::Array4< amrex::Real const > const &divu, amrex::Array4< amrex::Real const > const &force, int const *iconserv, amrex::Real dt, amrex::Vector< amrex::BCRec > const &h_bcrec, amrex::BCRec const *pbc, bool is_velocity)

Detailed Description

Collection of routines for the BDS (Bell-Dawson-Shubin) algorithm.

Function Documentation

ComputeConc()

void BDS::ComputeConc	(	amrex::Box const &	bx,
		const amrex::Geometry &	geom,
		int	icomp,
		amrex::Array4< amrex::Real const > const &	s,
		AMREX_D_DECL(amrex::Array4< amrex::Real > const &sedgex, amrex::Array4< amrex::Real > const &sedgey, amrex::Array4< amrex::Real > const &sedgez)	,
		amrex::Array4< amrex::Real const > const &	slopes,
		AMREX_D_DECL(amrex::Array4< amrex::Real const > const &umac, amrex::Array4< amrex::Real const > const &vmac, amrex::Array4< amrex::Real const > const &wmac)	,
		amrex::Array4< amrex::Real const > const &	divu,
		amrex::Array4< amrex::Real const > const &	force,
		int const *	iconserv,
		amrex::Real	dt,
		amrex::Vector< amrex::BCRec > const &	h_bcrec,
		amrex::BCRec const *	pbc,
		bool	is_velocity
	)

Compute Conc for BDS algorithm.

Parameters

[in]	bx	Current grid patch
[in]	geom	Level geometry.
[in]	icomp	Component of the Array4s.
[in]	s	Array4 of state.
[in,out]	sedgex	Array4 containing x-edges.
[in,out]	sedgey	Array4 containing y-edges.
[in,out]	sedgez	Array4 containing z-edges.
[in]	slopes	Array4 containing slope information.
[in]	umac	Array4 for u-face velocity.
[in]	vmac	Array4 for v-face velocity.
[in]	wmac	Array4 for z-face velocity.
[in]	force	Array4 for forces.
[in]	iconserv	If true, use conservative form, otherwise use convective.
[in]	dt	Time step.
[in]	h_bcrec	Boundary conditions (host).
[in]	pbc	Boundary conditions (device).
[in]	is_velocity	Indicates a component is velocity so boundary conditions can be properly addressed. The header hydro_constants.H defines the component position by [XYZ]VEL macro.

ComputeEdgeState()

void BDS::ComputeEdgeState	(	amrex::Box const &	bx,
		int	ncomp,
		amrex::Array4< amrex::Real const > const &	q,
		AMREX_D_DECL(amrex::Array4< amrex::Real > const &xedge, amrex::Array4< amrex::Real > const &yedge, amrex::Array4< amrex::Real > const &zedge)	,
		AMREX_D_DECL(amrex::Array4< amrex::Real const > const &umac, amrex::Array4< amrex::Real const > const &vmac, amrex::Array4< amrex::Real const > const &wmac)	,
		amrex::Array4< amrex::Real const > const &	divu,
		amrex::Array4< amrex::Real const > const &	fq,
		amrex::Geometry	geom,
		amrex::Real	l_dt,
		amrex::Vector< amrex::BCRec > const &	h_bcrec,
		amrex::BCRec const *	pbc,
		int const *	iconserv,
		bool	is_velocity
	)

Uses the Bell-Dawson-Shubin (BDS) algorithm, a higher order Godunov method for scalar conservation laws in three dimensions, to compute edge states.

Parameters

[in]	bx	Current grid patch
[in]	ncomp	Number of components to work on
[in]	q	Array4 of state, starting at component of interest
[in,out]	xedge	Array4 containing x-edges, starting at component of interest
[in,out]	yedge	Array4 containing y-edges, starting at component of interest
[in,out]	zedge	Array4 containing z-edges, starting at component of interest
[in]	umac	x-Face velocities.
[in]	vmac	y-Face velocities.
[in]	wmac	z-Face velocities.
[in]	fq	Array4 for forces, starting at component of interest
[in]	geom	Level geometry.
[in]	l_dt	Time step.
[in]	h_bcrec	Boundary conditions (host).
[in]	pbc	Boundary conditions (device).
[in]	iconserv	If true, use conservative form, otherwise use convective.
[in]	is_velocity	Indicates a component is velocity so boundary conditions can be properly addressed. The header hydro_constants.H defines the component position by [XYZ]VEL macro.

ComputeSlopes()

void BDS::ComputeSlopes	(	amrex::Box const &	bx,
		const amrex::Geometry &	geom,
		int	icomp,
		amrex::Array4< amrex::Real const > const &	s,
		amrex::Array4< amrex::Real > const &	slopes,
		amrex::Vector< amrex::BCRec > const &	h_bcrec
	)

Compute bilinear slopes for BDS algorithm.

Parameters

[in]	bx	Current grid patch
[in]	geom	Level geometry.
[in]	icomp	Component of the state Array4.
[in]	s	Array4<const> of state vector.
[out]	slopes	Array4 to store slope information.
[in]	h_bcrec	Boundary conditions (host).