上传者: 38524472
|
上传时间: 2022-09-25 23:52:13
|
文件大小: 42KB
|
文件类型: PDF
import numpy as np
import time
1.1 Jacobi迭代算法
def Jacobi_tensor_V2(A,b,Delta,m,n,M):
start=time.perf_counter()#开始计时
find=0#用于标记是否在规定步数内收敛
X=np.ones(n)#迭代起始点
x=np.ones(n)#用于存储迭代的中间结果
d=np.ones(n)#用于存储Ax**(m-2)的对角线部分
m1=m-1
m2=2-m
for i in range(M):
print('X',X)
a=np.copy(A)
#得Ax**(m-2)
for j in ran