Numero massimo di valori di array numpy

Question

Ho bisogno di un modo rapido per mantenere un massimo di un array Numpy. Ad esempio, se il mio array è stato:

x = numpy.array([11,12,13,20,19,18,17,18,23,21]) 

che vorrei:

numpy.array([11,12,13,20,20,20,20,20,23,23]) 

Ovviamente avrei potuto farlo con un po 'di ciclo:

def running_max(x): 
    result = [x[0]] 
    for val in x: 
     if val > result[-1]: 
      result.append(val) 
     else: 
      result.append(result[-1]) 
    return result 

Ma le mie matrici avere centinaia di migliaia di voci e ho bisogno di chiamare questo molte volte. Sembra che debba esserci un trucco per liberare il loop, ma non riesco a trovare nulla che funzioni. L'alternativa sarà scrivere come estensione C, ma sembra che reinventerò la ruota.

Answer 1

numpy.maximum.accumulate funziona per me.

>>> import numpy 
>>> numpy.maximum.accumulate(numpy.array([11,12,13,20,19,18,17,18,23,21])) 
array([11, 12, 13, 20, 20, 20, 20, 20, 23, 23]) 

Answer 2

Come suggerito, c'è scipy.maximum.accumulate:

In [9]: x 
Out[9]: [1, 3, 2, 5, 4] 

In [10]: scipy.maximum.accumulate(x) 
Out[10]: array([1, 3, 3, 5, 5])