June 24, 2019
12:05 p.m.
On 6/24/19 11:46 AM, Allan Haldane wrote:
A no-clobber guarantee makes your "iterative mask" example solvable in an efficient (no-copy) way:
mask, last_mask = False while True: dat_mean = np.mean(MaskedArray(data, mask)) mask, last_mask = np.abs(data - mask) > cutoff, mask if np.all(mask == last_mask): break
Whoops, that should read "np.abs(data - dat_mean)" in there. Allan