import sys
import sympy.galgebra.GA as GA
import sympy.galgebra.latex_ex as tex
GA.set_main(sys.modules[__name__])
if __name__ == '__main__':
tex.Format()
GA.make_symbols('xbm alpha_1 delta__nugamma_r')
x = alpha_1*xbm/delta__nugamma_r
print 'x =',x
tex.xdvi()
#!/usr/bin/python
#Dirac.py
import sympy.galgebra.GA as GA
import sympy.galgebra.latex_ex as tex
import sys
GA.set_main(sys.modules[__name__])
if __name__ == '__main__':
metric = '1 0 0 0,'+\
'0 -1 0 0,'+\
'0 0 -1 0,'+\
'0 0 0 -1'
vars = GA.make_symbols('t x y z')
GA.MV.setup('gamma_t gamma_x gamma_y gamma_z', metric, True, vars)
parms = GA.make_symbols('m e')
tex.Format()
I = GA.MV(GA.ONE, 'pseudo')
nvars = len(vars)
psi = GA.MV('psi', 'spinor', fct=True)
psi.convert_to_blades()
A = GA.MV('A', 'vector', fct=True)
sig_x = gamma_x * gamma_t
sig_y = gamma_y * gamma_t
sig_z = gamma_z * gamma_t
print '$A$ is 4-vector potential'
print A
