SpECTRE  v2024.12.16
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Modules Pages
Cce::Solutions::LinearizedBondiSachs Struct Reference

Computes the analytic data for a Linearized solution to the Bondi-Sachs equations described in . More...

#include <LinearizedBondiSachs.hpp>

Classes

struct  ExtractionRadius
 
struct  Frequency
 
struct  InitialModes
 

Public Types

using options = implementation defined
 
- Public Types inherited from Cce::Solutions::WorldtubeData
using creatable_classes = implementation defined
 
using tags = implementation defined
 The set of available tags provided by the analytic solution.
 

Public Member Functions

 WRAPPED_PUPable_decl_template (LinearizedBondiSachs)
 
 LinearizedBondiSachs (CkMigrateMessage *msg)
 
 LinearizedBondiSachs (const std::array< std::complex< double >, 2 > &mode_constants, double extraction_radius, double frequency)
 
std::unique_ptr< WorldtubeDataget_clone () const override
 
void pup (PUP::er &p) override
 
std::unique_ptr< Cce::InitializeJ::InitializeJ< false > > get_initialize_j (double start_time) const override
 
- Public Member Functions inherited from Cce::Solutions::SphericalMetricData
 WRAPPED_PUPable_abstract (SphericalMetricData)
 
 SphericalMetricData (CkMigrateMessage *msg)
 
 SphericalMetricData (const double extraction_radius)
 
void jacobian (gsl::not_null< SphericaliCartesianJ * > jacobian, size_t l_max) const
 
void inverse_jacobian (gsl::not_null< CartesianiSphericalJ * > inverse_jacobian, size_t l_max) const
 
void dr_inverse_jacobian (gsl::not_null< CartesianiSphericalJ * > dr_inverse_jacobian, size_t l_max) const
 
void pup (PUP::er &p) override
 
- Public Member Functions inherited from Cce::Solutions::WorldtubeData
 WRAPPED_PUPable_abstract (WorldtubeData)
 
 WorldtubeData (const double extraction_radius)
 
 WorldtubeData (CkMigrateMessage *msg)
 
virtual std::unique_ptr< WorldtubeDataget_clone () const =0
 
template<typename... Tags>
tuples::TaggedTuple< Tags... > variables (const size_t output_l_max, const double time, tmpl::list< Tags... >) const
 Retrieve worldtube data represented by the analytic solution, at boundary angular resolution l_max and time time More...
 
void pup (PUP::er &p) override
 
virtual std::unique_ptr< Cce::InitializeJ::InitializeJ< false > > get_initialize_j (const double) const
 
virtual bool use_noninertial_news () const
 

Static Public Attributes

static constexpr Options::String help
 

Protected Member Functions

void prepare_solution (const size_t, const double) const override
 A no-op as the linearized solution does not have substantial shared computation to prepare before the separate component calculations. More...
 
void linearized_bondi_j (gsl::not_null< SpinWeighted< ComplexDataVector, 2 > * > bondi_j, size_t l_max, double time) const
 Computes the linearized solution for J. More...
 
void linearized_bondi_u (gsl::not_null< SpinWeighted< ComplexDataVector, 1 > * > bondi_u, size_t l_max, double time) const
 Compute the linearized solution for U. More...
 
void linearized_bondi_w (gsl::not_null< SpinWeighted< ComplexDataVector, 0 > * > bondi_w, size_t l_max, double time) const
 Computes the linearized solution for W. More...
 
void linearized_dr_bondi_j (gsl::not_null< SpinWeighted< ComplexDataVector, 2 > * > dr_bondi_j, size_t l_max, double time) const
 Computes the linearized solution for rJ. More...
 
void linearized_dr_bondi_u (gsl::not_null< SpinWeighted< ComplexDataVector, 1 > * > dr_bondi_u, size_t l_max, double time) const
 Compute the linearized solution for rU. More...
 
void linearized_dr_bondi_w (gsl::not_null< SpinWeighted< ComplexDataVector, 0 > * > dr_bondi_w, size_t l_max, double time) const
 Computes the linearized solution for rW. More...
 
void linearized_du_bondi_j (gsl::not_null< SpinWeighted< ComplexDataVector, 2 > * > du_bondi_j, size_t l_max, double time) const
 Computes the linearized solution for uJ. More...
 
void linearized_du_bondi_u (gsl::not_null< SpinWeighted< ComplexDataVector, 1 > * > du_bondi_u, size_t l_max, double time) const
 Compute the linearized solution for uU. More...
 
void linearized_du_bondi_w (gsl::not_null< SpinWeighted< ComplexDataVector, 0 > * > du_bondi_w, size_t l_max, double time) const
 Computes the linearized solution for uW. More...
 
void spherical_metric (gsl::not_null< tnsr::aa< DataVector, 3, ::Frame::Spherical<::Frame::Inertial > > * > spherical_metric, size_t l_max, double time) const override
 Compute the spherical coordinate metric from the linearized Bondi-Sachs system. More...
 
void dr_spherical_metric (gsl::not_null< tnsr::aa< DataVector, 3, ::Frame::Spherical<::Frame::Inertial > > * > dr_spherical_metric, size_t l_max, double time) const override
 Compute the radial derivative of the spherical coordinate metric from the linearized Bondi-Sachs system. More...
 
