b1, b2, r = np.broadcast_arrays(b1, b2, r)