Line data    Source code

       1           0 : // Distributed under the MIT License.
       2             : // See LICENSE.txt for details.
       3             : 
       4             : #pragma once
       5             : 
       6             : #include <cstddef>
       7             : #include <optional>
       8             : #include <type_traits>
       9             : 
      10             : #include "DataStructures/DataBox/PrefixHelpers.hpp"
      11             : #include "DataStructures/DataBox/Tag.hpp"
      12             : #include "DataStructures/Tensor/Tensor.hpp"
      13             : #include "DataStructures/Variables.hpp"
      14             : #include "DataStructures/VariablesTag.hpp"
      15             : #include "Domain/Tags.hpp"
      16             : #include "ParallelAlgorithms/Events/Tags.hpp"
      17             : #include "PointwiseFunctions/AnalyticSolutions/AnalyticSolution.hpp"
      18             : #include "PointwiseFunctions/AnalyticSolutions/Tags.hpp"
      19             : #include "PointwiseFunctions/InitialDataUtilities/InitialData.hpp"
      20             : #include "PointwiseFunctions/InitialDataUtilities/Tags/InitialData.hpp"
      21             : #include "Utilities/CallWithDynamicType.hpp"
      22             : #include "Utilities/Gsl.hpp"
      23             : #include "Utilities/TMPL.hpp"
      24             : 
      25             : /// \cond
      26             : namespace Tags {
      27             : struct Time;
      28             : }  // namespace Tags
      29             : /// \endcond
      30             : 
      31             : namespace evolution::Tags {
      32             : /*!
      33             :  * \brief Computes the analytic solution and adds `::Tags::Analytic` of the
      34             :  * `std::optional<Tensor>`s to the DataBox.
      35             :  *
      36             :  * \note If `InitialDataList` is not an empty `tmpl::list`, then
      37             :  * `evolution::initial_data::Tags::InitialData` is retrieved and downcast to the
      38             :  * initial data for computing the errors.
      39             :  */
      40             : template <size_t Dim, typename AnalyticFieldsTagList, bool UsingDgSubcell,
      41             :           typename InitialDataList = tmpl::list<>>
      42           1 : struct AnalyticSolutionsCompute
      43             :     : ::Tags::AnalyticSolutions<AnalyticFieldsTagList>,
      44             :       db::ComputeTag {
      45           0 :   using field_tags = AnalyticFieldsTagList;
      46           0 :   using base = ::Tags::AnalyticSolutions<AnalyticFieldsTagList>;
      47           0 :   using return_type = typename base::type;
      48           0 :   using argument_tags =
      49             :       tmpl::list<evolution::initial_data::Tags::InitialData,
      50             :                  tmpl::conditional_t<
      51             :                      UsingDgSubcell,
      52             :                      Events::Tags::ObserverCoordinates<Dim, Frame::Inertial>,
      53             :                      domain::Tags::Coordinates<Dim, Frame::Inertial>>,
      54             :                  ::Tags::Time>;
      55             : 
      56           0 :   static void function(
      57             :       const gsl::not_null<return_type*> analytic_solution,
      58             :       const evolution::initial_data::InitialData& initial_data,
      59             :       const tnsr::I<DataVector, Dim, Frame::Inertial>& inertial_coords,
      60             :       const double time) {
      61             :     call_with_dynamic_type<void, InitialDataList>(
      62             :         &initial_data, [&analytic_solution, &inertial_coords,
      63             :                         time](const auto* const data_or_solution) {
      64             :           function(analytic_solution, *data_or_solution, inertial_coords, time);
      65             :         });
      66             :   }
      67             : 
      68             :  private:
      69             :   template <typename AnalyticSolution>
      70           0 :   static void function(
      71             :       const gsl::not_null<return_type*> analytic_solution,
      72             :       const AnalyticSolution& analytic_solution_computer,
      73             :       const tnsr::I<DataVector, Dim, Frame::Inertial>& inertial_coords,
      74             :       const double time) {
      75             :     if constexpr (::is_analytic_solution_v<AnalyticSolution>) {
      76             :       *analytic_solution =
      77             :           variables_from_tagged_tuple(analytic_solution_computer.variables(
      78             :               inertial_coords, time, AnalyticFieldsTagList{}));
      79             :     } else {
      80             :       (void)analytic_solution_computer;
      81             :       (void)inertial_coords;
      82             :       (void)time;
      83             :       *analytic_solution = std::nullopt;
      84             :     }
      85             :   }
      86             : };
      87             : }  // namespace evolution::Tags