Replace sympy.collect with pymbolic CoefficientCollector in cost model Closes #125 See merge request !182