Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include <memory>
       7              : #include <optional>
       8              : #include <string>
       9              : 
      10              : #include "DataStructures/DataVector.hpp"
      11              : #include "DataStructures/Tensor/Tensor.hpp"
      12              : #include "Evolution/BoundaryConditions/Type.hpp"
      13              : #include "Evolution/Systems/CurvedScalarWave/Tags.hpp"
      14              : #include "Evolution/Systems/ScalarTensor/BoundaryConditions/BoundaryCondition.hpp"
      15              : #include "Options/String.hpp"
      16              : #include "PointwiseFunctions/GeneralRelativity/GeneralizedHarmonic/ConstraintDampingTags.hpp"
      17              : #include "PointwiseFunctions/InitialDataUtilities/InitialData.hpp"
      18              : #include "Utilities/Gsl.hpp"
      19              : #include "Utilities/Serialization/CharmPupable.hpp"
      20              : #include "Utilities/TMPL.hpp"
      21              : 
      22              : /// \cond
      23              : namespace Tags {
      24              : struct Time;
      25              : }  // namespace Tags
      26              : namespace domain::Tags {
      27              : template <size_t Dim, typename Frame>
      28              : struct Coordinates;
      29              : }  // namespace domain::Tags
      30              : /// \endcond
      31              : 
      32              : namespace ScalarTensor::BoundaryConditions {
      33              : /*!
      34              :  * \brief Sets Dirichlet boundary conditions using the analytic solution or
      35              :  * analytic data.
      36              :  */
      37            1 : class DirichletAnalytic final : public BoundaryCondition {
      38              :  public:
      39              :   /// \brief What analytic solution/data to prescribe.
      40            1 :   struct AnalyticPrescription {
      41            0 :     static constexpr Options::String help =
      42              :         "What analytic solution/data to prescribe.";
      43            0 :     using type = std::unique_ptr<evolution::initial_data::InitialData>;
      44              :   };
      45            0 :   struct Amplitude {
      46            0 :     using type = double;
      47            0 :     static constexpr Options::String help = {
      48              :         "Amplitude of the scalar at the boundary"};
      49              :   };
      50              : 
      51            0 :   using options = tmpl::list<AnalyticPrescription, Amplitude>;
      52              : 
      53            0 :   static constexpr Options::String help{
      54              :       "DirichletAnalytic boundary conditions."};
      55              : 
      56            0 :   DirichletAnalytic() = default;
      57            0 :   DirichletAnalytic(DirichletAnalytic&&) = default;
      58            0 :   DirichletAnalytic& operator=(DirichletAnalytic&&) = default;
      59            0 :   DirichletAnalytic(const DirichletAnalytic&);
      60            0 :   DirichletAnalytic& operator=(const DirichletAnalytic&);
      61            0 :   ~DirichletAnalytic() override = default;
      62              : 
      63            0 :   DirichletAnalytic(std::unique_ptr<evolution::initial_data::InitialData>
      64              :                         analytic_prescription,
      65              :                     double amplitude);
      66              : 
      67            0 :   explicit DirichletAnalytic(CkMigrateMessage* msg);
      68              : 
      69            0 :   WRAPPED_PUPable_decl_base_template(
      70              :       domain::BoundaryConditions::BoundaryCondition, DirichletAnalytic);
      71              : 
      72            0 :   auto get_clone() const -> std::unique_ptr<
      73              :       domain::BoundaryConditions::BoundaryCondition> override;
      74              : 
      75            0 :   static constexpr evolution::BoundaryConditions::Type bc_type =
      76              :       evolution::BoundaryConditions::Type::Ghost;
      77              : 
      78            0 :   void pup(PUP::er& p) override;
      79              : 
      80            0 :   using dg_interior_evolved_variables_tags = tmpl::list<>;
      81            0 :   using dg_interior_temporary_tags =
      82              :       tmpl::list<domain::Tags::Coordinates<3, Frame::Inertial>,
      83              :                  ::gh::Tags::ConstraintGamma1, ::gh::Tags::ConstraintGamma2,
      84              :                  ::CurvedScalarWave::Tags::ConstraintGamma1,
      85              :                  ::CurvedScalarWave::Tags::ConstraintGamma2>;
      86            0 :   using dg_gridless_tags = tmpl::list<::Tags::Time>;
      87              : 
      88            0 :   std::optional<std::string> dg_ghost(
      89              :       // GH evolved variables
      90              :       gsl::not_null<tnsr::aa<DataVector, 3, Frame::Inertial>*> spacetime_metric,
      91              :       gsl::not_null<tnsr::aa<DataVector, 3, Frame::Inertial>*> pi,
      92              :       gsl::not_null<tnsr::iaa<DataVector, 3, Frame::Inertial>*> phi,
      93              :       // Scalar evolved variables
      94              :       gsl::not_null<Scalar<DataVector>*> psi_scalar,
      95              :       gsl::not_null<Scalar<DataVector>*> pi_scalar,
      96              :       gsl::not_null<tnsr::i<DataVector, 3, Frame::Inertial>*> phi_scalar,
      97              :       // GH temporary variables
      98              :       gsl::not_null<Scalar<DataVector>*> gamma1,
      99              :       gsl::not_null<Scalar<DataVector>*> gamma2,
     100              :       gsl::not_null<Scalar<DataVector>*> lapse,
     101              :       gsl::not_null<tnsr::I<DataVector, 3, Frame::Inertial>*> shift,
     102              :       // Scalar temporary variables
     103              :       gsl::not_null<Scalar<DataVector>*> gamma1_scalar,
     104              :       gsl::not_null<Scalar<DataVector>*> gamma2_scalar,
     105              :       // Inverse metric
     106              :       gsl::not_null<tnsr::II<DataVector, 3, Frame::Inertial>*>
     107              :           inv_spatial_metric,
     108              :       // Mesh variables
     109              :       const std::optional<tnsr::I<DataVector, 3, Frame::Inertial>>&
     110              :           face_mesh_velocity,
     111              :       const tnsr::i<DataVector, 3, Frame::Inertial>& normal_covector,
     112              :       const tnsr::I<DataVector, 3, Frame::Inertial>& normal_vector,
     113              :       // GH interior variables
     114              :       const tnsr::I<DataVector, 3, Frame::Inertial>& coords,
     115              :       const Scalar<DataVector>& interior_gamma1,
     116              :       const Scalar<DataVector>& interior_gamma2,
     117              :       // Scalar interior variables
     118              :       const Scalar<DataVector>& gamma1_interior_scalar,
     119              :       const Scalar<DataVector>& gamma2_interior_scalar, double time) const;
     120              : 
     121              :  private:
     122            0 :   std::unique_ptr<evolution::initial_data::InitialData> analytic_prescription_;
     123            0 :   double amplitude_ = std::numeric_limits<double>::signaling_NaN();
     124              : };
     125              : }  // namespace ScalarTensor::BoundaryConditions