def Test_Reciprocal_Frame(): Print_Function() Format() coords = symbols('x y z') (ex,ey,ez,grad) = MV.setup('e_x e_y e_z',metric='[1,1,1]',coords=coords) mfvar = (u,v) = symbols('u v') eu = ex+ey ev = ex-ey (eu_r,ev_r) = ReciprocalFrame([eu,ev]) oprint('\\mbox{Frame}',(eu,ev),'\\mbox{Reciprocal Frame}',(eu_r,ev_r)) print r'%\bm{e}_{u}\cdot\bm{e}^{u} =',(eu|eu_r) print r'%\bm{e}_{u}\cdot\bm{e}^{v} =',eu|ev_r print r'%\bm{e}_{v}\cdot\bm{e}^{u} =',ev|eu_r print r'%\bm{e}_{v}\cdot\bm{e}^{v} =',ev|ev_r eu = ex+ey+ez ev = ex-ey (eu_r,ev_r) = ReciprocalFrame([eu,ev]) oprint('\\mbox{Frame}',(eu,ev),'\\mbox{Reciprocal Frame}',(eu_r,ev_r)) print r'%\bm{e}_{u}\cdot\bm{e}^{u} =',eu|eu_r print r'%\bm{e}_{u}\cdot\bm{e}^{v} =',eu|ev_r print r'%\bm{e}_{v}\cdot\bm{e}^{u} =',ev|eu_r print r'%\bm{e}_{v}\cdot\bm{e}^{v} =',ev|ev_r return
def Test_Reciprocal_Frame(): Print_Function() coords = symbols('x y z') (ex,ey,ez,grad) = MV.setup('ex ey ez',metric='[1,1,1]',coords=coords) mfvar = (u,v) = symbols('u v') eu = ex+ey ev = ex-ey (eu_r,ev_r) = ReciprocalFrame([eu,ev]) oprint('Frame',(eu,ev),'Reciprocal Frame',(eu_r,ev_r)) print 'eu.eu_r =',eu|eu_r print 'eu.ev_r =',eu|ev_r print 'ev.eu_r =',ev|eu_r print 'ev.ev_r =',ev|ev_r eu = ex+ey+ez ev = ex-ey (eu_r,ev_r) = ReciprocalFrame([eu,ev]) oprint('Frame',(eu,ev),'Reciprocal Frame',(eu_r,ev_r)) print 'eu.eu_r =',eu|eu_r print 'eu.ev_r =',eu|ev_r print 'ev.eu_r =',ev|eu_r print 'ev.ev_r =',ev|ev_r return