Variables
	parser = exahype2.ArgumentParser("ExaHyPE 2 - Finite Volumes Coupling Testing Script")
	min_depth
	degrees_of_freedom
	end_time
	number_of_snapshots
	periodic_boundary_conditions_x
	periodic_boundary_conditions_y
	periodic_boundary_conditions_z
	args = parser.parse_args()
float	max_h = 1.1 / 3.0**args.min_depth
float	min_h = max_h * 3.0 ** (-args.amr_levels)
dict	euler_unknowns
	EULER #.
str	euler_initial_conditions
str	euler_boundary_conditions
str	compute_primitive_variables
str	euler_max_eigenvalue
str	euler_flux
	euler_solver
	initial_conditions
	boundary_conditions
	flux
	max_eigenvalue
dict	advection_unknowns = {"scalar": 1}
	ADVECTION #.
dict	advection_auxiliary_variables
str	advection_initial_conditions
str	advection_boundary_conditions
str	advection_flux
str	advection_max_eigenvalue
	advection_solver
	_compute_time_step_size
	COUPLING #.
	number_of_variables_euler = euler_solver.auxiliary_variables + euler_solver.unknowns
tuple	number_of_variables_advection
	last_solver
	project
	PROJECT #.
float	time_in_between_plots = 0.0
	dimensions
	size
	offset
	min_end_time
	max_end_time
	first_plot_time_stamp
	periodic_BC
	mode
	make
	True
	make_clean_first
	throw_away_data_after_build

Variable Documentation

◆ _compute_time_step_size

fv-coupling._compute_time_step_size

protected

COUPLING #.

Definition at line 194 of file fv-coupling.py.

◆ advection_auxiliary_variables

dict fv-coupling.advection_auxiliary_variables

Initial value:

=  {
    "rho": 1,
    "rhoU": 1,
    "rhoV": 1,
}

Definition at line 126 of file fv-coupling.py.

◆ advection_boundary_conditions

str fv-coupling.advection_boundary_conditions

Initial value:

=  """
  // Reflective boundary conditions
  Qoutside[Shortcuts::scalar] = Qinside[Shortcuts::scalar];
  Qoutside[Shortcuts::rho]    = Qinside[Shortcuts::rho];
  Qoutside[Shortcuts::rhoU]   = -Qinside[Shortcuts::rhoU];
  Qoutside[Shortcuts::rhoV]   = -Qinside[Shortcuts::rhoV];
#if DIMENSIONS == 3
  Qoutside[Shortcuts::rhoW]   = -Qinside[Shortcuts::rhoW];
#endif
"""

Definition at line 144 of file fv-coupling.py.

◆ advection_flux

str fv-coupling.advection_flux

Initial value:

=  """
  const auto irho = 1.0 / Q[Shortcuts::rho];
  const auto u_n = Q[Shortcuts::rhoU + normal] * irho;
  F[Shortcuts::scalar] = Q[Shortcuts::scalar] * u_n;
"""

Definition at line 155 of file fv-coupling.py.

◆ advection_initial_conditions

str fv-coupling.advection_initial_conditions

Initial value:

=  """
  Q[Shortcuts::scalar]  = ((x(0) < 0.5) ? (50) : 0.0);
  Q[Shortcuts::rho]     = 0.1;
  Q[Shortcuts::rhoU]    = 0.0;
  Q[Shortcuts::rhoV]    = 0.0;
#if DIMENSIONS == 3
  Q[Shortcuts::rhoW]    = 0.0;
#endif
"""

Definition at line 134 of file fv-coupling.py.

◆ advection_max_eigenvalue

str fv-coupling.advection_max_eigenvalue

Initial value:

=  """
  if (dt == 0) {
    // Return dummy value initially because all velocities are 0
    return 1.0;
  }
  // We could also return a dummy value here but this way, mistakes are more obvious.
  const auto irho = 1.0 / Q[Shortcuts::rho];
  const auto u_n = Q[Shortcuts::rhoU + normal] * irho;
  return u_n;
"""

Definition at line 161 of file fv-coupling.py.

◆ advection_solver

fv-coupling.advection_solver

Initial value:

=  exahype2.solvers.fv.godunov.GlobalAdaptiveTimeStep(
    name="AdvectionFVSolver",
    patch_size=args.degrees_of_freedom,
    unknowns=advection_unknowns,
    auxiliary_variables=advection_auxiliary_variables,
    min_volume_h=min_h,
    max_volume_h=max_h,
    time_step_relaxation=0.5,
    solvers_before_in_coupling=[euler_solver],
)

Definition at line 172 of file fv-coupling.py.

◆ advection_unknowns

dict fv-coupling.advection_unknowns = {"scalar": 1}

ADVECTION #.

Definition at line 125 of file fv-coupling.py.

◆ args

fv-coupling.args = parser.parse_args()

