You just write `np.full_like(condlist[0], default)` and it does the same thing as the multiplication, but about twice as fast. And works even for non-numeric types. The latter is probably a bug in the current code -- you might want to add a test that uses string arrays or something.