t, npu = np.broadcast_arrays(period, numevents_per_unit)