void dt_spherical_metric (gsl::not_null< tnsr::aa< DataVector, 3, ::Frame::Spherical<::Frame::Inertial > > * > dt_spherical_metric, size_t l_max, double time) const override
 Compute the time derivative of the spherical coordinate metric from the linearized Bondi-Sachs system. More...
 
void variables_impl (gsl::not_null< Scalar< SpinWeighted< ComplexDataVector, -2 > > * > news, size_t l_max, double time, tmpl::type_< Tags::News >) const override
 Determines the News function from the linearized solution parameters. More...
 
virtual void variables_impl (gsl::not_null< tnsr::i< DataVector, 3 > * > cartesian_coordinates, size_t output_l_max, double time, tmpl::type_< Tags::CauchyCartesianCoords >) const
 
virtual void variables_impl (gsl::not_null< tnsr::i< DataVector, 3 > * > dr_cartesian_coordinates, size_t output_l_max, double time, tmpl::type_< Tags::Dr< Tags::CauchyCartesianCoords > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::aa< DataVector, 3 > * > spacetime_metric, size_t output_l_max, double time, tmpl::type_< gr::Tags::SpacetimeMetric< DataVector, 3 > >) const=0
 Computes the Cartesian spacetime metric from the spherical solution provided by the derived classes. More...
 
virtual void variables_impl (gsl::not_null< tnsr::aa< DataVector, 3 > * > dt_spacetime_metric, size_t output_l_max, double time, tmpl::type_<::Tags::dt< gr::Tags::SpacetimeMetric< DataVector, 3 > > >) const=0
 
virtual void variables_impl (gsl::not_null< tnsr::aa< DataVector, 3 > * > pi, size_t output_l_max, double time, tmpl::type_< gh::Tags::Pi< DataVector, 3 > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::iaa< DataVector, 3 > * > d_spacetime_metric, size_t output_l_max, double time, tmpl::type_< gh::Tags::Phi< DataVector, 3 > >) const=0
 Computes the spatial derivatives of the Cartesian spacetime metric from the spherical solution provided by the derived classes. More...
 
virtual void variables_impl (gsl::not_null< tnsr::ii< DataVector, 3 > * > spatial_metric, size_t output_l_max, double time, tmpl::type_< gr::Tags::SpatialMetric< DataVector, 3 > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::ii< DataVector, 3 > * > dt_spatial_metric, size_t output_l_max, double time, tmpl::type_<::Tags::dt< gr::Tags::SpatialMetric< DataVector, 3 > > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::ii< DataVector, 3 > * > dr_spatial_metric, size_t output_l_max, double time, tmpl::type_< Tags::Dr< gr::Tags::SpatialMetric< DataVector, 3 > > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::I< DataVector, 3 > * > shift, size_t output_l_max, double time, tmpl::type_< gr::Tags::Shift< DataVector, 3 > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::I< DataVector, 3 > * > dt_shift, size_t output_l_max, double time, tmpl::type_<::Tags::dt< gr::Tags::Shift< DataVector, 3 > > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::I< DataVector, 3 > * > dr_shift, size_t output_l_max, double time, tmpl::type_< Tags::Dr< gr::Tags::Shift< DataVector, 3 > > >) const
 
virtual void variables_impl (gsl::not_null< Scalar< DataVector > * > lapse, size_t output_l_max, double time, tmpl::type_< gr::Tags::Lapse< DataVector > >) const
 
virtual void variables_impl (gsl::not_null< Scalar< DataVector > * > dt_lapse, size_t output_l_max, double time, tmpl::type_<::Tags::dt< gr::Tags::Lapse< DataVector > > >) const
 
virtual void variables_impl (gsl::not_null< Scalar< DataVector > * > dr_lapse, size_t output_l_max, double time, tmpl::type_< Tags::Dr< gr::Tags::Lapse< DataVector > > >) const
 
virtual void variables_impl (gsl::not_null< Scalar< SpinWeighted< ComplexDataVector, -2 > > * > news, size_t output_l_max, double time, tmpl::type_< Tags::News >) const=0
 
void variables_impl (gsl::not_null< tnsr::aa< DataVector, 3 > * > spacetime_metric, size_t l_max, double time, tmpl::type_< gr::Tags::SpacetimeMetric< DataVector, 3 > >) const override
 Computes the Cartesian spacetime metric from the spherical solution provided by the derived classes. More...
 
void variables_impl (gsl::not_null< tnsr::aa< DataVector, 3 > * > dt_spacetime_metric, size_t l_max, double time, tmpl::type_< ::Tags::dt< gr::Tags::SpacetimeMetric< DataVector, 3 > > >) const override
 Computes the time derivative of the Cartesian spacetime metric from the spherical solution provided by the derived classes. More...
 
void variables_impl (gsl::not_null< tnsr::iaa< DataVector, 3 > * > d_spacetime_metric, size_t l_max, double time, tmpl::type_< gh::Tags::Phi< DataVector, 3 > >) const override
 Computes the spatial derivatives of the Cartesian spacetime metric from the spherical solution provided by the derived classes. More...
 
- Protected Member Functions inherited from Cce::Solutions::SphericalMetricData
void variables_impl (gsl::not_null< tnsr::aa< DataVector, 3 > * > spacetime_metric, size_t l_max, double time, tmpl::type_< gr::Tags::SpacetimeMetric< DataVector, 3 > >) const override
 Computes the Cartesian spacetime metric from the spherical solution provided by the derived classes. More...
 
