0

我有一个应该使用 ode 求解器解决的函数。问题是它有很多常数（比如 100），同时我不想把它们放在一个数组中，因为很难理解公式（这是化学的东西）。我在这个函数中的某些部分代码如下所示：

 # Component: Calcium Fluxes

    J_rel =  v1*(P_O1+P_O2)*(Ca_JSR - Ca_ss)*P_RyR # (micromolar_per_millisecond)
    J_tr = (Ca_NSR - Ca_JSR)/tau_tr #  (micromolar_per_millisecond)
    J_leak =  v2*(Ca_NSR - Ca_i) #  (micromolar_per_millisecond)
    J_up = ( v3*(Ca_i ^ 2.00000))/((K_m_up ^ 2.00000)+(Ca_i ^ 2.00000)) #  (micromolar_per_millisecond)
    J_xfer = (Ca_ss - Ca_i)/tau_xfer #  (micromolar_per_millisecond)
    P_RyR_prime =  - 0.0400000*P_RyR -  (( 0.100000*I_Ca_channels)/i_CaL_max)*(exp(( - ((V_m - 5.00000) ^ 2.00000)/648.000))) #  (micromolar_per_millisecond)

如果我想把它们放在一个数组中，它会变成这样：

J_rel =  algebraic[2]*(constants[21]+constants[3])*(constants[4] - Ca_ss)*P_RyR # (micromolar_per_millisecond)
J_tr = (Ca_NSR - Ca_JSR)/tau_tr #  (micromolar_per_millisecond)
J_leak =  constants[31]*(Ca_NSR - Ca_i) #  (micromolar_per_millisecond)
J_up = ( constants[69]*(Ca_i ^ 2.00000))/((constants[17] ^ 2.00000)+(Ca_i ^ 2.00000)) #  (micromolar_per_millisecond)
J_xfer = (Ca_ss - Ca_i)/tau_xfer #  (micromolar_per_millisecond)
P_RyR_prime =  - 0.constants[5]*P_RyR - (( 0.100000*constants[65])/constants[43])*(exp(( - ((V_m - constants[4]) ^ 2.00000)/constants[34]))) #  (micromolar_per_millisecond)

很难“看到”实际的公式。我必须将常量定义为全局变量吗？此函数由 ode 求解器使用，并被迭代约 50000 次。所以我需要一些有效的东西。

PS我正在使用python（scipy库），但我认为其他语言也可能有同样的问题。

非常感谢

4

3 回答 3

2

再想一想，为什么不直接使用全局变量呢？它会给你一个漂亮、简单和快速的结果。在 python 中，全局变量将被限制在定义它们的模块中，因此您不必担心它们会像在 C 中那样污染整个程序。

于 2011-03-03T22:00:42.343 回答
1

您可以使用函数关键字参数将字典“解包”到本地命名空间中：

c = { 'Ca_NSR' : 1.5, ... }

def f(Ca_JSR, Ca_ss, ...):
  J_rel =  v1*(P_O1+P_O2)*(Ca_JSR - Ca_ss)*P_RyR

f(**c)

有关f(**c) 含义的详细信息，请参阅有关关键字参数的文档。

于 2011-03-03T21:53:06.280 回答
1

如果你想避免污染命名空间，你可以使用字典：

c = { 'Ca_NSR' : 1.5, ... }

J_rel =  v1*(P_O1+P_O2)*(c['Ca_JSR'] - c['Ca_ss'])*P_RyR

你也可以使用一个类：

def Constants(object):
    def __init__(self):
        self.Ca_NSR = 1.5
        # etc...

c = Constants()
J_rel =  v1*(P_O1+P_O2)*(c.Ca_JSR - c.Ca_ss)*P_RyR
于 2011-03-03T21:45:55.583 回答