Stack numpy array as output of complex function produces ValueError: order must be one of 'C', 'F', 'A', or 'K' (got 's0f')

I have found a workaround, using tt.set_subtensor() in a loop to assign output values to an empty array. Painfully slow but it worked.