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              : 
       8              : #include "Elliptic/BoundaryConditions/BoundaryCondition.hpp"
       9              : #include "Elliptic/Protocols/FirstOrderSystem.hpp"
      10              : #include "Elliptic/Systems/ScalarGaussBonnet/Equations.hpp"
      11              : #include "Elliptic/Systems/Xcts/Tags.hpp"
      12              : #include "Evolution/Systems/CurvedScalarWave/Tags.hpp"
      13              : #include "NumericalAlgorithms/LinearOperators/PartialDerivatives.hpp"
      14              : #include "Utilities/ProtocolHelpers.hpp"
      15              : #include "Utilities/TMPL.hpp"
      16              : 
      17              : namespace sgb {
      18              : 
      19              : /*!
      20              :  * \brief The scalar equation in sGB theories of gravity in the decoupled limit,
      21              :  * assuming quasi-stationarity. See \ref sgb for details on the
      22              :  * explicit equation.
      23              :  */
      24            1 : struct FirstOrderSystem
      25              :     : tt::ConformsTo<elliptic::protocols::FirstOrderSystem> {
      26              :  private:
      27            0 :   using field = ::sgb::Tags::Psi;
      28              : 
      29              :  public:
      30            0 :   static constexpr size_t volume_dim = 3;
      31              : 
      32            0 :   using primal_fields = tmpl::list<field>;
      33            0 :   using primal_fluxes =
      34              :       tmpl::list<::Tags::Flux<field, tmpl::size_t<3>, Frame::Inertial>>;
      35              :   // Note that there are many more background fields required for the elliptic
      36              :   // solve, however these are numerically imported and so are not listed here.
      37            0 :   using background_fields = tmpl::list<
      38              :       Xcts::Tags::InverseConformalMetric<DataVector, 3, Frame::Inertial>,
      39              :       Xcts::Tags::ConformalChristoffelContracted<DataVector, 3,
      40              :                                                  Frame::Inertial>>;
      41            0 :   using inv_metric_tag =
      42              :       Xcts::Tags::InverseConformalMetric<DataVector, 3, Frame::Inertial>;
      43              : 
      44            0 :   using fluxes_computer = Fluxes;
      45            0 :   using sources_computer = Sources;
      46            0 :   using sources_computer_linearized = LinearizedSources;
      47              : 
      48            0 :   using boundary_conditions_base =
      49              :       elliptic::BoundaryConditions::BoundaryCondition<3>;
      50            0 :   using modify_boundary_data = void;
      51              : };
      52              : 
      53              : }  // namespace sgb