<a href="https://github.com/numpy/numpy/pull/99">https://github.com/numpy/numpy/pull/99</a><div><br></div><div><div>In [24]: a = np.ones((1000,1000))</div><div>In [25]: b = np.ones((1000,1000))</div><div>In [26]: c = np.zeros((1000, 1000))</div>
<div>In [27]: m = np.random.rand(1000,1000) &gt; 0.5</div><div><br></div><div>In [28]: timeit c[m] = a[m] + b[m]</div><div>1 loops, best of 3: 246 ms per loop</div><div><br></div><div>In [29]: timeit np.add(a, b, out=c, where=m)</div>
<div>10 loops, best of 3: 20.3 ms per loop</div><div><br></div><div>-Mark</div></div>