问题描述
我正在尝试使用if语句检查近似值何时收敛，但是它一直返回在2次迭代后收敛，这是不应该的，因为删除了if-break语句并将其设置为1000次迭代就可以了值非常不同。
def Jacobi2D(U,xrange,yrange,order):
Up=sym.zeros(len(xrange),len(yrange))
for k in range(1,order):
if k>=2 and U==Up:
print('Approximation converges after',k,'iterations')
break
Up=U #Save the U of the previous iteration for use in the formula
for i in range(1,len(xrange)-1):
for j in range(1,len(yrange)-1):
U[i,j]=Up[i,j]+(1/4)*(Up[i+1,j]+Up[i-1,j]+Up[i,j+1]+Up[i,j-1]-4*Up[i,j])
continue
continue
continue
return U
该函数将U作为维度为（len（xrange），len（yrange））的矩阵，并且“ order”是您希望近似值经历的迭代次数。有谁知道为什么在重新定义Up = U之前触发条件？任何帮助将不胜感激，谢谢！
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）