nums.numpy.logaddexp2

nums.numpy.logaddexp2(x1, x2, out=None, where=True, **kwargs)[source]

Logarithm of the sum of exponentiations of the inputs in base-2.

This docstring was copied from numpy.logaddexp2.

Some inconsistencies with the NumS version may exist.

Calculates log2(2**x1 + 2**x2). This function is useful in machine learning when the calculated probabilities of events may be so small as to exceed the range of normal floating point numbers. In such cases the base-2 logarithm of the calculated probability can be used instead. This function allows adding probabilities stored in such a fashion.

Parameters
  • x1 (BlockArray) – Input values. If x1.shape != x2.shape, they must be broadcastable to a common shape (which becomes the shape of the output).

  • x2 (BlockArray) – Input values. If x1.shape != x2.shape, they must be broadcastable to a common shape (which becomes the shape of the output).

  • out (BlockArray, None, or optional) – A location into which the result is stored. If provided, it must have a shape that the inputs broadcast to. If not provided or None, a freshly-allocated array is returned. A tuple (possible only as a keyword argument) must have length equal to the number of outputs.

  • where (BlockArray, optional) – This condition is broadcast over the input. At locations where the condition is True, the out array will be set to the ufunc result. Elsewhere, the out array will retain its original value. Note that if an uninitialized out array is created via the default out=None, locations within it where the condition is False will remain uninitialized.

  • **kwargs – For other keyword-only arguments, see the ufunc docs.

Returns

result – Base-2 logarithm of 2**x1 + 2**x2.

Return type

BlockArray

See also

logaddexp

Logarithm of the sum of exponentiations of the inputs.

Examples

The doctests shown below are copied from NumPy. They won’t show the correct result until you operate get(). >>> prob1 = nps.log(nps.array(1e-50)) # doctest: +SKIP >>> prob2 = nps.log(nps.array(2.5e-50)) # doctest: +SKIP >>> prob12 = nps.logaddexp2(prob1, prob2) # doctest: +SKIP >>> prob1.get(), prob2.get(), prob12.get() # doctest: +SKIP (array(-115.12925465), array(-114.21296392), array(-113.59955523)) >>> 2**prob12 # doctest: +SKIP array(6.35515844e-35)