Definition at line 17 of file fv-coupling.py.

◆ boundary_conditions

fv-coupling.boundary_conditions

Definition at line 116 of file fv-coupling.py.

◆ compute_primitive_variables

str fv-coupling.compute_primitive_variables

Initial value:

=  r"""
  const auto irho = 1.0 / Q[Shortcuts::rho];
  const auto u0 = Q[Shortcuts::rhoU + 0] * irho;
  const auto u1 = Q[Shortcuts::rhoU + 1] * irho;
#if DIMENSIONS == 3
  const auto u2 = Q[Shortcuts::rhoU + 2] * irho;
#endif
 
#if DIMENSIONS == 3
  const auto uSq = u0 * u0 + u1 * u1 + u2 * u2;
  const auto u_n = (normal == 0) ? u0 : (normal == 1) ? u1 : u2;
#else
  const auto uSq = u0 * u0 + u1 * u1;
  const auto u_n = (normal == 0) ? u0 : u1;
#endif
 
  const auto internalE = Q[Shortcuts::rhoE] - 0.5 * Q[Shortcuts::rho] * uSq;
  const auto p = (GAMMA - 1.0) * internalE;
"""

Definition at line 56 of file fv-coupling.py.

◆ degrees_of_freedom

fv-coupling.degrees_of_freedom

Definition at line 10 of file fv-coupling.py.

◆ dimensions

fv-coupling.dimensions

Definition at line 288 of file fv-coupling.py.

◆ end_time

fv-coupling.end_time

Definition at line 11 of file fv-coupling.py.

◆ euler_boundary_conditions

str fv-coupling.euler_boundary_conditions

Initial value:

=  """
  // Reflective boundary conditions
  Qoutside[Shortcuts::rho]  = Qinside[Shortcuts::rho];
  Qoutside[Shortcuts::rhoU] = -Qinside[Shortcuts::rhoU];
  Qoutside[Shortcuts::rhoV] = -Qinside[Shortcuts::rhoV];
#if DIMENSIONS == 3
  Qoutside[Shortcuts::rhoW] = -Qinside[Shortcuts::rhoW];
#endif
  Qoutside[Shortcuts::rhoE] = Qinside[Shortcuts::rhoE];
"""

Definition at line 45 of file fv-coupling.py.

◆ euler_flux

str fv-coupling.euler_flux

Initial value:

=  r"""
  {compute_primitive_variables}
 
  F[Shortcuts::rho] = Q[Shortcuts::rhoU + normal];
 
  F[Shortcuts::rhoU + 0] = Q[Shortcuts::rhoU + 0] * u_n;
  F[Shortcuts::rhoU + 1] = Q[Shortcuts::rhoU + 1] * u_n;
#if DIMENSIONS == 3
  F[Shortcuts::rhoU + 2] = Q[Shortcuts::rhoU + 2] * u_n;
#endif
 
  F[Shortcuts::rhoU + normal] += p;
 
  F[Shortcuts::rhoE] = (Q[Shortcuts::rhoE] + p) * u_n;
""".format(
    compute_primitive_variables=compute_primitive_variables
)

Definition at line 86 of file fv-coupling.py.

◆ euler_initial_conditions

str fv-coupling.euler_initial_conditions

Initial value:

=  """
  Q[Shortcuts::rho]  = 0.1;
  Q[Shortcuts::rhoU] = 0.0;
  Q[Shortcuts::rhoV] = 0.0;
#if DIMENSIONS == 3
  Q[Shortcuts::rhoW] = 0.0;
#endif
  Q[Shortcuts::rhoE] = ((x(0) < 0.5) ? (1.01) : (1.0));
"""

Definition at line 35 of file fv-coupling.py.

◆ euler_max_eigenvalue

str fv-coupling.euler_max_eigenvalue

Initial value:

=  r"""
  {compute_primitive_variables}
  const auto speedOfSound = std::sqrt(GAMMA * p * irho);
  auto result = std::fmax(0.0, std::fabs(u_n - speedOfSound));
  result = std::fmax(result, std::fabs(u_n + speedOfSound));
  return result;
""".format(
    compute_primitive_variables=compute_primitive_variables
)

Definition at line 76 of file fv-coupling.py.

◆ euler_solver

fv-coupling.euler_solver

Initial value:

=  exahype2.solvers.fv.godunov.GlobalAdaptiveTimeStep(
    name="EulerFVSolver",
    patch_size=args.degrees_of_freedom,
    unknowns=euler_unknowns,
    auxiliary_variables=0,
    min_volume_h=min_h,
    max_volume_h=max_h,
    time_step_relaxation=0.5,
)

Definition at line 104 of file fv-coupling.py.

◆ euler_unknowns

dict fv-coupling.euler_unknowns

Initial value:

=  {
    "rho": 1,
    "rhoU": 1,
    "rhoV": 1,
}

EULER #.

