aboutsummaryrefslogtreecommitdiff
path: root/m
diff options
context:
space:
mode:
authorIan Hinder <ian.hinder@aei.mpg.de>2011-11-20 15:44:47 +0100
committerIan Hinder <ian.hinder@aei.mpg.de>2011-11-20 15:44:47 +0100
commit4baaafa6b0c2e5e05398bbefcf872374e03277dc (patch)
treea9f1756005666a5c33056f0882b99821bced9eff /m
parentb3a25f307c9bb55370380b97e6074df1a30f645d (diff)
McLachlan_ADM.m: Fix index errors in constraint computation
Diffstat (limited to 'm')
-rw-r--r--m/McLachlan_ADM.m2
1 files changed, 1 insertions, 1 deletions
diff --git a/m/McLachlan_ADM.m b/m/McLachlan_ADM.m
index 1f778b5..92f0590 100644
--- a/m/McLachlan_ADM.m
+++ b/m/McLachlan_ADM.m
@@ -271,7 +271,7 @@ constraintsCalc =
gu[ua,ub] -> 1/detg detgExpr MatrixInverse[g[ua,ub]],
G[ua,lb,lc] -> 1/2 gu[ua,ud]
(PD[g[lb,ld],lc] + PD[g[lc,ld],lb] - PD[g[lb,lc],ld]),
- R[la,lb] -> G[uc,ld,la] G[lc,lb,ud] - G[uc,la,lb] G[lc,ld,ud]
+ R[la,lb] -> gu[us,ur](G[um,la,lr] G[uk,ls,lb] g[lk,lm] - G[um,la,lb] G[uk,ls,lr] g[lk,lm])
+ 1/2 gu[uc,ud] (- PD[g[lc,ld],la,lb] + PD[g[lc,la],ld,lb]
- PD[g[la,lb],lc,ld] + PD[g[ld,lb],lc,la]),
trR -> R[la,lb] gu[ua,ub],