python实现gauss-seidel迭代公式_多线性方程组迭代算法——Gauss-Seidel迭代算法的Python实现...

多线性方程组(张量)迭代算法的原理请看这里：原理部分请留言，不方便公开分享

import numpy as np

import time

1.1 Gauss-Seidel迭代算法

def GaussSeidel_tensor_V2(A,b,Delta,m,n,M):

start=time.perf_counter()

find=0

X=np.ones(n)

d=np.ones(n)

m1=m-1

m2=2-m

for i in range(M):

print('X',X)

x=np.copy(X)

#迭代更新

for j in range(n):

a=np.copy(A)

for k in range(m-2):

a=np.dot(a,x)

for k in range(n):

d[k]=a[k,k]

a[k,k]=m2*a[k,k]

x[j]=(b[j]-np.dot(a[j],x))/(m1*d[j])

#判断是否满足精度要求

if np.max(np.fabs(X-x))

find=1

break

X=np.copy(x)

end=time.perf_counter()

print('时间：',end-start)

print('迭代',i)

return X,find,i,end-start

1.2张量A的生成函数和向