我正在尝试将参数传递给函数,我认为我正确地执行了它,但它仍然会给出错误: TypeError:p_vinet()需要2个位置参数,但是给出了4个
首先是我的代码片段然后我会给出产生错误的部分。
volumeMgO是先前计算的数组,包括:
array([ 7.64798549, 7.67153344, 7.67153344, 7.8068763 , 7.97288941,
8.14781986, 8.33321177, 8.53118834, 8.74433596, 8.97545339,
9.22826581, 9.50740563, 9.81962839])
params_MgO是这样的:
params_MgO = [11.244, 160., 4.0]
vinet功能是:
def p_vinet(v, params):
"""
This function will calculate pressure from Vinet equation.
Parameters
==========
v = volume
params = [V0, K0, K0']
Returns
=======
Pressure calculated from Vinet equation
"""
f_v = np.power( v/params[0], 1./3.)
eta = 1.5 * (params[2] - 1.)
P = 3 * params[1] * ((1 - f_v) / np.power(f_v, 2) ) * np.exp(eta * (1 -
f_v))
return P
最后,斜率函数只是获取衍生物的一种简单方法:
def slope(func, x, h, args=()):
"""
find a slope for function f at point x
Parameters
=========
f = function
x = independent variable for f
h = distance along x to the neighboring points
Returns
=======
slope
"""
rise = func(x+h, *args) - func(x-h, *args)
run = 2.*h
s = rise/run
return s
现在问题就出现了。当我输入:
BulkModulus_MgO = np.zeros(volumeMgO.size)
for i in range(volumeMgO.size):
BulkModulus_MgO[i] = slope(p_vinet, volumeMgO[i], volumeMgO[i]*0.0001,
args=(params_MgO))
我收到此错误:
---------------------------------------------------------------------------
TypeError Traceback (most recent call last)
<ipython-input-116-60467d989bbc> in <module>()
1 BulkModulus_MgO = np.zeros(volumeMgO.size)
2 for i in range(volumeMgO.size):
----> 3 BulkModulus_MgO[i] = slope(p_vinet, volumeMgO[i],
volumeMgO[i]*0.0001, args=(params_MgO))
<ipython-input-100-618f25e85d34> in slope(func, x, h, args)
15 """
16
---> 17 rise = func(x+h, *args) - func(x-h, *args)
18 run = 2.*h
19
TypeError: p_vinet() takes 2 positional arguments but 4 were given
我不明白。 p_vinet需要参数v和params,并且我通过斜率函数中的x + h和x-h提供v,并且params是包含p_vinet解包的3个条目的列表。这就是2个论点。为什么说我供应4?
很抱歉,如果我提交代码的方式有些混乱。我在jupyter笔记本中进行编码,所有功能都是分开的。 volumeMgO与以前的代码分开计算,没有任何问题。
答案 0 :(得分:0)
让我们看一下这行代码:
BulkModulus_MgO[i] = slope(p_vinet, volumeMgO[i], volumeMgO[i]*0.0001,
args=(params_MgO))
args应该是这样的:args=(params_MgO,)
(其中包含一个元素的元组)而不是args=(params_MgO)
(不是元组,只是3个元素的数组),因为在第二种情况下解包*args
slope()
函数为您提供了3个额外的参数(params_MgO
的每个元素)。这就是你在slope()
函数中得到4个参数的原因。在第一种情况下,解包为整个数组提供单个参数(如p_vinet()
函数中的params。)