numpy/numpy

@bla1089 Just as a general remark, it's often nicest to operate by default with the logarithms of the numbers of interest, rather than ever working in natural space for intermediate values at all