void variables_impl (gsl::not_null< tnsr::aa< DataVector, 3 > * > dt_spacetime_metric, size_t l_max, double time, tmpl::type_< ::Tags::dt< gr::Tags::SpacetimeMetric< DataVector, 3 > > >) const override
 Computes the time derivative of the Cartesian spacetime metric from the spherical solution provided by the derived classes. More...
 
void variables_impl (gsl::not_null< tnsr::iaa< DataVector, 3 > * > d_spacetime_metric, size_t l_max, double time, tmpl::type_< gh::Tags::Phi< DataVector, 3 > >) const override
 Computes the spatial derivatives of the Cartesian spacetime metric from the spherical solution provided by the derived classes. More...
 
virtual void spherical_metric (gsl::not_null< tnsr::aa< DataVector, 3, ::Frame::Spherical<::Frame::Inertial > > * > spherical_metric, size_t l_max, double time) const =0
 Must be overriden in the derived class; should compute the spacetime metric of the analytic solution in spherical coordinates. More...
 
virtual void dr_spherical_metric (gsl::not_null< tnsr::aa< DataVector, 3, ::Frame::Spherical<::Frame::Inertial > > * > dr_spherical_metric, size_t l_max, double time) const =0
 Must be overriden in the derived class; should compute the first radial derivative of the spacetime metric of the analytic solution in spherical coordinates. More...
 
virtual void dt_spherical_metric (gsl::not_null< tnsr::aa< DataVector, 3, ::Frame::Spherical<::Frame::Inertial > > * > dt_spherical_metric, size_t l_max, double time) const =0
 Must be overriden in the derived class; should compute the first time derivative of the spacetime metric of the analytic solution in spherical coordinates. More...
 
virtual void variables_impl (gsl::not_null< tnsr::i< DataVector, 3 > * > cartesian_coordinates, size_t output_l_max, double time, tmpl::type_< Tags::CauchyCartesianCoords >) const
 
virtual void variables_impl (gsl::not_null< tnsr::i< DataVector, 3 > * > dr_cartesian_coordinates, size_t output_l_max, double time, tmpl::type_< Tags::Dr< Tags::CauchyCartesianCoords > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::aa< DataVector, 3 > * > spacetime_metric, size_t output_l_max, double time, tmpl::type_< gr::Tags::SpacetimeMetric< DataVector, 3 > >) const=0
 
virtual void variables_impl (gsl::not_null< tnsr::aa< DataVector, 3 > * > dt_spacetime_metric, size_t output_l_max, double time, tmpl::type_<::Tags::dt< gr::Tags::SpacetimeMetric< DataVector, 3 > > >) const=0
 
virtual void variables_impl (gsl::not_null< tnsr::aa< DataVector, 3 > * > pi, size_t output_l_max, double time, tmpl::type_< gh::Tags::Pi< DataVector, 3 > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::iaa< DataVector, 3 > * > d_spacetime_metric, size_t output_l_max, double time, tmpl::type_< gh::Tags::Phi< DataVector, 3 > >) const=0
 
virtual void variables_impl (gsl::not_null< tnsr::ii< DataVector, 3 > * > spatial_metric, size_t output_l_max, double time, tmpl::type_< gr::Tags::SpatialMetric< DataVector, 3 > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::ii< DataVector, 3 > * > dt_spatial_metric, size_t output_l_max, double time, tmpl::type_<::Tags::dt< gr::Tags::SpatialMetric< DataVector, 3 > > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::ii< DataVector, 3 > * > dr_spatial_metric, size_t output_l_max, double time, tmpl::type_< Tags::Dr< gr::Tags::SpatialMetric< DataVector, 3 > > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::I< DataVector, 3 > * > shift, size_t output_l_max, double time, tmpl::type_< gr::Tags::Shift< DataVector, 3 > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::I< DataVector, 3 > * > dt_shift, size_t output_l_max, double time, tmpl::type_<::Tags::dt< gr::Tags::Shift< DataVector, 3 > > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::I< DataVector, 3 > * > dr_shift, size_t output_l_max, double time, tmpl::type_< Tags::Dr< gr::Tags::Shift< DataVector, 3 > > >) const
 
virtual void variables_impl (gsl::not_null< Scalar< DataVector > * > lapse, size_t output_l_max, double time, tmpl::type_< gr::Tags::Lapse< DataVector > >) const
 
virtual void variables_impl (gsl::not_null< Scalar< DataVector > * > dt_lapse, size_t output_l_max, double time, tmpl::type_<::Tags::dt< gr::Tags::Lapse< DataVector > > >) const
 
virtual void variables_impl (gsl::not_null< Scalar< DataVector > * > dr_lapse, size_t output_l_max, double time, tmpl::type_< Tags::Dr< gr::Tags::Lapse< DataVector > > >) const
 
virtual void variables_impl (gsl::not_null< Scalar< SpinWeighted< ComplexDataVector, -2 > > * > news, size_t output_l_max, double time, tmpl::type_< Tags::News >) const=0
 
- Protected Member Functions inherited from Cce::Solutions::WorldtubeData
template<typename Tag >
const auto & cache_or_compute (const size_t output_l_max, const double time) const
 
