if np.any(r12):
        out[r12] = betaln(a[r12], b[r12])