aboutsummaryrefslogtreecommitdiff
path: root/m
diff options
context:
space:
mode:
authorIan Hinder <ian.hinder@aei.mpg.de>2011-06-13 11:22:40 +0200
committerIan Hinder <ian.hinder@aei.mpg.de>2011-06-13 11:22:40 +0200
commit7394810db2087ecd54afcac7cb703b145a83e987 (patch)
treecabc39bebecf54d332040bd16e83498a5c0240cf /m
parenta8b69dc06e565427bc96d16c443bb4b3392c4365 (diff)
Revert "Rename tensor indices"
This reverts commit a8b69dc06e565427bc96d16c443bb4b3392c4365.
Diffstat (limited to 'm')
-rw-r--r--m/McLachlan_ADM.m12
-rw-r--r--m/McLachlan_ADMConstraints.m6
2 files changed, 9 insertions, 9 deletions
diff --git a/m/McLachlan_ADM.m b/m/McLachlan_ADM.m
index ead5522..1bde284 100644
--- a/m/McLachlan_ADM.m
+++ b/m/McLachlan_ADM.m
@@ -236,9 +236,9 @@ evolCalc =
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[ud,lc,lb] - G[uc,la,lb] G[ud,lc,ld]
- + 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]),
+ R[la,lb] -> G[u1,l2,la] G[u2,l1,lb] - G[u1,la,lb] G[u2,l1,l2]
+ + 1/2 gu[u1,u2] (- PD[g[l1,l2],la,lb] + PD[g[l1,la],l2,lb]
+ - PD[g[la,lb],l1,l2] + PD[g[l2,lb],l1,la]),
Km[ua,lb] -> gu[ua,uc] K[lc,lb],
trK -> Km[ua,la],
@@ -288,9 +288,9 @@ 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]
- + 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]),
+ R[la,lb] -> G[u1,l2,la] G[l1,lb,u2] - G[u1,la,lb] G[l1,l2,u2]
+ + 1/2 gu[u1,u2] (- PD[g[l1,l2],la,lb] + PD[g[l1,la],l2,lb]
+ - PD[g[la,lb],l1,l2] + PD[g[l2,lb],l1,la]),
trR -> R[la,lb] gu[ua,ub],
Km[ua,lb] -> gu[ua,uc] K[lc,lb],
trK -> Km[ua,la],
diff --git a/m/McLachlan_ADMConstraints.m b/m/McLachlan_ADMConstraints.m
index 3e3ad3a..2925af2 100644
--- a/m/McLachlan_ADMConstraints.m
+++ b/m/McLachlan_ADMConstraints.m
@@ -158,9 +158,9 @@ ADMConstraintsCalc =
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]
- + 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]),
+ R[la,lb] -> G[u1,l2,la] G[l1,lb,u2] - G[u1,la,lb] G[l1,l2,u2]
+ + 1/2 gu[u1,u2] (- PD[g[l1,l2],la,lb] + PD[g[l1,la],l2,lb]
+ - PD[g[la,lb],l1,l2] + PD[g[l2,lb],l1,la]),
trR -> R[la,lb] gu[ua,ub],
Km[ua,lb] -> gu[ua,uc] K[lc,lb],