virtual void prepare_solution (size_t output_l_max, double time) const =0
 
virtual void variables_impl (gsl::not_null< tnsr::i< DataVector, 3 > * > cartesian_coordinates, size_t output_l_max, double time, tmpl::type_< Tags::CauchyCartesianCoords >) const
 
virtual void variables_impl (gsl::not_null< tnsr::i< DataVector, 3 > * > dr_cartesian_coordinates, size_t output_l_max, double time, tmpl::type_< Tags::Dr< Tags::CauchyCartesianCoords > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::aa< DataVector, 3 > * > spacetime_metric, size_t output_l_max, double time, tmpl::type_< gr::Tags::SpacetimeMetric< DataVector, 3 > >) const =0
 
virtual void variables_impl (gsl::not_null< tnsr::aa< DataVector, 3 > * > dt_spacetime_metric, size_t output_l_max, double time, tmpl::type_<::Tags::dt< gr::Tags::SpacetimeMetric< DataVector, 3 > > >) const =0
 
virtual void variables_impl (gsl::not_null< tnsr::aa< DataVector, 3 > * > pi, size_t output_l_max, double time, tmpl::type_< gh::Tags::Pi< DataVector, 3 > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::iaa< DataVector, 3 > * > d_spacetime_metric, size_t output_l_max, double time, tmpl::type_< gh::Tags::Phi< DataVector, 3 > >) const =0
 
virtual void variables_impl (gsl::not_null< tnsr::ii< DataVector, 3 > * > spatial_metric, size_t output_l_max, double time, tmpl::type_< gr::Tags::SpatialMetric< DataVector, 3 > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::ii< DataVector, 3 > * > dt_spatial_metric, size_t output_l_max, double time, tmpl::type_<::Tags::dt< gr::Tags::SpatialMetric< DataVector, 3 > > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::ii< DataVector, 3 > * > dr_spatial_metric, size_t output_l_max, double time, tmpl::type_< Tags::Dr< gr::Tags::SpatialMetric< DataVector, 3 > > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::I< DataVector, 3 > * > shift, size_t output_l_max, double time, tmpl::type_< gr::Tags::Shift< DataVector, 3 > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::I< DataVector, 3 > * > dt_shift, size_t output_l_max, double time, tmpl::type_<::Tags::dt< gr::Tags::Shift< DataVector, 3 > > >) const
 
virtual void variables_impl (gsl::not_null< tnsr::I< DataVector, 3 > * > dr_shift, size_t output_l_max, double time, tmpl::type_< Tags::Dr< gr::Tags::Shift< DataVector, 3 > > >) const
 
virtual void variables_impl (gsl::not_null< Scalar< DataVector > * > lapse, size_t output_l_max, double time, tmpl::type_< gr::Tags::Lapse< DataVector > >) const
 
virtual void variables_impl (gsl::not_null< Scalar< DataVector > * > dt_lapse, size_t output_l_max, double time, tmpl::type_<::Tags::dt< gr::Tags::Lapse< DataVector > > >) const
 
virtual void variables_impl (gsl::not_null< Scalar< DataVector > * > dr_lapse, size_t output_l_max, double time, tmpl::type_< Tags::Dr< gr::Tags::Lapse< DataVector > > >) const
 
virtual void variables_impl (gsl::not_null< Scalar< SpinWeighted< ComplexDataVector, -2 > > * > news, size_t output_l_max, double time, tmpl::type_< Tags::News >) const =0
 

Protected Attributes

std::complex< double > c_2a_ = std::numeric_limits<double>::signaling_NaN()
 
std::complex< double > c_3a_ = std::numeric_limits<double>::signaling_NaN()
 
std::complex< double > c_2b_ = std::numeric_limits<double>::signaling_NaN()
 
std::complex< double > c_3b_ = std::numeric_limits<double>::signaling_NaN()
 
double frequency_ = 0.0
 
- Protected Attributes inherited from Cce::Solutions::WorldtubeData
IntermediateCacheTuple intermediate_cache_
 
double extraction_radius_ = std::numeric_limits<double>::quiet_NaN()
 

Additional Inherited Members

- Static Public Member Functions inherited from Cce::Solutions::SphericalMetricData
static void dr_jacobian (gsl::not_null< SphericaliCartesianJ * > dr_jacobian, size_t l_max)
 
- Protected Types inherited from Cce::Solutions::WorldtubeData
using IntermediateCacheTuple = tuples::tagged_tuple_from_typelist< tmpl::transform< tmpl::list< Tags::CauchyCartesianCoords, Tags::Dr< Tags::CauchyCartesianCoords >, gr::Tags::SpacetimeMetric< DataVector, 3 >, gh::Tags::Pi< DataVector, 3 >, gh::Tags::Phi< DataVector, 3 >, gr::Tags::SpatialMetric< DataVector, 3 >, gr::Tags::Shift< DataVector, 3 >, gr::Tags::Lapse< DataVector >, ::Tags::dt< gr::Tags::SpacetimeMetric< DataVector, 3 > >, ::Tags::dt< gr::Tags::SpatialMetric< DataVector, 3 > >, ::Tags::dt< gr::Tags::Shift< DataVector, 3 > >, ::Tags::dt< gr::Tags::Lapse< DataVector > >, Tags::Dr< gr::Tags::SpatialMetric< DataVector, 3 > >, Tags::Dr< gr::Tags::Shift< DataVector, 3 > >, Tags::Dr< gr::Tags::Lapse< DataVector > >, Tags::News >, tmpl::bind< IntermediateCacheTag, tmpl::_1 > > >
 

