Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include <complex>
       7              : #include <cstddef>
       8              : #include <limits>
       9              : #include <memory>
      10              : 
      11              : #include "DataStructures/SpinWeighted.hpp"
      12              : #include "DataStructures/Tensor/TypeAliases.hpp"
      13              : #include "Evolution/Systems/Cce/Initialize/InitializeJ.hpp"
      14              : #include "Parallel/NodeLock.hpp"
      15              : #include "Utilities/Gsl.hpp"
      16              : #include "Utilities/Serialization/CharmPupable.hpp"
      17              : 
      18              : /// \cond
      19              : class DataVector;
      20              : class ComplexDataVector;
      21              : /// \endcond
      22              : 
      23              : namespace Cce::Solutions::LinearizedBondiSachs_detail::InitializeJ {
      24              : // First hypersurface Initialization for the
      25              : // `Cce::Solutions::LinearizedBondiSachs` analytic solution.
      26              : //
      27              : // This initialization procedure should not be used except when the
      28              : // `Cce::Solutions::LinearizedBondiSachs` analytic solution is used,
      29              : // as a consequence, this initial data generator is deliberately not
      30              : // option-creatable; it should only be obtained from the `get_initialize_j`
      31              : // function of `Cce::InitializeJ::LinearizedBondiSachs`.
      32              : //
      33              : // It lives in its own lightweight header (depending only on `InitializeJ.hpp`,
      34              : // not the heavy `LinearizedBondiSachs.hpp`/`SphericalMetricData.hpp`) so that
      35              : // `InitializeJ.hpp` can bottom-include it for the `call_with_dynamic_type`
      36              : // dispatch over `InitializeJ<false>::creatable_classes` without an include
      37              : // cycle.
      38              : struct LinearizedBondiSachs : ::Cce::InitializeJ::InitializeJ<false> {
      39              :   WRAPPED_PUPable_decl_template(LinearizedBondiSachs);  // NOLINT
      40              :   explicit LinearizedBondiSachs(CkMigrateMessage* /*unused*/) {}
      41              : 
      42              :   LinearizedBondiSachs() = default;
      43              : 
      44              :   LinearizedBondiSachs(double start_time, double frequency,
      45              :                        std::complex<double> c_2a, std::complex<double> c_2b,
      46              :                        std::complex<double> c_3a, std::complex<double> c_3b);
      47              : 
      48              :   // Deliberately not option-creatable (see note above); obtained only via
      49              :   // `Cce::Solutions::LinearizedBondiSachs::get_initialize_j`. It is still in
      50              :   // `InitializeJ<false>::creatable_classes` so the dispatch and charm
      51              :   // registration see it; `factory_creatable = false` keeps it out of the
      52              :   // option factory.
      53              :   static constexpr bool factory_creatable = false;
      54              : 
      55              :   std::unique_ptr<InitializeJ> get_clone() const override;
      56              : 
      57              :   void operator()(
      58              :       gsl::not_null<Scalar<SpinWeighted<ComplexDataVector, 2>>*> j,
      59              :       gsl::not_null<tnsr::i<DataVector, 3>*> cartesian_cauchy_coordinates,
      60              :       gsl::not_null<
      61              :           tnsr::i<DataVector, 2, ::Frame::Spherical<::Frame::Inertial>>*>
      62              :           angular_cauchy_coordinates,
      63              :       const Scalar<SpinWeighted<ComplexDataVector, 2>>& boundary_j,
      64              :       const Scalar<SpinWeighted<ComplexDataVector, 2>>& boundary_dr_j,
      65              :       const Scalar<SpinWeighted<ComplexDataVector, 0>>& r,
      66              :       const Scalar<SpinWeighted<ComplexDataVector, 0>>& beta, size_t l_max,
      67              :       size_t number_of_radial_points,
      68              :       gsl::not_null<Parallel::NodeLock*> hdf5_lock) const;
      69              : 
      70              :   void pup(PUP::er& /*p*/) override;
      71              : 
      72              :  private:
      73              :   std::complex<double> c_2a_ = std::numeric_limits<double>::signaling_NaN();
      74              :   std::complex<double> c_2b_ = std::numeric_limits<double>::signaling_NaN();
      75              :   std::complex<double> c_3a_ = std::numeric_limits<double>::signaling_NaN();
      76              :   std::complex<double> c_3b_ = std::numeric_limits<double>::signaling_NaN();
      77              :   double frequency_ = std::numeric_limits<double>::signaling_NaN();
      78              :   double time_ = std::numeric_limits<double>::signaling_NaN();
      79              : };
      80              : }  // namespace Cce::Solutions::LinearizedBondiSachs_detail::InitializeJ