initial_state_contrast.h

Go to the documentation of this file.

//
// SPDX-License-Identifier: Apache-2.0
// [LANL Copyright Statement]
// Copyright (C) 2022 - 2026 by the ryujin authors
// Copyright (C) 2023 - 2024 by Triad National Security, LLC
//
 
#pragma once
 
#include <ryujin/base/compile_time_options.h>
#include <ryujin/typed/initial_state.h>
 
namespace ryujin
{
  namespace ShallowWaterInitialStates
  {
    template <typename Description, int dim, typename Number>
    class Contrast : public InitialState<Description, dim, Number>
    {
    public:
      using HyperbolicSystem = typename Description::HyperbolicSystem;
      using View = typename HyperbolicSystem::template View<dim, Number>;
      using state_type = typename View::state_type;
 
      Contrast(const HyperbolicSystem &hyperbolic_system,
               const std::string subsection)
          : InitialState<Description, dim, Number>("contrast", subsection)
          , hyperbolic_system_(hyperbolic_system)
      {
        primitive_left_[0] = 1.;
        primitive_left_[1] = 0.0;
        this->add_parameter("primitive state left",
                            primitive_left_,
                            "Initial 1d primitive state (h, u) on the left");
 
        primitive_right_[0] = 1.;
        primitive_right_[1] = 0.0;
        this->add_parameter("primitive state right",
                            primitive_right_,
                            "Initial 1d primitive state (h, u) on the right");
      }
 
      state_type compute(const dealii::Point<dim> &point, Number /*t*/) final
      {
        const auto view = hyperbolic_system_.template view<dim, Number>();
        const auto temp = point[0] > 0. ? primitive_right_ : primitive_left_;
        return view.from_initial_state(temp);
      }
 
    private:
      const HyperbolicSystem &hyperbolic_system_;
 
      dealii::Tensor<1, 2, Number> primitive_left_;
      dealii::Tensor<1, 2, Number> primitive_right_;
    };
  } // namespace ShallowWaterInitialStates
} // namespace ryujin