我正在尝试使用NumPy模拟任意精度的二进制算术。举一个简单的例子，我有一段代码（在基本的非NumPy Python中），它将一个二进制数加一，其中二进制数表示为从最低到最高有效位的0和1的列表（因此它们从右到左读取为二进制数）：

def increment_bits(bits):
    """
    Returns binary representation corresponding to bits + 1
    where bits is a list of 0s and 1s.
    
    >>> increment_bits([1, 1, 1, 0, 1])   # 23 + 1 == 24
    [0, 0, 0, 1, 1]
    >>> increment_bits([1, 1, 1])         #  7 + 1 ==  8   <- an extra bit now needed
    [0, 0, 0, 1]
    """
    new_bits = bits[:]
    for i, v in enumerate(new_bits):
        if v:
            new_bits[i] = 0
        else:
            new_bits[i] = 1
            return new_bits
    # if we have made it here, then there is a final "carry"
    new_bits.append(1)
    return new_bits

字符串
我的目标是使用NumPy实现相同的效果，但速度更快，但我是NumPy的新手，不确定迭代NumPy数组以实现所需效果的最快（或最NumPy）方法。