computing Einstein Tensor

79 views
Skip to first unread message

M E C

unread,
Nov 11, 2024, 10:15:41 AMNov 11
to xAct Tensor Computer Algebra
 Hello, I am trying to calculate the Einstein tensor for the metric.

CTensor{{(-1 + (m r n)/(a^2 Cos[χ]^2 + r^2))/(1 + k Cos[χ] r)^2, 0, 0, -((a m r n Sin[χ]^2)/((1 + k Cos[χ] r)^2 (a^2 Cos[χ]^2 + r^2))), 0}, {0, (a^2 Cos[χ]^2 + r^2)/((1 + k Cos[χ] r)^2 (a^2 + r^2 - m r n)), 0, 0, 0}, {0, 0, (a^2 Cos[χ]^2 + r^2)/(1 + k Cos[χ] r)^2, 0, 0}, {-((a m r n Sin[χ]^2)/((1 + k Cos[χ] r)^2 (a^2 Cos[χ]^2 + r^2))), 0, 0, (Sin[χ]^2 (a^2 + r^2 + (a^2 m r n Sin[χ]^2)/(a^2 Cos[χ]^2 + r^2)))/(1 + k Cos[χ] r)^2, 0}, {0, 0, 0, 0, (Cos[χ]^2 r^2)/(1 + k Cos[χ] r)^2}}, {-B, -B}, 0] 

I have followed the instructions:

<< xAct`xTras`
<< xAct`xTensor`
<< xAct`xCoba`
DefConstantSymbol[{a, m, k, n}]
....
DefManifold[M, 5, {α, β, γ, δ, ι, κ, λ, μ, o, Σ, ρ, σ, τ, Υ, ω, ν}] 
DefChart[B, M, {0, 1, 2, 3, 4}, {t, r, χ, θ, φ}]    
\[Rho]2 = r[]^2 + a^2 Cos[\[Chi][]]^2;
\[CapitalDelta] = r[]^2 - m*r[]^n + a^2;
schw1 = CTensor[....
SetCMetric[schw1, B, SignatureOfMetric -> {3, 1, 0}]
MetricCompute[schw1, B, All]
schw1[-\[Alpha], -\[Beta]]
CD = LC[schw1]
Einstein[CD][\[Alpha], -\[Beta]]

  The system is unable to calculate the Einstein tensor. I also tried using MetricCompute[schw1, B, Einstein, All -> False]. Any recommendations? Thank you very much for reading. Atteched de file
atempt.nb

Leo Stein

unread,
Nov 11, 2024, 12:47:52 PMNov 11
to M E C, xAct Tensor Computer Algebra
Dear Milko,

xCoba can calculate the Einstein tensor with no problem; the slowest parts for me are computing Riemann with all indices down (about 45 sec), raising an index to make the (1,3) form of Riemann (about 1 minute), and tracing to compute Ricci (about 20 sec). All of these times include applying Simplify with some specific assumptions at each stage, to try to keep the size of the expressions small. If you give Simplify more information, and can sometimes find simpler expressions. This is done with the CVSimplify option to MetricCompute.

Let me note also the following changes to your notebook: for SignatureOfMetric, you probably mea {4,1,0}, not {3,1,0}; and for CD, you probably want to use CovDOfMetric[schw1] instead of LC (they end up giving the same result, but LC recomputes whereas CovDOfMetric uses the cached value that was computed by MetricCompute).

There are some more hopefully useful usage patterns in the attached notebook.

Best
Leo

--
You received this message because you are subscribed to the Google Groups "xAct Tensor Computer Algebra" group.
To unsubscribe from this group and stop receiving emails from it, send an email to xact+uns...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/xact/2cf4b61e-4dd8-4df5-8728-98c4353684afn%40googlegroups.com.
atempt.nb

M E C

unread,
Nov 11, 2024, 2:33:12 PMNov 11
to xAct Tensor Computer Algebra
  Thank you very much for the help. I will review it in detail.  

M E C

unread,
Dec 10, 2024, 9:39:20 AMDec 10
to xAct Tensor Computer Algebra
  Hello. I have the following question. The following routine is not able to run. Would there be any suggestions to make it run?  
preguntar.nb

M E C

unread,
Dec 10, 2024, 10:02:52 AMDec 10
to xAct Tensor Computer Algebra
  Apologies for my mistake, in the previous email I attached the wrong file. I am now attaching the correct file. My question was: "I have the following question. The following routine is not able to run. Would there be any suggestions to make it run? " thank you very much for reading 
preguntar 1.nb
Reply all
Reply to author
Forward
0 new messages