Viktor,
Thanks for your quick reply. Can't give you a Maple code at the
moment, will do it later. Please take a look at the original paper I'm
trying to repeat:
http://articles.adsabs.harvard.edu/cgi-bin/nph-iarticle_query?1977ApJ...214..550P&data_type=PDF_HIGH&whole_paper=YES&type=PRINTER&filetype=.pdf
The link points to a PDF file; dual basis is given by expression 12 at
page 4; resulting coefficients are given by expression 15 at p.4.
Ivan
On Sat, Nov 29, 2008 at 12:20 AM, Viktor T. Toth <vttoth at vttoth.com> wrote:
> 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
>
>