Рекурсивное перемножение матриц списка с условием

def mult_matr(a,b):
    # Количество столбцов a должно быть равно к-ву строк b
    ra=len(a)
    ca=len(a[0])
    rb=len(b)
    cb=len(b[0])
    if ca == rb:
        # создадим матрицу-результат
        r = [[0 for _ in range(cb)] for _ in range(ra)]
        for i in range(ra):
            for j in range(cb):
                s=0
                for k in range(ca):
                    s=s+a[i][k]*b[k][j]
                r[i][j]=s
        return r        
    else:
        return None

# Умножение матриц
 
def mult_matr(a,b):
    # Количество столбцов a должно быть равно к-ву строк b
    ra=len(a)
    ca=len(a[0])
    rb=len(b)
    cb=len(b[0])
    if ca == rb:
        # создадим матрицу-результат
        r = [[0 for _ in range(cb)] for _ in range(ra)]
        for i in range(ra):
            for j in range(cb):
                s=0
                for k in range(ca):
                    s=s+a[i][k]*b[k][j]
                r[i][j]=s
        return r        
    else:
        return None
        
# Сравнение матриц        
        
def comp_matr(a,b,eps):
    n=len(a)
    for row in range(n):
        for col in range(n):
            if abs(a[row][col]-b[row][col]) > eps:
                return False
    return True            
 
# Главный цикл 
    
def main_loop(matr_list,n):
    
    sz=len(matr_list)
 
    if sz == 60:
        return matr_list
        
    if n >= sz:
        return main_loop(matr_list,0)
        
    q=matr_list[n]
    tmp=[]
    
    for v in matr_list:
        
        w=mult_matr(v,q)
        
        p=0
        
        for u in matr_list:
            if comp_matr(u,w,1.0e-5):
                p=1
                break
        
        if p==0:
            tmp.append(w)
            
    return main_loop(tmp,n+1)        
    
    
def start():
    
    a=[ [[1,2,3],[4,5,6],[7,8,9]] , [[-1,-2,-3],[-4,-5,-6],[-7,-8,-9]] , [[0,2,3],[4,0,6],[7,8,0]] ]
    
    print(main_loop(a,0))
    
start()

from random import randint
 
def get_rand_matr(n):
    return [[randint(1, 9) for j in range(n)] for i in range(n)]
 
def mult_matrices(a, b):
    n = len(a)
    res = [[0] * n for _ in range(n)]
    for i in range(n):
        for j in range(n):
            for k in range(n):
                res[i][j] += a[i][k] * b[k][j]
    return res
 
def eq_matr(a, b):
    for i in range(len(a)):
        if a[i] != b[i]:
            return False
    return True
 
def matr_in_arr(a, arr):
    for k in range(len(arr)):
        if eq_matr(arr[k], a):
            return True
    return False
 
arr = [get_rand_matr(3) for _ in range(3)]
for i in range(len(arr)):
    print(f'# {i+1}')
    print(*arr[i], sep='\n')
    print()
 
try:
    i = 0
    while True:
        for j in range(len(arr)):
            a = mult_matrices(arr[i], arr[j])
            if not matr_in_arr(a, arr):
                arr.append(a)
                print(f'# {len(arr)}, i = {i}, j = {j}')
                print(*arr[-1], sep='\n')
                print()
                if len(arr) >= 60:
                    raise Exception
        i += 1
except:
    print('Конец.')

from math import isclose
 
def mult_matrices(a, b):
    n=len(a)
    res=[[0.0] * n for _ in range(n)]
    for i in range(n):
        for j in range(n):
            for k in range(n):
                res[i][j]+=a[i][k] * b[k][j]
    return res
 
def eq_matr(a, b):
    for i in range(len(a)):
        for j in range(len(a)):
            if not isclose(a[i][j], b[i][j]):
                return False
    return True
 
def matr_in_arr(a, arr):
    for k in range(len(arr)):
        if eq_matr(arr[k], a):
            return True
    return False
 
def input_matr(n, name):
    print(f'\nВведите матрицу {name} порядка {n}:')
    res = []
    for i in range(n):
        res.append([float(x) for x in input(f'\tстрока {i+1} (значения через пробел): ').split()])
    return res
 
arr = []
for i in range(3):
    arr.append(input_matr(3, f'{i+1}'))
 
for i in range(len(arr)):
    print(f'# {i + 1}')
    print(*arr[i], sep='\n')
    print()
 
try:
    i = 0
    while True:
        for j in range(len(arr)):
            a = mult_matrices(arr[i], arr[j])
            if not matr_in_arr(a, arr):
                arr.append(a)
                print(f'# {len(arr)}, i = {i}, j = {j}')
                print(*arr[-1], sep='\n')
                print()
                if len(arr) >= 60:
                    raise Exception
        i += 1
except:
    print('Конец.')