[ET Trac] [Einstein Toolkit] #1646: McLachlan does not give a reference for the gauge evolution equations

[Einstein Toolkit]

#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.

-- 
Ticket URL: <https://trac.einsteintoolkit.org/ticket/1646#comment:2>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit