Is it possible to vectorize recursive calculation of a NumPy array where each element depends on the previous one?

by

You might think this would work: 

import numpy as np
n = len(Tm)
t = np.empty(n)

t[0] = 0  # or whatever the initial condition is 
t[1:] = Tm[1:] + (t[0:n-1] - Tm[1:])**(-tau[1:])

but it doesn’t: you can’t actually do recursion in numpy this way (since numpy calculates the whole RHS and then assigns it to the LHS).

So unless you can come up with a non-recursive version of this formula, you’re stuck with an explicit loop:

tt = np.empty(n)
tt[0] = 0.
for i in range(1,n):
    tt[i] = Tm[i] + (tt[i-1] - Tm[i])**(-tau[i])

More Related Contents:

  1. Taking subarrays from numpy array with given stride/stepsize
  2. What is vectorization? [closed]
  3. Bin elements per row – Vectorized 2D Bincount for NumPy
  4. Vectorized searchsorted numpy
  5. Vectorized NumPy linspace for multiple start and stop values
  6. Using Numpy Vectorize on Functions that Return Vectors
  7. Vectorizing Haversine distance calculation in Python
  8. Why np.hypot and np.subtract.outer very fast compared to vanilla broadcast ? Using Numba for speedup numpy in parallel for distance matrix calculation
  9. Vectorized look-up of values in Pandas dataframe
  10. Python vectorizing nested for loops
  11. Efficient evaluation of a function at every cell of a NumPy array
  12. python recursive vectorization with timeseries
  13. How to invert a permutation array in numpy
  14. vectorize conditional assignment in pandas dataframe
  15. Subtracting numpy arrays of different shape efficiently
  16. NumPy version of “Exponential weighted moving average”, equivalent to pandas.ewm().mean()
  17. Efficiently replace elements in array based on dictionary – NumPy / Python
  18. Concatenate range arrays given start, stop numbers in a vectorized way – NumPy
  19. numpy array TypeError: only integer scalar arrays can be converted to a scalar index
  20. Numpy Broadcast to perform euclidean distance vectorized
  21. Vectorized moving window on 2D array in numpy
  22. Vectorizing `numpy.random.choice` for given 2D array of probabilities along an axis
  23. Numpy: argmax over multiple axes without loop
  24. How can I tell if NumPy creates a view or a copy?
  25. Convert array of indices to one-hot encoded array in NumPy
  26. Py_initialize / Py_Finalize not working twice with numpy
  27. “In” operator for numpy arrays?
  28. Problem in combining bar plot and line plot (python)
  29. PIP Install Numpy throws an error “ascii codec can’t decode byte 0xe2”
  30. Screen Capture with OpenCV and Python-2.7

Leave a Comment