DownsampleFactorMax to perform max pooling over the trailing two dimensions.

Added tensor.prod function to go with it.

Update the tutorial example on random numbers to account for the default_update mechanism. Update the corresponding test.