From a8b69dc06e565427bc96d16c443bb4b3392c4365 Mon Sep 17 00:00:00 2001 From: Erik Schnetter Date: Sun, 12 Jun 2011 22:22:13 -0400 Subject: Rename tensor indices Tensor indices with numbers are not supported any more in Kranc; rename them to use letters instead. --- m/McLachlan_ADM.m | 12 ++++++------ m/McLachlan_ADMConstraints.m | 6 +++--- 2 files changed, 9 insertions(+), 9 deletions(-) (limited to 'm') 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], -- cgit v1.2.3