Re: [Numpy-discussion] best way of speeding up a filtering-like algorithm

2018-03-28 Thread Eric Wieser
Well, one tip to start with: numpy.where(some_comparison, True, False) is the same as but slower than some_comparison Eric On Wed, 28 Mar 2018 at 18:36 Moroney, Catherine M (398E) < catherine.m.moro...@jpl.nasa.gov> wrote: > Hello, > > > > I have the following sample code (pretty simple algor

Re: [Numpy-discussion] best way of speeding up a filtering-like algorithm

2018-03-28 Thread Joseph Fox-Rabinovitz
It looks like you are creating a coastline mask (or a coastline mask + some other mask), and computing the ratio of two quantities in a particular window around each point. If your coastline covers a sufficiently large portion of the image, you may get quite a bit of mileage using an efficient conv

[Numpy-discussion] PR adding support for object arrays to np.isinf, np.isnan, np.isfinite

2018-03-28 Thread Joseph Fox-Rabinovitz
I have opened PR #10820 to add support for `dtype=object` to `np.isinf`, `np.isnan`, `np.isfinite`. The PR is a fairly minor change, but I would like to make sure that I understand at least the basics of ufuncs before I start adding support for datetimes and timedeltas to `np.isfinite` and eventual