>>671
mask = (a!=0).sum(axis=2) >= 1
mask = np.expand_dims(mask, 2)
org = mask*a + (~mask)*org