Function update_column

utils/creatematrix.py:298–299

(col)

        print("updating A")

        def update_column(col):
            A[1:, col] -= v2 * (a0.T @ A[1:, col]) / S

        if S == 0:
            raise ValueError("S is 0")