Detailed Description

Computes the analytic data for a Linearized solution to the Bondi-Sachs equations described in .

Details

The solution represented by this function is generated with only (2,±2) and (3,±3) modes, and is constructed according to the linearized solution documented in Section VI of . For this solution, we impose additional restrictions that the linearized solution be asymptotically flat so that it is compatible with the gauge transformations performed in the SpECTRE regularity-preserving CCE. Using the notation of , we set:

B2=B3=0C2b=3C2a/ν2C3b=3iC3a/ν3

where C2a and C3a may be specified freely and are taken via input option InitialModes.

Member Function Documentation

◆ dr_spherical_metric()

void Cce::Solutions::LinearizedBondiSachs::dr_spherical_metric ( gsl::not_null< tnsr::aa< DataVector, 3, ::Frame::Spherical<::Frame::Inertial > > * >  dr_spherical_metric,
size_t  l_max,
double  time 
) const
overrideprotectedvirtual

Compute the radial derivative of the spherical coordinate metric from the linearized Bondi-Sachs system.

Details

This function dispatches to the individual computations in this class to determine the Bondi-Sachs scalars for the linearized solution. Once the scalars are determined, the radial derivative of the metric is assembled via (note β=0 in this solution)

rgabdxadxb=(W+rrW2rhABUAUBr2(rhAB)UAUB2r2hABUArUB)(dtdr)2(4rhABUB+2r2((rhAB)UB+hABrUB))(dtdr)dxA+(2rhAB+r2rhAB)dxAdxB,

where indices with capital letters refer to angular coordinates and the angular tensors may be written in terms of spin-weighted scalars. Doing so gives the metric components,

rgtt=(W+rrW2rRe(J¯U2+KUU¯)r2rRe(J¯U2+KUU¯))rgtr=rgttrgtθ=2rRe(KU+JU¯)+r2rRe(KU+JU¯)rgtϕ=2rIm(KU+JU¯)+r2rIm(KU+JU¯)rgrr=rgttrgrθ=rgtθrgrϕ=rgtϕrgθθ=2rRe(J+K)+r2Re(rJ+rK)rgθϕ=2rIm(J)+r2Im(rJ)rgϕϕ=2rRe(KJ)+r2Re(rKrJ),

and all other components are zero.

Implements Cce::Solutions::SphericalMetricData.

◆ dt_spherical_metric()

void Cce::Solutions::LinearizedBondiSachs::dt_spherical_metric ( gsl::not_null< tnsr::aa< DataVector, 3, ::Frame::Spherical<::Frame::Inertial > > * >  dt_spherical_metric,
size_t  l_max,
double  time 
) const
overrideprotectedvirtual

Compute the time derivative of the spherical coordinate metric from the linearized Bondi-Sachs system.

Details

This function dispatches to the individual computations in this class to determine the Bondi-Sachs scalars for the linearized solution. Once the scalars are determined, the metric is assembled via (note β=0 in this solution, and note that we take coordinate t=u in converting to the Cartesian coordinates)

tgabdxadxb=(ruWr2uhABUAUB2r2hABUBuUA)(dtdr)22r2(uhABUB+hABuUB)(dtdr)dxA+r2uhABdxAdxB,

where indices with capital letters refer to angular coordinates and the angular tensors may be written in terms of spin-weighted scalars. Doing so gives the metric components,

tgtt=(ruWr2uRe(J¯U2+KUU¯))tgtr=tgtttgtθ=r2uRe(KU+JU¯)tgtϕ=r2uIm(KU+JU¯)tgrr=tgtttgrθ=tgtθtgrϕ=tgtϕtgθθ=r2Re(uJ+uK)tgθϕ=r2Im(uJ)tgϕϕ=r2Re(uKuJ),

and all other components are zero.

Implements Cce::Solutions::SphericalMetricData.

◆ get_clone()

std::unique_ptr< WorldtubeData > Cce::Solutions::LinearizedBondiSachs::get_clone ( ) const
overridevirtual

◆ get_initialize_j()

std::unique_ptr< Cce::InitializeJ::InitializeJ< false > > Cce::Solutions::LinearizedBondiSachs::get_initialize_j ( double  start_time) const
overridevirtual

Reimplemented from Cce::Solutions::WorldtubeData.

◆ linearized_bondi_j()

void Cce::Solutions::LinearizedBondiSachs::linearized_bondi_j ( gsl::not_null< SpinWeighted< ComplexDataVector, 2 > * >  bondi_j,
size_t  l_max,
double  time 
) const
protected

Computes the linearized solution for J.

Details

The linearized solution for J is given by ,

J=12(2Y22+2Y22)Re(J2(r)eiνu)+60(2Y332Y33)Re(J3(r)eiνu),

where

J2(r)=C2a4rC2b12r3,J3(r)=C3a10riνC3b6r3C3b4r4.

◆ linearized_bondi_u()

