Peano
Loading...
Searching...
No Matches
gaussian-bell Namespace Reference

Variables

str initial_conditions
 
str analytical_solution
 
 parser = exahype2.ArgumentParser("ExaHyPE 2 - Finite Volumes Testing Script")
 
 min_depth
 
 degrees_of_freedom
 
 number_of_snapshots
 
 periodic_boundary_conditions_x
 
 periodic_boundary_conditions_y
 
 periodic_boundary_conditions_z
 
 args = parser.parse_args()
 
list size = [1.0, 1.0, 1.0]
 
float max_h = 1.1 * min(size) / (3.0**args.min_depth)
 
float min_h = max_h * 3.0 ** (-args.amr_levels)
 
 fv_solver
 
 flux
 
 max_eigenvalue
 
 limiter
 
 error_measurement_implementation
 
 output_file_name
 
 project
 
float time_in_between_plots = 0.0
 
 dimensions
 
 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

◆ analytical_solution

str gaussian-bell.analytical_solution
Initial value:
1= """
2 // The initial condition is transported without deformation in a periodic
3 // domain [-.5, +.5], i.e. the value at a given point is the value at the
4 // position x - v*t but accounting for periodicity.
5 // Therefore we find the point x - v*t, and then shift this by instances
6 // of 1.0 until we find the point within the domain.
7#if DIMENSIONS == 2
8 tarch::la::Vector<DIMENSIONS, double> pos = { (x[0] - t) + (int)(t + .5 - x[0]),
9 (x[1] - t) + (int)(t + .5 - x[1]) };
10#else
11 tarch::la::Vector<DIMENSIONS, double> pos = { (x[0] - t) + (int)(t + .5 - x[0]),
12 (x[1] - t) + (int)(t + .5 - x[1]),
13 (x[2] - t) + (int)(t + .5 - x[2]) };
14#endif
15
16 solution[Shortcuts::rho] = 0.02 * (1.0 + std::exp(-0.5 * tarch::la::norm2(pos) * tarch::la::norm2(pos) / 0.01));
17 solution[Shortcuts::rhoU + 0] = 1.0 * solution[Shortcuts::rho];
18 solution[Shortcuts::rhoU + 1] = 1.0 * solution[Shortcuts::rho];
19#if DIMENSIONS == 3
20 solution[Shortcuts::rhoU + 2] = 1.0 * solution[Shortcuts::rho];
21#endif
22
23 const tarch::la::Vector<DIMENSIONS, double> momenta(1.0 * solution[Shortcuts::rho]);
24
25 // Should lead to a constant p over the domain given the initial conditions
26 solution[Shortcuts::rhoE] = PRESSURE / (GAMMA - 1.0) + 0.5 * (momenta * momenta) / solution[Shortcuts::rho];
27"""

Definition at line 23 of file gaussian-bell.py.

◆ args

gaussian-bell.args = parser.parse_args()

Definition at line 61 of file gaussian-bell.py.

◆ degrees_of_freedom

gaussian-bell.degrees_of_freedom

Definition at line 55 of file gaussian-bell.py.

◆ dimensions

gaussian-bell.dimensions

Definition at line 108 of file gaussian-bell.py.

◆ error_measurement_implementation

gaussian-bell.error_measurement_implementation

Definition at line 88 of file gaussian-bell.py.

◆ first_plot_time_stamp

gaussian-bell.first_plot_time_stamp

Definition at line 113 of file gaussian-bell.py.

◆ flux

gaussian-bell.flux

Definition at line 79 of file gaussian-bell.py.

◆ fv_solver

gaussian-bell.fv_solver
Initial value:
1= exahype2.solvers.fv.musclhancock.GlobalAdaptiveTimeStep(
2 name="MUSCLHancockSolver",
3 patch_size=args.degrees_of_freedom,
4 unknowns={"rho": 1, "rhoU": args.dimensions, "rhoE": 1},
5 auxiliary_variables=0,
6 min_volume_h=min_h,
7 max_volume_h=max_h,
8 time_step_relaxation=0.5,
9)

Definition at line 67 of file gaussian-bell.py.

