Ошибка при использовании лямбда-функции с Cython - PullRequest
Новая
       5

Ошибка при использовании лямбда-функции с Cython

6 голосов
/ 28 ноября 2010

Я пытаюсь использовать Cython для ускорения кода. Cython выдает ошибку, которая гласит «Ожидается идентификатор или литерал», когда я использую лямбда-функции. Насколько я могу судить, лямбда-функции предназначены для поддержки в Cython 0.13. Я не прав в этом вопросе? Если они действительно поддерживаются, нужно ли мне что-то делать, кроме того, что я имею здесь, чтобы реализовать их? 

def f(e_1, e_2, rho):
    """Bivariate Normal pdf with mean zero, unit variances, and correlation coefficient rho."""
    return (1.0 / (2.0 * pi * sqrt(1 - rho**2))) * exp(-(1.0 / (2*(1 - rho**2))) * (e_1**2 + e_2**2 - 2*rho*e_1*e_2))

def P_zero(b_10, b_11, b_20, b_21, rho, gamma, x):
    """Returns the probability of observing zero entrants in a market by numerically
    integrating out the unobserved firm-specific profit shocks."""
    h_z = lambda e_1: -inf
    g_z = lambda e_1: -b_10 - b_11*x[0] - gamma*x[1]
    I   = lambda e_1, e_2: f(e_1, e_2, rho)
    return dblquad(I, -inf, (-b_20 - b_21*x[0] - gamma*x[2]), h_z, g_z)[0]

Ответы [ 2 ]

2 голосов
/ 28 ноября 2010

По моему мнению, вы должны изменить h_z = lambda e_1: -inf на h_z = lambda e_1: -float('inf'), если вы не определили inf где-то еще.

0 голосов
/ 19 февраля 2011

Я могу скомпилировать приведенный ниже упрощенный код Cython и нормально его запустить, используя Cython 0.14.1 на OS X 10.6.6.Я не знаю деталей относительно того, почему это не работает на 0.13.Самое простое решение - обновить Cython, если это возможно. 

def f(e_1, e_2, rho):
    return e_1 + e_2 + rho

def dummy(a, b, c, d, e):
    return [a(1,2) + b + c + d(1) + e(3)]

def P_zero(b_10, b_11, b_20, b_21, rho, gamma, x):
    h_z = lambda e_1: -1000
    g_z = lambda e_1: -b_10 - b_11 * x[0] - gamma * x[1]
    I   = lambda e_1, e_2: f(e_1, e_2, rho)
    return dummy(I, -1000, (-b_20 - b_21 * x[0] - gamma * x[2]), h_z, g_z)[0]

print P_zero(1, 2, 3, 4, 5, 6, [6, 7, 8]) 
# outputs "-2122"
...