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 <string>
       8             : #include <vector>
       9             : 
      10             : #include "DataStructures/Tensor/TypeAliases.hpp"
      11             : #include "Domain/Tags.hpp"
      12             : #include "Elliptic/BoundaryConditions/BoundaryCondition.hpp"
      13             : #include "Elliptic/BoundaryConditions/BoundaryConditionType.hpp"
      14             : #include "Options/String.hpp"
      15             : #include "Utilities/Gsl.hpp"
      16             : #include "Utilities/TMPL.hpp"
      17             : 
      18             : /// \cond
      19             : class DataVector;
      20             : /// \endcond
      21             : 
      22             : namespace Punctures::BoundaryConditions {
      23             : 
      24             : /// Impose asymptotic flatness boundary conditions $\partial_r(ru)=0$
      25           1 : class Flatness : public elliptic::BoundaryConditions::BoundaryCondition<3> {
      26             :  private:
      27           0 :   using Base = elliptic::BoundaryConditions::BoundaryCondition<3>;
      28             : 
      29             :  public:
      30           0 :   static constexpr Options::String help = "Asymptotic flatness d_r(ru)=0";
      31           0 :   using options = tmpl::list<>;
      32             : 
      33           0 :   Flatness() = default;
      34           0 :   Flatness(const Flatness&) = default;
      35           0 :   Flatness& operator=(const Flatness&) = default;
      36           0 :   Flatness(Flatness&&) = default;
      37           0 :   Flatness& operator=(Flatness&&) = default;
      38           0 :   ~Flatness() = default;
      39             : 
      40             :   /// \cond
      41             :   explicit Flatness(CkMigrateMessage* m) : Base(m) {}
      42             :   using PUP::able::register_constructor;
      43             :   WRAPPED_PUPable_decl_template(Flatness);
      44             :   /// \endcond
      45             : 
      46           0 :   std::unique_ptr<domain::BoundaryConditions::BoundaryCondition> get_clone()
      47             :       const override {
      48             :     return std::make_unique<Flatness>(*this);
      49             :   }
      50             : 
      51           0 :   std::vector<elliptic::BoundaryConditionType> boundary_condition_types()
      52             :       const override {
      53             :     return {elliptic::BoundaryConditionType::Neumann};
      54             :   }
      55             : 
      56           0 :   using argument_tags =
      57             :       tmpl::list<domain::Tags::Coordinates<3, Frame::Inertial>>;
      58           0 :   using volume_tags = tmpl::list<>;
      59             : 
      60           0 :   static void apply(gsl::not_null<Scalar<DataVector>*> field,
      61             :                     gsl::not_null<Scalar<DataVector>*> n_dot_field_gradient,
      62             :                     const tnsr::I<DataVector, 3>& x);
      63             : 
      64           0 :   using argument_tags_linearized =
      65             :       tmpl::list<domain::Tags::Coordinates<3, Frame::Inertial>>;
      66           0 :   using volume_tags_linearized = tmpl::list<>;
      67             : 
      68           0 :   static void apply_linearized(
      69             :       gsl::not_null<Scalar<DataVector>*> field_correction,
      70             :       gsl::not_null<Scalar<DataVector>*> n_dot_field_gradient_correction,
      71             :       const tnsr::I<DataVector, 3>& x);
      72             : };
      73             : 
      74           0 : bool operator==(const Flatness& lhs, const Flatness& rhs);
      75           0 : bool operator!=(const Flatness& lhs, const Flatness& rhs);
      76             : 
      77             : }  // namespace Punctures::BoundaryConditions