itensor and ctensor: indexed objects in moving frames in specific coordinate systems and frame basis
Subject: itensor and ctensor: indexed objects in moving frames in specific coordinate systems and frame basis
From: Valery Pipin
Date: Wed, 09 Feb 2011 14:58:38 -0800
09.02.2011 15:52, Yaroslav Alekseev ?????:
> Viktor T. Toth <vttoth at vttoth.com> ?????(?) ? ????? ?????? Tue, 08 Feb
> 2011 22:39:04 +0300:
>
>> I think the only step missing is an evaluation of ifb before
>> converting to a
>> ctensor expression. Try this:
>>
>> load(itensor);
>> load(ctensor);
>> iframe_flag:true;
>> cframe_flag:true;
>> ct_coordsys(polar)$
>> depends([V,ifg,ifr,ifri],ct_coords)$
>> ishow(covdiff(V([],[j]),j))$
>> ishow(ev(%,icc2,ifc2,ifc1,ifb))$
>> eqn:eq=%$
>> ceqn:ic_convert(eqn)$
>> ev(ceqn);
>>
>> A word of warning: Presently, ic_convert has no explicit knowledge
>> about the
>> meaning of the symbols ifg, ifr, and ifri. It just treats them as it
>> would
>> treat any other indexed objects. In the present example, it is not
>> relevant
>> as the symbols cancel; in other cases, however, you may need to take
>> steps
>> to identify the itensor symbols ifr, ifri, and ifg with the
>> corresponding
>> ctensor matrices fri, ifri, and lfg/ufg as appropriate.
>>
>>
>> Viktor
>
> Ok, thank you for response.
> I have to define symbols ifr, ifri, and ifg. So, I've tried this:
>
> load(itensor);
> load(ctensor);
> iframe_flag:true;
> cframe_flag:true;
> ct_coordsys(polar)$
> depends(V,ct_coords)$
> ishow(covdiff(V([],[j]),j))$
> ishow(ev(%,icc2,ifc2,ifc1,ifb))$
> eqn:eq=%$
> ceqn:ic_convert(eqn)$
> ifg:matrix([1,0],[0,1]);
> fri:matrix([1,0],[0,1/r]);
> ifri:matrix([1,0],[0,r]);
> ev(ceqn);
>
> But I've got very strange answer:
>
>
> (%i14) ev(ceqn);
> d d
> (%o14) ---- (V ) + -- (V )
> dphi 2 dr 1
>
> It seems to me that the answer should be something like this (ordinary
> divergence):
>
>
> dV2
> ----
> dphi dV1 V1
> (%o21) ---- + --- + --
> r dr r
>
>
Seems you are mixing the co-variant (or con-travariant) components with
ones for orthonormal
but non-coordinate basis. (I mean you can not get the basis component
via differentiation)
Rgds