initial_state_uniform.h

Go to the documentation of this file.

//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
// Copyright (C) 2022 - 2026 by the ryujin authors
//
 
#pragma once
 
#include <ryujin/base/compile_time_options.h>
#include <ryujin/scalar_conservation/hyperbolic_system.h>
#include <ryujin/typed/initial_state.h>
 
namespace ryujin
{
  namespace ScalarConservation
  {
    struct HyperbolicDescription;
 
    template <int dim, typename Number>
    class Uniform : public InitialState<HyperbolicDescription, dim, Number>
    {
    public:
      using View = HyperbolicSystemView<dim, Number>;
      using state_type = typename View::state_type;
 
      Uniform(const HyperbolicSystem &hyperbolic_system,
              const std::string subsection)
          : InitialState<HyperbolicDescription, dim, Number>("uniform",
                                                             subsection)
          , hyperbolic_system(hyperbolic_system)
      {
        primitive_[0] = 1.0;
        this->add_parameter(
            "primitive state", primitive_, "Initial 1d primitive state");
      }
 
      state_type compute(const dealii::Point<dim> & /*point*/,
                         Number /*t*/) final
      {
        const auto view = hyperbolic_system.view<dim, Number>();
        return view.from_primitive_state(view.expand_state(primitive_));
      }
 
    private:
      const HyperbolicSystem &hyperbolic_system;
 
      state_type primitive_;
    };
  } // namespace ScalarConservation
} // namespace ryujin