dot (outer product) + reduce fusion operation

[rjpower]

This depends on #18; we should have an optimization which fuses dot product with a reduction; this will allow us to write applications like k-means in a natural way:

closest = argmin(dot(pts, clusters) - sum(clusters**2, axis=1))