diff options
author | Peter Diener <diener@amarok.cct.lsu.edu> | 2007-11-09 19:35:45 -0600 |
---|---|---|
committer | Peter Diener <diener@amarok.cct.lsu.edu> | 2007-11-09 19:35:45 -0600 |
commit | 4a30a43b9295d01cc2a3eed40f67adce88d33d61 (patch) | |
tree | 138095f54ea740a4ab44c6f83fd896a58a1ac265 | |
parent | b08330b05518b8083dc368d012010868786e813c (diff) |
Fix RHS for phi
Turns out equation 23 in the shift puncture paper is not a typo. The RHS
phi now converges.
Signed-off-by: Peter Diener <diener@amarok.cct.lsu.edu>
-rw-r--r-- | m/McLachlan.m | 5 |
1 files changed, 3 insertions, 2 deletions
diff --git a/m/McLachlan.m b/m/McLachlan.m index 010d503..d3ccf68 100644 --- a/m/McLachlan.m +++ b/m/McLachlan.m @@ -448,8 +448,9 @@ evolCalcBSSN = R[la,lb] -> Rt[la,lb] + Rphi[la,lb], (* PRD 62, 044034 (2000), eqn. (10) *) - dot[phi] -> - (1/6) alpha trK - + Lie[phi, beta] + (1/6) phi PD[beta[ua],la], + (* PRD 67 084023 (2003), eqn. (16) and (23) *) + dot[phi] -> - (1/6) alpha trK + + Lie[phi, beta] + (1/6) PD[beta[ua],la], (* PRD 62, 044034 (2000), eqn. (9) *) dot[gt[la,lb]] -> - 2 alpha At[la,lb] + Lie[gt[la,lb], beta] - (2/3) gt[la,lb] PD[beta[uc],lc], |