p_rightvars = p_start.subs( theta, s/r ).subs( r, 1/la )
+def dprint(*args):
+ if options.quiet: return
+ print(*args)
+
def dbg(*args):
+ if options.quiet: return
for vn in args:
print('\n ' + vn + '\n')
pprint(eval(vn))
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('----------------------------------------')
+dprint('----------------------------------------')
dbg('q_sqparm', 'q_dirn_sqparm')
-print('----------------------------------------')
+dprint('----------------------------------------')
for v in 'sh','th','la','mu':
dbg('diff(q_sqparm,%s)' % v)
dbg('diff(q_dirn_sqparm,%s)' % v)
-print('----------------------------------------')
+dprint('----------------------------------------')
gamma = symbols('gamma')