o*******w
发帖数: 349 | 1 我用的是 Maxima 5.17.1 不知是不是个bug. 算出的 Einstein tensor 给出的不是
symmetric. 这不对呀。Ricci tensor 结果是对称的。
有谁帮忙在你的类似软件上run一下。我可以贴我的代码，和结果你 run 一下就行。就
几行。 | m*****n
发帖数: 3575 | | o*******w
发帖数: 349 | 3 如果你方便的话，不知能不能帮我 run 一下 (不才对这两个软件知之甚少) ?
这是metric tensor,
( [1 - t2*C2/(M2*r2), C/M2/r*t*m , 0],
[C/M2/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 - t2*C2/(M2*r2), C/M2/r*t*m , 0],
[C/M2/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)
【在 m*****n 的大作中提到】
: R
: Octave
| m*****n
发帖数: 3575 | 4 你到统计版问问
【在 o*******w 的大作中提到】
: 如果你方便的话，不知能不能帮我 run 一下 (不才对这两个软件知之甚少) ?
: 这是metric tensor,
: ( [1 - t2*C2/(M2*r2), C/M2/r*t*m , 0],
: [C/M2/r*t*m, -(m/M)2, 0],
: [0, 0, -1]
: );
: 我的 Maxima 代码是这个, 您参考：
: load(ctensor);
: (ratwtlvl:false,ratfac:true);
: dim:3;
| o*******w
发帖数: 349 | 5 那的人主要用 sas 和 R 吧。据愚所知这两个都没有算tensor 的 package （或叫
module）
【在 m*****n 的大作中提到】
: 你到统计版问问
| o*******w
发帖数: 349 | 6 我用的是 Maxima 5.17.1 不知是不是个bug. 算出的 Einstein tensor 给出的不是
symmetric. 这不对呀。Ricci tensor 结果是对称的。
有谁帮忙在你的类似软件上run一下。我可以贴我的代码和结果。你 run 一下就行。就
几行。 | m*****n
发帖数: 3575 | | o*******w
发帖数: 349 | 8 如果你方便的话，不知能不能帮我 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)
【在 m*****n 的大作中提到】
: R
: Octave
| m*****n
发帖数: 3575 | 9 你到统计版问问
【在 o*******w 的大作中提到】
: 如果你方便的话，不知能不能帮我 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;
| o*******w
发帖数: 349 | 10 那的人主要用 sas 和 R 吧。据愚所知这两个都没有算tensor 的 package （或叫
module）
【在 m*****n 的大作中提到】
: 你到统计版问问
|
|
