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 <memory>
       8              : #include <string>
       9              : 
      10              : #include "DataStructures/SpinWeighted.hpp"
      11              : #include "DataStructures/Tensor/TypeAliases.hpp"
      12              : #include "Evolution/Systems/Cce/Initialize/InitializeJ.hpp"
      13              : #include "Utilities/Gsl.hpp"
      14              : #include "Utilities/Serialization/CharmPupable.hpp"
      15              : #include "Utilities/TMPL.hpp"
      16              : 
      17              : /// \cond
      18              : class ComplexDataVector;
      19              : /// \endcond
      20              : 
      21              : namespace Cce {
      22              : namespace InitializeJ {
      23              : 
      24              : /*!
      25              :  * \brief Initialize \f$J\f$ on the first hypersurface from provided boundary
      26              :  * values of \f$J\f$, \f$R\f$, and \f$\partial_r J\f$.
      27              :  *
      28              :  * \details This initial data is chosen to take the function:
      29              :  *
      30              :  * \f[ J = \frac{A}{r} + \frac{B}{r^3},\f]
      31              :  *
      32              :  * where
      33              :  *
      34              :  * \f{align*}{
      35              :  * A &= R \left( \frac{3}{2} J|_{r = R} + \frac{1}{2} R \partial_r J|_{r =
      36              :  * R}\right) \notag\\
      37              :  * B &= - \frac{1}{2} R^3 (J|_{r = R} + R \partial_r J|_{r = R})
      38              :  * \f}
      39              :  */
      40              : template <>
      41            1 : struct InverseCubic<true> : InitializeJ<true> {
      42            0 :   using options = tmpl::list<>;
      43            0 :   static constexpr Options::String help = {
      44              :       "Initialization process where J is set to a simple Ansatz with a\n"
      45              :       " A/r + B/r^3 piece such that it is smooth with the Cauchy data at the \n"
      46              :       "worldtube"};
      47              : 
      48            0 :   WRAPPED_PUPable_decl_template(InverseCubic);  // NOLINT
      49            0 :   explicit InverseCubic(CkMigrateMessage* /*unused*/) {}
      50              : 
      51            0 :   InverseCubic() = default;
      52              : 
      53            0 :   std::unique_ptr<InitializeJ> get_clone() const override;
      54              : 
      55            0 :   void operator()(
      56              :       gsl::not_null<Scalar<SpinWeighted<ComplexDataVector, 2>>*> j,
      57              :       gsl::not_null<tnsr::i<DataVector, 3>*> cartesian_cauchy_coordinates,
      58              :       gsl::not_null<
      59              :           tnsr::i<DataVector, 2, ::Frame::Spherical<::Frame::Inertial>>*>
      60              :           angular_cauchy_coordinates,
      61              :       gsl::not_null<tnsr::i<DataVector, 3>*> cartesian_inertial_coordinates,
      62              :       gsl::not_null<
      63              :           tnsr::i<DataVector, 2, ::Frame::Spherical<::Frame::Inertial>>*>
      64              :           angular_inertial_coordinates,
      65              :       const Scalar<SpinWeighted<ComplexDataVector, 2>>& boundary_j,
      66              :       const Scalar<SpinWeighted<ComplexDataVector, 2>>& boundary_dr_j,
      67              :       const Scalar<SpinWeighted<ComplexDataVector, 0>>& r,
      68              :       const Scalar<SpinWeighted<ComplexDataVector, 0>>& beta, size_t l_max,
      69              :       size_t number_of_radial_points,
      70              :       gsl::not_null<Parallel::NodeLock*> hdf5_lock) const;
      71              : 
      72            0 :   void pup(PUP::er& /*p*/) override;
      73              : };
      74              : 
      75              : /*!
      76              :  * \brief Initialize \f$J\f$ on the first hypersurface from provided boundary
      77              :  * values of \f$J\f$, \f$R\f$, and \f$\partial_r J\f$.
      78              :  *
      79              :  * \details This initial data is chosen to take the function:
      80              :  *
      81              :  * \f[ J = \frac{A}{r} + \frac{B}{r^3},\f]
      82              :  *
      83              :  * where
      84              :  *
      85              :  * \f{align*}{
      86              :  * A &= R \left( \frac{3}{2} J|_{r = R} + \frac{1}{2} R \partial_r J|_{r =
      87              :  * R}\right) \notag\\
      88              :  * B &= - \frac{1}{2} R^3 (J|_{r = R} + R \partial_r J|_{r = R})
      89              :  * \f}
      90              :  */
      91              : template <>
      92            1 : struct InverseCubic<false> : InitializeJ<false> {
      93            0 :   using options = tmpl::list<>;
      94            0 :   static constexpr Options::String help = {
      95              :       "Initialization process where J is set to a simple Ansatz with a\n"
      96              :       " A/r + B/r^3 piece such that it is smooth with the Cauchy data at the \n"
      97              :       "worldtube"};
      98              : 
      99            0 :   WRAPPED_PUPable_decl_template(InverseCubic);  // NOLINT
     100            0 :   explicit InverseCubic(CkMigrateMessage* /*unused*/) {}
     101              : 
     102            0 :   InverseCubic() = default;
     103              : 
     104            0 :   std::unique_ptr<InitializeJ> get_clone() const override;
     105              : 
     106            0 :   void operator()(
     107              :       gsl::not_null<Scalar<SpinWeighted<ComplexDataVector, 2>>*> j,
     108              :       gsl::not_null<tnsr::i<DataVector, 3>*> cartesian_cauchy_coordinates,
     109              :       gsl::not_null<
     110              :           tnsr::i<DataVector, 2, ::Frame::Spherical<::Frame::Inertial>>*>
     111              :           angular_cauchy_coordinates,
     112              :       const Scalar<SpinWeighted<ComplexDataVector, 2>>& boundary_j,
     113              :       const Scalar<SpinWeighted<ComplexDataVector, 2>>& boundary_dr_j,
     114              :       const Scalar<SpinWeighted<ComplexDataVector, 0>>& r,
     115              :       const Scalar<SpinWeighted<ComplexDataVector, 0>>& beta, size_t l_max,
     116              :       size_t number_of_radial_points,
     117              :       gsl::not_null<Parallel::NodeLock*> hdf5_lock) const;
     118              : 
     119            0 :   void pup(PUP::er& /*p*/) override;
     120              : };
     121              : }  // namespace InitializeJ
     122              : }  // namespace Cce