python-Lambda函数未返回正确的值
内容导读
互联网集市收集整理的这篇技术教程文章主要介绍了python-Lambda函数未返回正确的值,小编现在分享给大家,供广大互联网技能从业者学习和参考。文章包含2828字,纯文字阅读大概需要5分钟。
内容图文
我试图使吉莱斯皮算法变型,并确定反应倾向.我试图使用lambda表达式自动生成倾向向量.但是,在创建SSA.P时,一切都会出错.代码块中的最后一个循环PROPLOOP返回两个倾向,其中使用P_alternative生成的倾向是正确的.问题是:如何获得与SSA.P_alternative相同的SSA.P值?
import numpy as np
from numpy.random import uniform
class Markov:
def __init__(self,z0,t0,tf,rates,stoich):
self.S=stoich
self.z0=z0
self.rates=rates
self.P=self.propensities()
self.P_alternative=[
lambda z,rate:(0.5*rate[0]*z[0]*(z[0]-1)),
lambda z,rate:rate[1]*np.prod(z[0]),
lambda z,rate:rate[2]*np.prod(z[1]),
lambda z,rate:rate[3]*np.prod(z[1]),
lambda z,rate:rate[4]*np.prod(z[np.array([0,1])]),
lambda z,rate:rate[5]]
self.t0=t0
self.tf=tf
def propensities(self):
prop=[]
for i,reac in enumerate(self.S.T):
if all(z>=0 for z in reac):
prop.append(lambda z,rate:rate[i])
if any(z==-1 for z in reac):
j=np.where(reac==-1)[0]
prop.append(lambda z,rate:rate[i]*np.prod(z[j]))
if any(z==-2 for z in reac):
j=np.where(reac==-2)[0][0]
prop.append(lambda z,rate:(0.5*rate[i]*z[j]*(z[j]-1))[0])
return prop
stoich=np.array([
[-2, -1, 2, 0, -1, 0],
[ 1, 0, -1, -1, -1, 1],
[ 0, 0, 0, 1, 1, 0]])
rates=np.array([1.0,0.02,200.0,0.0004,0.9,0.9])
z0=np.array([540,730,0])
SSA=Markov(z0=z0,t0=0,tf=100,rates=rates,stoich=stoich)
#PROPLOOP; the values should be equal for both SSA.P and SSA.P_alternative, where SSA.P_alternative is the correct one
for i in xrange(len(SSA.P)):
print "Inexplicably wrong",SSA.P[i](z0,rates)
print "Correct answer",SSA.P_alternative[i](z0,rates), "\n"
输出为:
Inexplicably wrong 130977.0
Correct answer 145530.0
Inexplicably wrong 354780.0
Correct answer 10.8
Inexplicably wrong 354780.0
Correct answer 146000.0
Inexplicably wrong 354780.0
Correct answer 0.292
Correct answer 354780.0
Correct answer 354780.0
Inexplicably wrong 0.9
Correct answer 0.9
解决方法:
问题是您要在循环中创建lambda函数,并且它们引用的变量i和j可能随着循环的进行而改变.
lambda在创建时不会复制i或j的值,它只是保留对定义它们的名称空间的引用.当它在稍后调用时使用变量时,它将在该名称空间中查找它们.由于您的lambda在循环(实际上是整个函数)结束后被调用,因此它们都可以看到变量给出的最终值,这不是您想要的.这就解释了为什么两个版本的代码在上一次迭代中给出相同的输出. i和j的最终值是最后一个lambda函数的期望值.
您可以通过使lambda在定义时保留i和j当前值的副本来解决此问题.最简单的方法是使用默认参数:
for i,reac in enumerate(self.S.T):
if all(z>=0 for z in reac):
prop.append(lambda z, rate, i=i: rate[i]) # add i=i here and further down
if any(z==-1 for z in reac):
j=np.where(reac==-1)[0]
prop.append(lambda z, rate, i=i, j=j: rate[i]*np.prod(z[j]))
if any(z==-2 for z in reac):
j=np.where(reac==-2)[0][0]
prop.append(lambda z, rate, i=i, j=j: (0.5*rate[i]*z[j]*(z[j]-1))[0])
lambda定义中的i = i(必要时还有j = j)使lambda函数的变量参数具有默认值,该默认值是外部命名空间中i(和j)的当前值.由于调用lambda函数时仅传递两个参数,因此将使用保存的默认值.
内容总结
以上是互联网集市为您收集整理的python-Lambda函数未返回正确的值全部内容,希望文章能够帮你解决python-Lambda函数未返回正确的值所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错,欢迎将互联网集市网站推荐给程序员好友。
内容备注
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 gblab@vip.qq.com 举报,一经查实,本站将立刻删除。
内容手机端
扫描二维码推送至手机访问。