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 <pup.h>
       8             : #include <string>
       9             : #include <vector>
      10             : 
      11             : #include "DataStructures/Tensor/TypeAliases.hpp"
      12             : #include "Elliptic/BoundaryConditions/BoundaryCondition.hpp"
      13             : #include "Elliptic/BoundaryConditions/BoundaryConditionType.hpp"
      14             : #include "Options/Context.hpp"
      15             : #include "Options/String.hpp"
      16             : #include "Utilities/Gsl.hpp"
      17             : #include "Utilities/Serialization/CharmPupable.hpp"
      18             : #include "Utilities/TMPL.hpp"
      19             : 
      20             : /// \cond
      21             : class DataVector;
      22             : /// \endcond
      23             : 
      24           0 : namespace sgb::BoundaryConditions {
      25             : 
      26             : /// Do not apply a boundary condition, used exclusively for singular boundary
      27             : /// value problems.
      28             : 
      29           1 : class DoNothing : public elliptic::BoundaryConditions::BoundaryCondition<3> {
      30             :  private:
      31           0 :   using Base = elliptic::BoundaryConditions::BoundaryCondition<3>;
      32             : 
      33             :  public:
      34           0 :   static constexpr Options::String help =
      35             :       "Do not apply a boundary condition, used exclusively for singular "
      36             :       "boundary value problems.";
      37             : 
      38           0 :   using options = tmpl::list<>;
      39             : 
      40           0 :   DoNothing() = default;
      41           0 :   DoNothing(const DoNothing&) = default;
      42           0 :   DoNothing& operator=(const DoNothing&) = default;
      43           0 :   DoNothing(DoNothing&&) = default;
      44           0 :   DoNothing& operator=(DoNothing&&) = default;
      45           0 :   ~DoNothing() override = default;
      46             : 
      47             :   /// \cond
      48             :   explicit DoNothing(CkMigrateMessage* m) : Base(m) {}
      49             :   using PUP::able::register_constructor;
      50             :   WRAPPED_PUPable_decl_template(DoNothing);
      51             :   /// \endcond
      52             : 
      53           0 :   std::unique_ptr<domain::BoundaryConditions::BoundaryCondition> get_clone()
      54             :       const override {
      55             :     return std::make_unique<DoNothing>(*this);
      56             :   }
      57             : 
      58           0 :   std::vector<elliptic::BoundaryConditionType> boundary_condition_types()
      59             :       const override {
      60             :     return {1, elliptic::BoundaryConditionType::Dirichlet};
      61             :   }
      62             : 
      63           0 :   using argument_tags = tmpl::list<>;
      64           0 :   using volume_tags = tmpl::list<>;
      65             : 
      66           0 :   void apply(gsl::not_null<Scalar<DataVector>*> /*field*/,
      67             :              gsl::not_null<Scalar<DataVector>*> /*n_dot_field_gradient*/,
      68             :              const tnsr::i<DataVector, 3>& /*deriv_field*/) const {};
      69             : 
      70           0 :   using argument_tags_linearized = tmpl::list<>;
      71           0 :   using volume_tags_linearized = tmpl::list<>;
      72             : 
      73           0 :   void apply_linearized(
      74             :       gsl::not_null<Scalar<DataVector>*> /*field_correction*/,
      75             :       gsl::not_null<Scalar<DataVector>*>
      76             :       /*n_dot_field_gradient_correction*/,
      77             :       const tnsr::i<DataVector, 3>& /*deriv_field_correction*/) const {};
      78             : };
      79             : 
      80           0 : inline bool operator==(const DoNothing& /*lhs*/, const DoNothing& /*rhs*/) {
      81             :   return true;
      82             : }
      83             : 
      84           0 : inline bool operator!=(const DoNothing& /*lhs*/, const DoNothing& /*rhs*/) {
      85             :   return false;
      86             : }
      87             : 
      88             : /// \cond
      89             : PUP::able::PUP_ID DoNothing::my_PUP_ID = 0;  // NOLINT
      90             : /// \endcond
      91             : 
      92             : }  // namespace sgb::BoundaryConditions