◆ initial_conditions

gaussian-bell.initial_conditions
Initial value:
1= """
2 Q[Shortcuts::rho] = 0.02 * (1.0 + std::exp(-0.5 * tarch::la::norm2(x) * tarch::la::norm2(x) / 0.01));
3
4 Q[Shortcuts::rhoU + 0] = 1.0 * Q[Shortcuts::rho];
5 Q[Shortcuts::rhoU + 1] = 1.0 * Q[Shortcuts::rho];
6#if DIMENSIONS == 3
7 Q[Shortcuts::rhoU + 2] = 1.0 * Q[Shortcuts::rho];
8#endif
9
10 const tarch::la::Vector<DIMENSIONS, double> momenta(1.0 * Q[Shortcuts::rho]);
11
12 // Should lead to a constant pressure over the domain given the initial conditions
13 Q[Shortcuts::rhoE] = PRESSURE / (GAMMA - 1.0) + 0.5 * (momenta * momenta) / Q[Shortcuts::rho];
14"""

Definition at line 8 of file gaussian-bell.py.

◆ limiter

gaussian-bell.limiter

Definition at line 81 of file gaussian-bell.py.

◆ make

gaussian-bell.make

Definition at line 130 of file gaussian-bell.py.

◆ make_clean_first

gaussian-bell.make_clean_first

Definition at line 130 of file gaussian-bell.py.

◆ max_eigenvalue

gaussian-bell.max_eigenvalue

Definition at line 80 of file gaussian-bell.py.

◆ max_end_time

gaussian-bell.max_end_time

Definition at line 112 of file gaussian-bell.py.

◆ max_h

float gaussian-bell.max_h = 1.1 * min(size) / (3.0**args.min_depth)

Definition at line 64 of file gaussian-bell.py.

◆ min_depth

gaussian-bell.min_depth

Definition at line 53 of file gaussian-bell.py.

◆ min_end_time

gaussian-bell.min_end_time

Definition at line 111 of file gaussian-bell.py.

◆ min_h

float gaussian-bell.min_h = max_h * 3.0 ** (-args.amr_levels)

Definition at line 65 of file gaussian-bell.py.

◆ mode

gaussian-bell.mode

Definition at line 125 of file gaussian-bell.py.

◆ number_of_snapshots

gaussian-bell.number_of_snapshots

Definition at line 56 of file gaussian-bell.py.

◆ offset

gaussian-bell.offset

Definition at line 110 of file gaussian-bell.py.

◆ output_file_name

gaussian-bell.output_file_name

Definition at line 89 of file gaussian-bell.py.

◆ parser

gaussian-bell.parser = exahype2.ArgumentParser("ExaHyPE 2 - Finite Volumes Testing Script")

Definition at line 51 of file gaussian-bell.py.

◆ periodic_BC

gaussian-bell.periodic_BC

Definition at line 115 of file gaussian-bell.py.

◆ periodic_boundary_conditions_x

gaussian-bell.periodic_boundary_conditions_x

Definition at line 57 of file gaussian-bell.py.

◆ periodic_boundary_conditions_y

gaussian-bell.periodic_boundary_conditions_y

Definition at line 58 of file gaussian-bell.py.

◆ periodic_boundary_conditions_z

gaussian-bell.periodic_boundary_conditions_z

Definition at line 59 of file gaussian-bell.py.

◆ project

gaussian-bell.project
Initial value:
1= exahype2.Project(
2 namespace=["tests", "exahype2", "fv"],
3 project_name="GaussianBell",
4 directory=".",
5 executable="ExaHyPE",
6)

Definition at line 92 of file gaussian-bell.py.

◆ size

gaussian-bell.size = [1.0, 1.0, 1.0]

Definition at line 63 of file gaussian-bell.py.

◆ throw_away_data_after_build

gaussian-bell.throw_away_data_after_build

Definition at line 130 of file gaussian-bell.py.

◆ time_in_between_plots

gaussian-bell.time_in_between_plots = 0.0

Definition at line 102 of file gaussian-bell.py.

◆ True

gaussian-bell.True

Definition at line 130 of file gaussian-bell.py.