Initializes the constraint damping parameters \(\gamma_1\) and \(\gamma_2\). More...

#include <InitializeConstraintGammas.hpp>

Public Types
using	return_tags
using	argument_tags
using	simple_tags = return_tags
using	compute_tags = tmpl::list<>
using	simple_tags_from_options = tmpl::list<>
using	const_global_cache_tags = tmpl::list<>
using	mutable_global_cache_tags = tmpl::list<>

Static Public Member Functions
static void	apply (const gsl::not_null< Scalar< DataVector > * > gamma1, const gsl::not_null< Scalar< DataVector > * > gamma2, const tnsr::I< DataVector, Dim, Frame::Inertial > &inertial_coords)

Detailed Description

template<size_t Dim>
struct CurvedScalarWave::Worldtube::Initialization::InitializeConstraintDampingGammas< Dim >

Initializes the constraint damping parameters \(\gamma_1\) and \(\gamma_2\).

Details

The constraint damping parameters of the CurvedScalarWave system are initialized according to a Gaussian

\begin{align*} \gamma_1 &= 0 \\ \gamma_2 &= A e^{-(\sigma r)^2 } + c, \end{align*}

where \(r = \sqrt{\delta_{ij} x^i x^j}\) and with \(A = 10\), \(\sigma = 10^{-1}\), \(c = 10^{-4}\). These values were experimentally found to ensure a stable evolution when using the worldtube scheme.

DataBox changes: