Signed-off-by: Ian Jackson <ijackson@chiark.greenend.org.uk>
#!/usr/bin/python
from sympy import *
+import itertools
import sys
)
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')
+ (defs, v) = cse(d)
+ for vn, val in defs:
+ print ccode(val, assign_to=vn)
+ print ccode(v[0], assign_to='BAR')
#bad = q_orgcoords[0]
#badd = diff(bad, la)