void Cce::Solutions::LinearizedBondiSachs::linearized_bondi_u ( gsl::not_null< SpinWeighted< ComplexDataVector, 1 > * >  bondi_u,
size_t  l_max,
double  time 
) const
protected

Compute the linearized solution for U.

Details

The linearized solution for U is given by ,

U=3(1Y22+1Y22)Re(U2(r)eiνu)+6(1Y331Y33)Re(U3(r)eiνu),

where

U2(r)=C2a2r2+iνC2b3r3+C2b4r4U3(r)=C3a2r22ν2C3b3r3+5iνC3b4r4+C3br5

◆ linearized_bondi_w()

void Cce::Solutions::LinearizedBondiSachs::linearized_bondi_w ( gsl::not_null< SpinWeighted< ComplexDataVector, 0 > * >  bondi_w,
size_t  l_max,
double  time 
) const
protected

Computes the linearized solution for W.

Details

The linearized solution for W is given by ,

W=12(0Y22+0Y22)Re(W2(r)eiνu)+12(0Y330Y33)Re(W3(r)eiνu),

where

W2(r)=ν2C2br2+iνC2br3+C2b2r4,W3(r)=2iν3C3br24iν2C3br3+5νC3b2r4+3C3br5.

◆ linearized_dr_bondi_j()

void Cce::Solutions::LinearizedBondiSachs::linearized_dr_bondi_j ( gsl::not_null< SpinWeighted< ComplexDataVector, 2 > * >  dr_bondi_j,
size_t  l_max,
double  time 
) const
protected

Computes the linearized solution for rJ.

Details

The linearized solution for rJ is given by ,

rJ=12(2Y22+2Y22)Re(rJ2(r)eiνu)+60(2Y332Y33)Re(rJ3(r)eiνu),

where

rJ2(r)=C2a4r2+C2b4r4,rJ3(r)=C3a10r2+iνC3b2r4+C3br5.

◆ linearized_dr_bondi_u()

void Cce::Solutions::LinearizedBondiSachs::linearized_dr_bondi_u ( gsl::not_null< SpinWeighted< ComplexDataVector, 1 > * >  dr_bondi_u,
size_t  l_max,
double  time 
) const
protected

Compute the linearized solution for rU.

Details

The linearized solution for rU is given by ,

rU=3(1Y22+1Y22)Re(rU2(r)eiνu)+6(1Y331Y33)Re(rU3(r)eiνu),

where

rU2(r)=C2ar3iνC2br4C2br5rU3(r)=C3ar3+2ν2C3br45iνC3br55C3br6

◆ linearized_dr_bondi_w()

void Cce::Solutions::LinearizedBondiSachs::linearized_dr_bondi_w ( gsl::not_null< SpinWeighted< ComplexDataVector, 0 > * >  dr_bondi_w,
size_t  l_max,
double  time 
) const
protected

Computes the linearized solution for rW.

Details

The linearized solution for W is given by ,

rW=12(0Y22+0Y22)Re(rW2(r)eiνu)+12(0Y330Y33)Re(rW3(r)eiνu),

where

rW2(r)=2ν2C2br33iνC2br42C2br5,rW3(r)=4iν3C3br3+12iν2C3br410νC3br515C3br6.

◆ linearized_du_bondi_j()

void Cce::Solutions::LinearizedBondiSachs::linearized_du_bondi_j ( gsl::not_null< SpinWeighted< ComplexDataVector, 2 > * >  du_bondi_j,
size_t  l_max,
double  time 
) const
protected

Computes the linearized solution for uJ.

Details

The linearized solution for uJ is given by ,

uJ=12(2Y22+2Y22)Re(iνJ2(r)eiνu)+60(2Y332Y33)Re(iνJ3(r)eiνu),

where

J2(r)=C2a4rC2b12r3,J3(r)=C3a10riνC3b6r3C3b4r4.

◆ linearized_du_bondi_u()

void Cce::Solutions::LinearizedBondiSachs::linearized_du_bondi_u ( gsl::not_null< SpinWeighted< ComplexDataVector, 1 > * >  du_bondi_u,
size_t  l_max,
double  time 
) const
protected

Compute the linearized solution for uU.

Details

The linearized solution for U is given by ,

uU=3(2Y22+2Y22)Re(iνU2(r)eiνu)+6(2Y332Y33)Re(iνU3(r)eiνu),

where

U2(r)=C2a2r2+iνC2b3r3+C2b4r4U3(r)=C3a2r22ν2C3b3r3+5iνC3b4r4+C3br5

◆ linearized_du_bondi_w()

void Cce::Solutions::LinearizedBondiSachs::linearized_du_bondi_w ( gsl::not_null< SpinWeighted< ComplexDataVector, 0 > * >  du_bondi_w,
size_t  l_max,
double  time 
) const
protected

Computes the linearized solution for uW.

Details

The linearized solution for uW is given by ,

uW=12(1Y22+1Y22)Re(iνW2(r)eiνu)+12(1Y331Y33)Re(iνW3(r)eiνu),

where

W2(r)=ν2C2br2+iνC2br3+C2b2r4,W3(r)=2iν3C3br24iν2C3br3+5νC3b2r4+3C3br5.

◆ prepare_solution()

void Cce::Solutions::LinearizedBondiSachs::prepare_solution ( const  size_t,
const double   
) const
inlineoverrideprotectedvirtual

