Given the generalized harmonic variables in the volume, computes the quantities that will be interpolated onto an apparent horizon. More...

#include <ComputeHorizonVolumeQuantities.hpp>

Public Types
using	allowed_src_tags = tmpl::list< gr::Tags::SpacetimeMetric< 3, Frame::Inertial >, GeneralizedHarmonic::Tags::Pi< 3, Frame::Inertial >, GeneralizedHarmonic::Tags::Phi< 3, Frame::Inertial >, Tags::deriv< GeneralizedHarmonic::Tags::Phi< 3, Frame::Inertial >, tmpl::size_t< 3 >, Frame::Inertial > >

using	required_src_tags = tmpl::list< gr::Tags::SpacetimeMetric< 3, Frame::Inertial >, GeneralizedHarmonic::Tags::Pi< 3, Frame::Inertial >, GeneralizedHarmonic::Tags::Phi< 3, Frame::Inertial > >

Static Public Member Functions
template<typename SrcTagList , typename DestTagList >
static void	apply (const gsl::not_null< Variables< DestTagList > * > target_vars, const Variables< SrcTagList > &src_vars, const Mesh< 3 > &mesh) noexcept
	Single-frame case.

template<typename SrcTagList , typename DestTagList , typename TargetFrame >
static void	apply (const gsl::not_null< Variables< DestTagList > * > target_vars, const Variables< SrcTagList > &src_vars, const Mesh< 3 > &mesh, const Jacobian< DataVector, 3, TargetFrame, Frame::Inertial > &jacobian, const InverseJacobian< DataVector, 3, Frame::Logical, TargetFrame > &inverse_jacobian) noexcept
	Dual-frame case.

Detailed Description

Given the generalized harmonic variables in the volume, computes the quantities that will be interpolated onto an apparent horizon.

This is meant to be the primary compute_vars_to_interpolate for the horizon finder.

SrcTagList and DestTagList have limited flexibility, and their restrictions are static_asserted inside the apply functions. The lack of complete flexibility is intentional, because most computations (e.g. for observers) should be done only on the horizon surface (i.e. after interpolation) as opposed to in the volume; only those computations that require data in the volume (e.g. volume numerical derivatives) should be done here.

For the dual-frame case, takes the Jacobians and does numerical derivatives in order to avoid Hessians.

Public Types

Static Public Member Functions

Detailed Description