Add constant folding and propagation patterns

[eb8680]

In #161 and #157 we've seen test failures caused by a failure to reduce expressions like 0 + Gaussian(...) or 1 * Gaussian(...). We should add any missing patterns for constant folding and propagation and move them into the normalize interpreter in #156.

Partly handled in #165 (units only)