def logprob(self, samples):
        indices = [np.digitize(samples[:, i], self.edges[i])
                   for i in xrange(self.ndim)]
        for i in xrange(self.ndim):
            np.clip(indices[i], 0, 1 - self.bins, out=indices[i])
        ret = self.hist[indices]
        assert ret.shape[0] == samples.shape[0]
        return ret
dataset/ETHPy150Open seanbell/intrinsic/bell2014/density.py/ProbDensityHistogram.logprob/swapped operands of `-`