code: def logregularize(self, epsilon=2**-1074):
          
          self.numerator.array[self.denominator.array == 0] = epsilon
          self.denominator.array[self.denominator.array == 0] = 1
          return self
docs: Find bins in the denominator that are 0, and set them to 1,
          while setting the corresponding bin in the numerator to
          float epsilon.  This has the effect of allowing the
          logarithm of the ratio array to be evaluated without error.