Please cast to correct type. You probably need something like `dtype = np.result_type(default, *choicelist)` and then use that dtype to create all result arrays. This probably solves already the bug Julian found.