Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include "DataStructures/DataBox/Protocols/Mutator.hpp"
       7              : #include "Evolution/Systems/Ccz4/FiniteDifference/System.hpp"
       8              : #include "Evolution/Systems/Ccz4/FiniteDifference/Tags.hpp"
       9              : #include "Evolution/Systems/Ccz4/Tags.hpp"
      10              : #include "Utilities/Gsl.hpp"
      11              : #include "Utilities/ProtocolHelpers.hpp"
      12              : #include "Utilities/TMPL.hpp"
      13              : 
      14              : namespace Ccz4::fd {
      15              : /*!
      16              :  * \brief Mutator to enforce constrained evolution after every time step
      17              :  *
      18              :  * Togglable option in input file via option tag defined in Ccz4::fd::Tags.
      19              :  * If constrained_evolution is true, enforce
      20              :  * the unit determinant constraint on the conformal spatial metric
      21              :  * and traceless constraint on the ATilde
      22              :  */
      23            1 : struct EnforceConstrainedEvolution : tt::ConformsTo<db::protocols::Mutator> {
      24            0 :   static constexpr size_t dim = System::volume_dim;
      25            0 :   using return_tags = tmpl::list<::Ccz4::Tags::ConformalMetric<DataVector, dim>,
      26              :                                  ::Ccz4::Tags::ATilde<DataVector, dim>>;
      27            0 :   using argument_tags = tmpl::list<::Ccz4::fd::Tags::ConstrainedEvolution>;
      28              : 
      29            0 :   static void apply(
      30              :       gsl::not_null<tnsr::ii<DataVector, dim>*> conformal_spatial_metric,
      31              :       gsl::not_null<tnsr::ii<DataVector, dim>*> a_tilde,
      32              :       bool constrained_evolution);
      33              : };
      34              : }  // namespace Ccz4::fd