numpy.ufunc.reduceat

ufunc.reduceat(a, indices, axis=0, dtype=None, out=None)

Выполняет (локальное) сокращение с указанными срезами по одному осевому.

Для i в range(len(indices)), reduceat вычисляет ufunc.reduce(a[indices[i]:indices[i+1]]), что становится i-ой обобщенной «строкой» параллельно axis в окончательном результате (например, в двумерном массиве, если axis = 0, это становится i-ой строкой, но если axis = 1, это становится i-ым столбцом). Есть три исключения:

• когда i = len(indices) - 1 (так для последнего индекса), indices[i+1] = a.shape[axis].
• если indices[i] >= indices[i + 1], i-ая обобщенная «строка» просто a[indices[i]].
• если indices[i] >= len(a) или indices[i] < 0, возникает ошибка.

Форма выходных данных зависит от размера indices и может быть больше, чем a (это происходит, если len(indices) > a.shape[axis]).

Параметры:

a : array_like Массив, над которым выполняется действие. indices : array_like Связанные индексы, разделенные запятыми (не двоеточием), определяющие срезы для сокращения. axis : int, необязательно Ось, вдоль которой применять reduceat. dtype : код типа данных, необязательно Тип, используемый для представления промежуточных результатов. По умолчанию совпадает с типом данных выходного массива, если он предоставлен, или типом данных входного массива, если выходной массив не предоставлен. out : ndarray, None или кортеж из ndarray и None, необязательно Местоположение, в которое сохраняется результат. Если не предоставлено или None, возвращается только что выделенный массив. Для согласованности с ufunc.__call__, если задано в качестве ключевого аргумента, это может быть заключено в кортеж из одного элемента. Изменено в версии 1.13.0: Разрешены кортежи в качестве ключевого аргумента.

Возвращает:

r : ndarray Значения, полученные после сокращения. Если out был предоставлен, r является ссылкой на out.

Примечания

Дескриптивный пример:

Если a является 1-мерным, функция ufunc.accumulate(a) эквивалентна ufunc.reduceat(a, indices)[::2] где indices является range(len(array) - 1) с нулем, помещенным в каждый другой элемент: indices = zeros(2 * len(a) - 1), indices[1::2] = range(1, len(a)).

Не путайте название этого атрибута: reduceat(a) не обязательно меньше, чем a.

Примеры

Чтобы взять сумму последовательных четырёх значений:

>>> np.add.reduceat(np.arange(8),[0,4, 1,5, 2,6, 3,7])[::2]
array([ 6, 10, 14, 18])

Пример с 2-мерным массивом:

>>> x = np.linspace(0, 15, 16).reshape(4,4)
>>> x
array([[  0.,   1.,   2.,   3.],
       [  4.,   5.,   6.,   7.],
       [  8.,   9.,  10.,  11.],
       [ 12.,  13.,  14.,  15.]])

# reduce such that the result has the following five rows:
# [row1 + row2 + row3]
# [row4]
# [row2]
# [row3]
# [row1 + row2 + row3 + row4]

>>> np.add.reduceat(x, [0, 3, 1, 2, 0])
array([[ 12.,  15.,  18.,  21.],
       [ 12.,  13.,  14.,  15.],
       [  4.,   5.,   6.,   7.],
       [  8.,   9.,  10.,  11.],
       [ 24.,  28.,  32.,  36.]])

# reduce such that result has the following two columns:
# [col1 * col2 * col3, col4]

>>> np.multiply.reduceat(x, [0, 3], 1)
array([[    0.,     3.],
       [  120.,     7.],
       [  720.,    11.],
       [ 2184.,    15.]])