A no-op as the linearized solution does not have substantial shared computation to prepare before the separate component calculations.

Implements Cce::Solutions::WorldtubeData.

◆ spherical_metric()

void Cce::Solutions::LinearizedBondiSachs::spherical_metric ( gsl::not_null< tnsr::aa< DataVector, 3, ::Frame::Spherical<::Frame::Inertial > > * >  spherical_metric,
size_t  l_max,
double  time 
) const
overrideprotectedvirtual

Compute the spherical coordinate metric from the linearized Bondi-Sachs system.

Details

This function dispatches to the individual computations in this class to determine the Bondi-Sachs scalars for the linearized solution. Once the scalars are determined, the metric is assembled via (note β=0 in this solution)

ds2=((1+rW)r2hABUAUB)(dtdr)22(dtdr)dr2r2hABUB(dtdr)dxA+r2hABdxAdxB,

where indices with capital letters refer to angular coordinates and the angular tensors may be written in terms of spin-weighted scalars. Doing so gives the metric components,

gtt=(1+rWr2Re(J¯U2+KUU¯))gtr=1gttgrr=2+gttgtθ=r2Re(KU+JU¯)gtϕ=r2Im(KU+JU¯)grθ=gtθgrϕ=gtϕgθθ=r2Re(J+K)gθϕ=r2Im(J)gϕϕ=r2Re(KJ),

and all other components are zero.

Implements Cce::Solutions::SphericalMetricData.

◆ variables_impl() [1/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< Scalar< DataVector > * >  dr_lapse,
size_t  output_l_max,
double  time,
tmpl::type_< Tags::Dr< gr::Tags::Lapse< DataVector > > >   
) const
protectedvirtual

Reimplemented from Cce::Solutions::WorldtubeData.

◆ variables_impl() [2/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< Scalar< DataVector > * >  dt_lapse,
size_t  output_l_max,
double  time,
tmpl::type_<::Tags::dt< gr::Tags::Lapse< DataVector > > >   
) const
protectedvirtual

Reimplemented from Cce::Solutions::WorldtubeData.

◆ variables_impl() [3/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< Scalar< DataVector > * >  lapse,
size_t  output_l_max,
double  time,
tmpl::type_< gr::Tags::Lapse< DataVector > >   
) const
protectedvirtual

Reimplemented from Cce::Solutions::WorldtubeData.

◆ variables_impl() [4/20]

void Cce::Solutions::LinearizedBondiSachs::variables_impl ( gsl::not_null< Scalar< SpinWeighted< ComplexDataVector, -2 > > * >  news,
size_t  l_max,
double  time,
tmpl::type_< Tags::News  
) const
overrideprotectedvirtual

Determines the News function from the linearized solution parameters.

Details

The News is determined from the formula given in ,

N=123(2Y22+2Y22)Re(iν3C2beiνu)+115(2Y332Y33)Re(ν4C3beiνu)

Implements Cce::Solutions::WorldtubeData.

◆ variables_impl() [5/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< Scalar< SpinWeighted< ComplexDataVector, -2 > > * >  news,
size_t  output_l_max,
double  time,
tmpl::type_< Tags::News  
) const
protectedvirtual

◆ variables_impl() [6/20]

void Cce::Solutions::SphericalMetricData::variables_impl ( gsl::not_null< tnsr::aa< DataVector, 3 > * >  dt_spacetime_metric,
size_t  l_max,
double  time,
tmpl::type_< ::Tags::dt< gr::Tags::SpacetimeMetric< DataVector, 3 > > >   
) const
overrideprotected

Computes the time derivative of the Cartesian spacetime metric from the spherical solution provided by the derived classes.

Details

The derived classes provide the time derivative of the spherical metric data via the virtual function SphericalMetricData::dt_spherical_metric() at a resolution determined by the l_max argument. This function performs the coordinate transformation using the Jacobian computed from SphericalMetricData::inverse_jacobian().

◆ variables_impl() [7/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< tnsr::aa< DataVector, 3 > * >  dt_spacetime_metric,
size_t  output_l_max,
double  time,
tmpl::type_<::Tags::dt< gr::Tags::SpacetimeMetric< DataVector, 3 > > >   
) const
protectedvirtual

◆ variables_impl() [8/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< tnsr::aa< DataVector, 3 > * >  pi,
size_t  output_l_max,
double  time,
tmpl::type_< gh::Tags::Pi< DataVector, 3 > >   
) const
protectedvirtual

Reimplemented from Cce::Solutions::WorldtubeData.

◆ variables_impl() [9/20]

void Cce::Solutions::SphericalMetricData::variables_impl ( gsl::not_null< tnsr::aa< DataVector, 3 > * >  spacetime_metric,
size_t  l_max,
double  time,
tmpl::type_< gr::Tags::SpacetimeMetric< DataVector, 3 > >   
) const
overrideprotectedvirtual

Computes the Cartesian spacetime metric from the spherical solution provided by the derived classes.

Details

The derived classes provide spherical metric data via the virtual function SphericalMetricData::spherical_metric() at a resolution determined by the l_max argument. This function performs the coordinate transformation using the Jacobian computed from SphericalMetricData::inverse_jacobian().

Reimplemented from Cce::Solutions::SphericalMetricData.

