二元有限域矩阵的秩的计算 python语言实现

import numpy as np


def gfrank（data）:
    b = np.array（data）
    b = b[b[: 0].argsort（）]
    M = len（b）
    print（b[::-1]）
    b = b[::-1]
    l = 0
    k = 0
    while l < M:
        t = 0
        if b[k][l] == 0:
            for i in range（k + 1 M）:
                if b[i][l] == 1:
                    b[[k i] :] = b[[i k] :]
                    t = t + 1