o*******w
发帖数: 349 | 1 我用Maxima （5.17.1） 算出的 Einstein tensor 给出的不是symmetric. 这不对呀。
不知是不是个bug. Ricci tensor 结果是对称的。
有谁能帮忙在你的 R（或类似软件)上 run一下。
这是我的 metric tensor ,
( [1 - t^2*C^2/(M^2*r^2), C/M^2/r*t*m , 0],
[C/M^2/r*t*m, -(m/M)^2, 0],
[0, 0, -1]
);
我在 Maxima 上的代码是这个, 您参考：
load(ctensor);
(ratwtlvl:false,ratfac:true);
dim:3;
ct_coords:[t,x,r];
(ct_coords) [t,x,r];
lg:matrix([1 - t^2*C^2/(M^2*r^2), C/M^2/r*t*m , 0],
[C/M^2/r*t*m, -(m/M)^2, 0],
[0, 0, -1]
);
cmetric();
ug;
("christof() computes the Christoffel-symbols of the 1st and 2nd kind")$
christof(mcs);
ricci(true);
einstein(true);
最后的结果是这，
ein1,1 = - 5C^2*t^2/(4*M^2*r^4)
ein1,2 = C*t(2*C^2t^2 + M^2*r2)/(M^2*m*r^5)
ein1,3 = C^2*t/(2*M^2*r^3)
ein2,1 = C*m*t/(M^2*r^3)
ein2,2 = 7C^2*t^2/(4*M^2*r^4)
ein2,3 = C*m/(2*M^2*r^2)
ein3,2 = C/(2*m*r^2)
ein3,3 = C2*t^2/(4*M^2*r^4) |
|
