PacletInstall["BowenPing`STensor`"]
In[]:=
Needs["BowenPing`STensor`"]
In[]:=
g=CreateTensor"g",{t,r,θ,ϕ},DiagonalMatrix-1-,,,,IsMetric->True;
2M
r
-1
1-
2M
r
2
r
2
r
2
Sin[θ]
In[]:=
Rdown=Head[g["ce"]RiemannTensor[g]["e","abd"]]
Out[]=
STensor
In[]:=
Rup=Head[RiemannTensor[g]["a","efg"]Inverse[g]["be"]Inverse[g]["cf"]Inverse[g]["dg"]]
Out[]=
STensor
In[]:=
Rup["abcd"]Rdown["abcd"]
Out[]=
-
24
2
M
6
r