Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include <utility>
       7              : 
       8              : #include "Evolution/TypeTraits.hpp"
       9              : #include "PointwiseFunctions/AnalyticData/AnalyticData.hpp"
      10              : #include "PointwiseFunctions/AnalyticSolutions/AnalyticSolution.hpp"
      11              : #include "Utilities/Requires.hpp"
      12              : #include "Utilities/TMPL.hpp"
      13              : 
      14            0 : namespace evolution::Initialization {
      15              : /// Extract initial data either from analytic data or from an analytic
      16              : /// solution at a specified time.
      17              : template <typename SolutionOrData, typename Coordinates, typename... Tags,
      18              :           Requires<is_analytic_solution_v<SolutionOrData>> = nullptr>
      19            1 : decltype(auto) initial_data(const SolutionOrData& solution_or_data,
      20              :                             Coordinates&& coordinates, const double time,
      21              :                             const tmpl::list<Tags...> tags) {
      22              :   return solution_or_data.variables(std::forward<Coordinates>(coordinates),
      23              :                                     time, tags);
      24              : }
      25              : 
      26              : /// \cond
      27              : template <typename SolutionOrData, typename Coordinates, typename... Tags,
      28              :           Requires<is_analytic_data_v<SolutionOrData> or
      29              :                    is_numeric_initial_data_v<SolutionOrData>> = nullptr>
      30              : decltype(auto) initial_data(const SolutionOrData& solution_or_data,
      31              :                             Coordinates&& coordinates, const double /*time*/,
      32              :                             const tmpl::list<Tags...> tags) {
      33              :   return solution_or_data.variables(std::forward<Coordinates>(coordinates),
      34              :                                     tags);
      35              : }
      36              : /// \endcond
      37              : }  // namespace evolution::Initialization