SciPy函数调用之间的根值不一致

问题描述

使用opt.root,我想存储调用函数的参数。

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()位)