what is the quickest way to iterate through a numpy array
Solution 1
These are my timings on a slower machine
In [1034]: timeit [i for i in np.arange(10000000)]
1 loop, best of 3: 2.16 s per loop
If I generate the range directly (Py3 so this is a genertor) times are much better. Take this a baseline for a list comprehension of this size.
In [1035]: timeit [i for i in range(10000000)]
1 loop, best of 3: 1.26 s per loop
tolist
converts the arange to a list first; takes a bit longer, but the iteration is still on a list
In [1036]: timeit [i for i in np.arange(10000000).tolist()]
1 loop, best of 3: 1.6 s per loop
Using list()
- same time as direct iteration on the array; that suggests that the direct iteration first does this.
In [1037]: timeit [i for i in list(np.arange(10000000))]
1 loop, best of 3: 2.18 s per loop
In [1038]: timeit np.arange(10000000).tolist()
1 loop, best of 3: 927 ms per loop
same times a iterating on the .tolist
In [1039]: timeit list(np.arange(10000000))
1 loop, best of 3: 1.55 s per loop
In general if you must loop, working on a list is faster. Access to elements of a list is simpler.
Look at the elements returned by indexing.
a[0]
is another numpy
object; it is constructed from the values in a
, but not simply a fetched value
list(a)[0]
is the same type; the list is just [a[0], a[1], a[2]]]
In [1043]: a = np.arange(3)
In [1044]: type(a[0])
Out[1044]: numpy.int32
In [1045]: ll=list(a)
In [1046]: type(ll[0])
Out[1046]: numpy.int32
but tolist
converts the array into a pure list, in this case, as list of ints. It does more work than list()
, but does it in compiled code.
In [1047]: ll=a.tolist()
In [1048]: type(ll[0])
Out[1048]: int
In general don't use list(anarray)
. It rarely does anything useful, and is not as powerful as tolist()
.
What's the fastest way to iterate through array - None. At least not in Python; in c code there are fast ways.
a.tolist()
is the fastest, vectorized way of creating a list integers from an array. It iterates, but does so in compiled code.
But what is your real goal?
Solution 2
This is actually not surprising. Let's examine the methods one a time starting with the slowest.
[i for i in np.arange(10000000)]
This method asks python to reach into the numpy array (stored in the C memory scope), one element at a time, allocate a Python object in memory, and create a pointer to that object in the list. Each time you pipe between the numpy array stored in the C backend and pull it into pure python, there is an overhead cost. This method adds in that cost 10,000,000 times.
Next:
[i for i in np.arange(10000000).tolist()]
In this case, using .tolist()
makes a single call to the numpy C backend and allocates all of the elements in one shot to a list. You then are using python to iterate over that list.
Finally:
list(np.arange(10000000))
This basically does the same thing as above, but it creates a list of numpy's native type objects (e.g. np.int64
). Using list(np.arange(10000000))
and np.arange(10000000).tolist()
should be about the same time.
So, in terms of iteration, the primary advantage of using numpy
is that you don't need to iterate. Operation are applied in an vectorized fashion over the array. Iteration just slows it down. If you find yourself iterating over array elements, you should look into finding a way to restructure the algorithm you are attempting, in such a way that is uses only numpy operations (it has soooo many built-in!) or if really necessary you can use np.apply_along_axis
, np.apply_over_axis
, or np.vectorize
.
piRSquared
Finance professional turned python/pandas enthusiast. Canonicals How to Pivot a Pandas DataFrame Q|A Pandas concat Q|A Operate with a Series on every column of a DataFrame Q|A How Does Numpy's axis argument work with sum A Support If you've found my questions and/or answers helpful, feel free to show your support or appreciation by up-voting. If you are looking for more of my content, you can explore below or visit a few of my selected answers. Renaming Columns Via Pipeline Selecting Rows From a DataFrame Editorial How to write a good question Underrated https://stackoverflow.com/a/57014212/2336654 #SOreadytohelp
Updated on October 06, 2021Comments
-
piRSquared over 2 years
I noticed a meaningful difference between iterating through a numpy array "directly" versus iterating through via the
tolist
method. See timing below:directly
[i for i in np.arange(10000000)]
viatolist
[i for i in np.arange(10000000).tolist()]
considering I've discovered one way to go faster. I wanted to ask what else might make it go faster?
what is fastest way to iterate through a numpy array?
-
Ébe Isaac over 7 yearsThat is odd. I tried it myself several times and it seems that converting it to list does make it faster all the time. Thanks for bringing this to the light.
-
Divakar over 7 yearsJust iterate and get the list or do some processing too? Using just
list(np.arange(1000000))
looks quite fast. -
piRSquared over 7 years@Divakar see stackoverflow.com/a/40575522/2336654
-
Divakar over 7 yearsBut then
np.arange(1000000).tolist()
gives the same thing aslist(np.arange(1000000))
, so maybe my earlier comment isn't quite the expected thing I guess. -
piRSquared over 7 yearsI rolled back my edit because that is very fast to get at the list. But I still have to iterate through it and do processing.
-
Ignacio Vergara Kausel over 7 yearsMy question is why would you want to iterate over a numpy array instead of using vectorized functions.
-
hpaulj over 7 years
list()
produces a list ofnp.int32
objects;tolist
produces a list ofint
. They are not the same. -
piRSquared over 7 years@IgnacioVergaraKausel because I can't figure out a fast vectorized O(n) method. I'll post a question about it later today.
-
hpaulj over 7 yearsWhat's the goal this iteration? Just generating a list of integers?
tolist
is the fastest way. Applying some scalar function to each element of the array? -
piRSquared over 7 years@hpaulj I'm trying to calculate the cumulative count by unique value. See my answer to another question here stackoverflow.com/a/40575522/2336654. In this problem, I iterate through the array and track how many times I've seen an item, returning the count with every iteration. I've been trying to vectorize this. In fact, one of my answers in that link is an O(n^2) vectorized solution. I've explored the different return options of
np.unique
and haven't come up with a satisfactory answer. All this is information I'd include in another question. This was a by product of that one. -
hpaulj over 7 yearsIf you are doing something complicated at each step, the outer iteration mechanism doesn't make much difference in the time. Regardless of how you iterate you are repeating that costly step N-thousand times.
-
piRSquared over 7 years@hpaulj it's not at all complicated.
-
-
MaxNoe over 7 yearsBut there is a subtle difference between
list(np.arange(10))
andnp.arange(10).tolist()
: the first will result in a list ofnp.int64
the second in a list of pythonint
s. The first can be problematic for doing stuff like serialisation, e.g. using json. json will error on the first because it cannot handlenp.int64
-
piRSquared over 7 yearsThis is very useful and is why I've upvoted it and I hope others do to. I'm leaving the question open for now as I'm still left wanting to see other options of iteration through the array.
-
piRSquared over 7 yearsThanks @hpaulj this comes very close to actually answering my question in that you stated... "What's the fastest way to iterate through array - None." I'll likely be selecting this as my answer, but I'm leaving it open for a bit.
-
minhle_r7 over 2 years
v2
is unused, why would you need it? I think the difference you get is a flux, plus most of the time is spent on printing, not accessing the data.