10#include <compile_time_options.h>
18 namespace ShallowWaterInitialStates
27 template <
typename Description,
int dim,
typename Number>
37 const std::string subsection)
40 , hyperbolic_system_(hyperbolic_system)
46 "The (positive) slope of the inclined plane used to "
47 "define the bathymetry");
50 this->add_parameter(
"initial discharge",
52 "The initial (unit) discharge in [m^2 / s]");
58 const auto view = hyperbolic_system_.template view<dim, Number>();
65 const Number n = view.manning_friction_coefficient();
66 const Number exponent = 1. / (2. + 4. / 3.);
68 Number profile = n * n * q_0_ * q_0_ / slope_;
76 initial_precomputed_type
final
79 return {compute_bathymetry(point)};
85 DEAL_II_ALWAYS_INLINE
inline Number
86 compute_bathymetry(
const dealii::Point<dim> &point)
const
88 return -slope_ * point[0];
typename Description::HyperbolicSystem HyperbolicSystem
SlopingFriction(const HyperbolicSystem &hyperbolic_system, const std::string subsection)
typename Description::template HyperbolicSystemView< dim, Number > View
typename View::state_type state_type
auto initial_precomputations(const dealii::Point< dim > &point) -> typename InitialState< Description, dim, Number >::initial_precomputed_type final
state_type compute(const dealii::Point< dim > &, Number) final
T pow(const T x, const T b)
Euler::HyperbolicSystem HyperbolicSystem