问题描述
import scipy.optimize as opt
x0=[0,0]
prev_x=x0
count=0
def fun(x):
global prev_x,count
print('---')
print('iteration',count)
count+=1
print('last iteration',prev_x)
print('this iteration',x)
prev_x=x
return x[0]**2*x[1]+1,x[1]-x[0]+2
result=opt.root(fun,x0,method='lm')
我希望对于每个迭代步骤,保存在last_x中的值都应与上一个迭代的当前值匹配。 运行上面的代码,我得到例如之间的不同值迭代8和9。
我是否缺少有关Python和/或SciPy行为的基本知识,或者这是一个错误吗?
[edit]我把问题简化到了重点。
解决方法
分配列表不会在整个值之间复制值,只是给它提供了不同的名称,因此通过一个名称对列表所做的更改会反映在另一个名称上。要获得您期望的行为,请将相关行更改为
prev_x=x.copy()
(请注意copy()
位)