Dear Ivan,
The frame-base code is one of the more recent additions to ctensor, so it is
quite possible that it has bugs.
Can you send me your Maple code (or your hand calculations) so that I can do
a comparison?
Viktor
-----Original Message-----
From: maxima-bounces at math.utexas.edu [mailto:maxima-bounces at math.utexas.edu]
On Behalf Of Ivan Zolotukhin
Sent: Friday, November 28, 2008 6:30 AM
To: maxima at math.utexas.edu
Subject: ctensor: compute Ricci rotation coefficients
Hello,
I'm trying to check with Maxima some paper that calculates Ricci
rotation coefficients in a custom base. Result by Maxima is wrong
since I'm able to repeat calculations made in paper in Maple. Could
you please check my Maxima session log below and tell me what I'm
doing wrong? Sorry if there are blunders, I'm new to Maxima.
load(ctensor);
cframe_flag:true;
dim:3;
ct_coords:[r,p,z];
depends(D, [r, p, z]);
depends(g, [r]);
depends(b, [r]);
depends(W, [r,p]);
lfg:ident(3);
fri:matrix( [ D, -D * ( diff(g, r) * cos(b) - z/r * diff(W, p) ), -D
* r * W ], [0, 1/r, 0], [0, 0, 1]);
cmetric();
trrc(all);
There should be only 4 nonvanishing components, but Maxima produces
more and they all are incorrect, e.g. lcs[2,1,2] = 1 / rD, whereas it
should be 1/r.
_______________________________________________
Maxima mailing list
Maxima at math.utexas.edu
http://www.math.utexas.edu/mailman/listinfo/maxima