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.