Replace sympy.collect with pymbolic CoefficientCollector in cost model?

The following discussion from !141 (merged) should be addressed:

• @inducer started a discussion:

  Is there a reason why pymbolic's CoefficientCollector couldn't have done the job?