aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorIan Hinder <ian.hinder@aei.mpg.de>2011-06-13 11:36:30 +0200
committerIan Hinder <ian.hinder@aei.mpg.de>2011-06-13 11:36:30 +0200
commitec24d427970e8d107a89324bce46238677f49810 (patch)
tree3fd3638a788df18c9a1db83bd4d26e667baac1b4
parent7322f58be5dd767c71438116615eae1e1e3da0d0 (diff)
Rename tensor indices
Tensor indices with numbers are not supported any more in Kranc; rename them to use letters instead. (originally by Erik Schnetter)
-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 1bde284..ead5522 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[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]),
+ 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]),
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[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]),
+ 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]),
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 2925af2..3e3ad3a 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[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]),
+ 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]),
trR -> R[la,lb] gu[ua,ub],
Km[ua,lb] -> gu[ua,uc] K[lc,lb],