Ответ @ Скотта правильный. Однако я искал инкрементную версию, которую, я думаю, я нашел:
Определите results = np.empty((0, output_shape))
, а затем обновите его в цикле, используя results = np.concatenate((results, some_func(x)))
Я не уверен, как мне следует думать о размере 0 в numpy, но это работает.