Issue #2520
new

The current gauge setting in the TOV example match what was used in the ET paper so that one produces the same plot. However that choice is somewhat unusual:

ML_BSSN::harmonicN = 1 # these set: ML_BSSN::harmonicF = 1.0 # dot[alpha] = - alpha A ML_BSSN::evolveA = 1 # dot[A] = dot[trK] ML_BSSN::evolveB = 1 # these set: dot[beta] = dot[dot[beta]] = 0 ML_BSSN::ShiftGammaCoeff = 0.0 # by choosing coeffs to zero the RHS ML_BSSN::AlphaDriver = 0.0 ML_BSSN::BetaDriver = 0.0 ML_BSSN::advectLapse = 0 ML_BSSN::advectShift = 0

and it may make more sense to use settings that are closer to BNS setup so that they can be generalized eg:

ML_BSSN::harmonicN = 1 # 1+log ML_BSSN::harmonicF = 2.0 # 1+log ML_BSSN::ShiftGammaCoeff = 0.75 ML_BSSN::AlphaDriver = 0.0 ML_BSSN::BetaDriver = 1.0 ML_BSSN::LapseAdvectionCoeff = 1.0 ML_BSSN::ShiftAdvectionCoeff = 1.0