Definition at line 26 of file fv-coupling.py.

◆ first_plot_time_stamp

fv-coupling.first_plot_time_stamp

Definition at line 293 of file fv-coupling.py.

◆ flux

fv-coupling.flux

Definition at line 117 of file fv-coupling.py.

◆ initial_conditions

fv-coupling.initial_conditions

Definition at line 115 of file fv-coupling.py.

◆ last_solver

fv-coupling.last_solver

Definition at line 265 of file fv-coupling.py.

◆ make

fv-coupling.make

Definition at line 309 of file fv-coupling.py.

◆ make_clean_first

fv-coupling.make_clean_first

Definition at line 309 of file fv-coupling.py.

◆ max_eigenvalue

fv-coupling.max_eigenvalue

Definition at line 118 of file fv-coupling.py.

◆ max_end_time

fv-coupling.max_end_time

Definition at line 292 of file fv-coupling.py.

◆ max_h

float fv-coupling.max_h = 1.1 / 3.0**args.min_depth

Definition at line 19 of file fv-coupling.py.

◆ min_depth

fv-coupling.min_depth

Definition at line 9 of file fv-coupling.py.

◆ min_end_time

fv-coupling.min_end_time

Definition at line 291 of file fv-coupling.py.

◆ min_h

float fv-coupling.min_h = max_h * 3.0 ** (-args.amr_levels)

Definition at line 20 of file fv-coupling.py.

◆ mode

fv-coupling.mode

Definition at line 304 of file fv-coupling.py.

◆ number_of_snapshots

fv-coupling.number_of_snapshots

Definition at line 12 of file fv-coupling.py.

◆ number_of_variables_advection

tuple fv-coupling.number_of_variables_advection

Initial value:

=  (
    advection_solver.auxiliary_variables + advection_solver.unknowns
)

Definition at line 203 of file fv-coupling.py.

◆ number_of_variables_euler

fv-coupling.number_of_variables_euler = euler_solver.auxiliary_variables + euler_solver.unknowns

Definition at line 202 of file fv-coupling.py.

◆ offset

fv-coupling.offset

Definition at line 290 of file fv-coupling.py.

◆ parser

fv-coupling.parser = exahype2.ArgumentParser("ExaHyPE 2 - Finite Volumes Coupling Testing Script")

Definition at line 6 of file fv-coupling.py.

◆ periodic_BC

fv-coupling.periodic_BC

Definition at line 295 of file fv-coupling.py.

◆ periodic_boundary_conditions_x

fv-coupling.periodic_boundary_conditions_x

Definition at line 13 of file fv-coupling.py.

◆ periodic_boundary_conditions_y

fv-coupling.periodic_boundary_conditions_y

Definition at line 14 of file fv-coupling.py.

◆ periodic_boundary_conditions_z

fv-coupling.periodic_boundary_conditions_z

Definition at line 15 of file fv-coupling.py.

◆ project

fv-coupling.project

Initial value:

=  exahype2.Project(
    namespace=["tests", "exahype2", "fv"],
    project_name=".",
    directory=".",
    executable="ExaHyPE",
)

PROJECT #.

Definition at line 271 of file fv-coupling.py.

◆ size

fv-coupling.size

Definition at line 289 of file fv-coupling.py.

◆ throw_away_data_after_build

fv-coupling.throw_away_data_after_build

Definition at line 309 of file fv-coupling.py.

◆ time_in_between_plots

fv-coupling.time_in_between_plots = 0.0

Definition at line 282 of file fv-coupling.py.

◆ True

fv-coupling.True

Definition at line 309 of file fv-coupling.py.

Variables

Variable Documentation

◆ _compute_time_step_size

◆ advection_auxiliary_variables

◆ advection_boundary_conditions

◆ advection_flux

◆ advection_initial_conditions

◆ advection_max_eigenvalue

◆ advection_solver

◆ advection_unknowns

◆ args

◆ boundary_conditions

◆ compute_primitive_variables

◆ degrees_of_freedom

◆ dimensions

◆ end_time

◆ euler_boundary_conditions

◆ euler_flux

◆ euler_initial_conditions

◆ euler_max_eigenvalue

◆ euler_solver

◆ euler_unknowns

◆ first_plot_time_stamp

◆ flux

◆ initial_conditions

◆ last_solver

◆ make

◆ make_clean_first

◆ max_eigenvalue

◆ max_end_time

◆ max_h

◆ min_depth

◆ min_end_time

◆ min_h

◆ mode

◆ number_of_snapshots

◆ number_of_variables_advection

◆ number_of_variables_euler

◆ offset

◆ parser

◆ periodic_BC

◆ periodic_boundary_conditions_x

◆ periodic_boundary_conditions_y

◆ periodic_boundary_conditions_z

◆ project

◆ size

◆ throw_away_data_after_build

◆ time_in_between_plots

◆ True