-#!/usr/bin/python
+#!/usr/bin/python3
from sympy import *
import itertools
q_dirn_sqparm = q_dirn_orgcoords.replace(s, sh**2).replace(t, th**2)
q_sqparm = q_orgcoords .replace(s, sh**2).replace(t, th**2)
-print '----------------------------------------'
+print('----------------------------------------')
dbg('q_sqparm', 'q_dirn_sqparm')
-print '----------------------------------------'
+print('----------------------------------------')
for v in 'sh','th','la','mu':
dbg('diff(q_sqparm,%s)' % v)
dbg('diff(q_dirn_sqparm,%s)' % v)
-print '----------------------------------------'
+print('----------------------------------------')
gamma = symbols('gamma')
sinc(small),
1 - small*small/factorial(3) - small**4/factorial(5)
)
- print smalls
+ print(smalls)
dbg('result_fordiff');
for v in 'sh','th','la','mu','gamma','kappa':
dbg('diff(result_fordiff,%s)' % v)
d = diff(result_fordiff,eval(v))
- print ccode(d, assign_to='FOO')
+ print(ccode(d, assign_to='FOO'))
(defs, v) = cse(d)
for vn, val in defs:
- print ccode(val, assign_to=vn)
- print ccode(v[0], assign_to='BAR')
+ print(ccode(val, assign_to=vn))
+ print(ccode(v[0], assign_to='BAR'))
#bad = q_orgcoords[0]
#badd = diff(bad, la)