Библиотека Python для итеративного решателя Гаусса-Зейделя?

Я знаю, что это старо, но я не нашел какой-либо существующей библиотеки в python для gauss-seidel. Вместо этого я создал свою маленькую функцию, которая с помощью матрицы перестановок, как показано в другом моем ответе матрица перестановки создаст решение (вектор x) для любой квадратной матрицы, в том числе с нулями по диагонали.

def gauss_seidel(A, b, tolerance, max_iterations, x):
#x is the initial condition


iter1 = 0
#Iterate
for k in range(max_iterations):
    iter1 = iter1 + 1
    print ("The solution vector in iteration", iter1, "is:", x)    
    x_old  = x.copy()
    
    #Loop over rows
    for i in range(A.shape[0]):
        x[i] = (b[i] - np.dot(A[i,:i], x[:i]) - np.dot(A[i,(i+1):], x_old[(i+1):])) / A[i ,i]
        
    #Stop condition 
    #LnormInf corresponds to the absolute value of the greatest element of the vector.
    
    LnormInf = max(abs((x - x_old)))/max(abs(x_old))   
    print ("The L infinity norm in iteration", iter1,"is:", LnormInf)
    if  LnormInf < tolerance:
        break

       
return x

После повторного поиска я нашел следующее, что можно было бы использовать в качестве готового пакета, но сам я его не использовал числовой PyPI