#1646: McLachlan does not give a reference for the gauge evolution equations ------------------------------------+--------------------------------------- Reporter: rhaas | Owner: Type: defect | Status: confirmed Priority: minor | Milestone: Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: McLachlan ------------------------------------+---------------------------------------
Comment (by rhaas):
No references for the gauge evolution are given in the rewrite code (as of git hash 387e14a) where the relevant Kranc code still reads: {{{
dotalpha -> - (harmonicF alpha^harmonicN IfA[A, + trK - IfCCZ4[2 Theta, 0] + alphaDriver (alpha - 1)]), dot[alpha] -> (+ dotalpha + IfThen[advectLapse!=0, Upwind[beta[ua], alpha, la], 0] + Dissipation[alpha]),
... (* TODO: add not only advection terms, but full Lie derivatives! *) dot[beta[ua]] -> (+ dotbeta[ua] + IfThen[advectShift!=0, Upwind[beta[ub], beta[ua], lb], 0] + Dissipation[beta[ua]]),
dot[B[ua]] -> IfB[+ dotXt[ua] + IfThen[fixAdvectionTerms!=0, Upwind[beta[ub], Xt[ua], lb], 0] - (betaDriverValue (+ B[ua] + IfThen[fixAdvectionTerms!=0 && advectShift!=0, Upwind[beta[ub], beta[ua], lb] / (shiftGammaCoeffValue alpha^shiftAlphaPower), 0])) + IfThen[fixAdvectionTerms==0 && advectShift!=0, Upwind[beta[ub], B[ua], lb], 0] + Dissipation[B[ua]], 0],
}}} ie not reference is given at all.