Histogram back-end

Currently, histograms (in PartialHistogram) are computed with numpy.histogram which is relatively slow.
If it becomes a bottleneck for processing, we might consider implementing a OpenCL/Cuda backend.