◆ variables_impl() [10/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< tnsr::aa< DataVector, 3 > * >  spacetime_metric,
size_t  l_max,
double  time,
tmpl::type_< gr::Tags::SpacetimeMetric< DataVector, 3 > >   
) const
protectedvirtual

Computes the Cartesian spacetime metric from the spherical solution provided by the derived classes.

Details

The derived classes provide spherical metric data via the virtual function SphericalMetricData::spherical_metric() at a resolution determined by the l_max argument. This function performs the coordinate transformation using the Jacobian computed from SphericalMetricData::inverse_jacobian().

Reimplemented from Cce::Solutions::SphericalMetricData.

◆ variables_impl() [11/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< tnsr::i< DataVector, 3 > * >  cartesian_coordinates,
size_t  output_l_max,
double  time,
tmpl::type_< Tags::CauchyCartesianCoords  
) const
protectedvirtual

Reimplemented from Cce::Solutions::WorldtubeData.

◆ variables_impl() [12/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< tnsr::i< DataVector, 3 > * >  dr_cartesian_coordinates,
size_t  output_l_max,
double  time,
tmpl::type_< Tags::Dr< Tags::CauchyCartesianCoords > >   
) const
protectedvirtual

Reimplemented from Cce::Solutions::WorldtubeData.

◆ variables_impl() [13/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< tnsr::I< DataVector, 3 > * >  dr_shift,
size_t  output_l_max,
double  time,
tmpl::type_< Tags::Dr< gr::Tags::Shift< DataVector, 3 > > >   
) const
protectedvirtual

Reimplemented from Cce::Solutions::WorldtubeData.

◆ variables_impl() [14/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< tnsr::I< DataVector, 3 > * >  dt_shift,
size_t  output_l_max,
double  time,
tmpl::type_<::Tags::dt< gr::Tags::Shift< DataVector, 3 > > >   
) const
protectedvirtual

Reimplemented from Cce::Solutions::WorldtubeData.

◆ variables_impl() [15/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< tnsr::I< DataVector, 3 > * >  shift,
size_t  output_l_max,
double  time,
tmpl::type_< gr::Tags::Shift< DataVector, 3 > >   
) const
protectedvirtual

Reimplemented from Cce::Solutions::WorldtubeData.

◆ variables_impl() [16/20]

void Cce::Solutions::SphericalMetricData::variables_impl ( gsl::not_null< tnsr::iaa< DataVector, 3 > * >  d_spacetime_metric,
size_t  l_max,
double  time,
tmpl::type_< gh::Tags::Phi< DataVector, 3 > >   
) const
overrideprotectedvirtual

Computes the spatial derivatives of the Cartesian spacetime metric from the spherical solution provided by the derived classes.

Details

The derived classes provide the radial derivative of the spherical metric data via the virtual function SphericalMetricData::dr_spherical_metric() at a resolution determined by the l_max_ argument. This function performs the additional angular derivatives necessary to assemble the full spatial derivative and performs the coordinate transformation to Cartesian coordinates via the Jacobians computed in SphericalMetricData::inverse_jacobian() and SphericalMetricData::inverse_jacobian().

Reimplemented from Cce::Solutions::SphericalMetricData.

◆ variables_impl() [17/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< tnsr::iaa< DataVector, 3 > * >  d_spacetime_metric,
size_t  l_max,
double  time,
tmpl::type_< gh::Tags::Phi< DataVector, 3 > >   
) const
protectedvirtual

Computes the spatial derivatives of the Cartesian spacetime metric from the spherical solution provided by the derived classes.

Details

The derived classes provide the radial derivative of the spherical metric data via the virtual function SphericalMetricData::dr_spherical_metric() at a resolution determined by the l_max_ argument. This function performs the additional angular derivatives necessary to assemble the full spatial derivative and performs the coordinate transformation to Cartesian coordinates via the Jacobians computed in SphericalMetricData::inverse_jacobian() and SphericalMetricData::inverse_jacobian().

Reimplemented from Cce::Solutions::SphericalMetricData.

◆ variables_impl() [18/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< tnsr::ii< DataVector, 3 > * >  dr_spatial_metric,
size_t  output_l_max,
double  time,
tmpl::type_< Tags::Dr< gr::Tags::SpatialMetric< DataVector, 3 > > >   
) const
protectedvirtual

Reimplemented from Cce::Solutions::WorldtubeData.

◆ variables_impl() [19/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< tnsr::ii< DataVector, 3 > * >  dt_spatial_metric,
size_t  output_l_max,
double  time,
tmpl::type_<::Tags::dt< gr::Tags::SpatialMetric< DataVector, 3 > > >   
) const
protectedvirtual

Reimplemented from Cce::Solutions::WorldtubeData.

◆ variables_impl() [20/20]

virtual void Cce::Solutions::WorldtubeData::variables_impl ( gsl::not_null< tnsr::ii< DataVector, 3 > * >  spatial_metric,
size_t  output_l_max,
double  time,
tmpl::type_< gr::Tags::SpatialMetric< DataVector, 3 > >   
) const
protectedvirtual

Reimplemented from Cce::Solutions::WorldtubeData.

Member Data Documentation

◆ help

constexpr Options::String Cce::Solutions::LinearizedBondiSachs::help
staticconstexpr
Initial value:
{
"A linearized Bondi-Sachs analytic solution"}

The documentation for this struct was generated from the following file: