Implement optimization pass to automatically vectorize functions

[VMatthijs]

For example, to automate things like this:

for (n in 1:240)
  y[n] ~ normal(mu, sigma);

becomes

y[1:240] ~ normal(mu, sigma);