Ejecución del máximo de los valores numpy de la matriz

Question

Necesito una forma rápida de mantener un máximo de ejecución de una matriz numpy. Por ejemplo, si mi arsenal era:

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

me gustaría:

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

Obviamente, podría hacer esto con un poco de bucle:

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 

Pero mis matrices tener cientos de miles de entradas y necesito llamar esto muchas veces. Parece que debe haber un truco numpy para eliminar el loop, pero parece que no puedo encontrar nada que funcione. La alternativa será escribir esto como una extensión C, pero parece que estaría reinventando la rueda.

Answer 1

numpy.maximum.accumulate funciona para mí.

>>> 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

Como se ha sugerido, hay 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])