2013-10-28 5 views
3

Я понимаю, что (одно использование) vectorize numpy позволяет мне отправлять массив функции, которая обычно принимает только скаляры, вместо использования встроенной функции отображения (в сочетании с лямбда-функции или тому подобное). Однако в следующем сценарии у меня возникают разные результаты, когда я использую map vs numpy.vectorize, и я не могу понять, почему.Различные результаты между python map и numpy vectorize

import numpy as np 

def basis2(dim, k, x): 
    y = np.array([-0.2, -0.13, -0.06, 0, 0.02, 0.06, 0.15, 0.3, 0.8, 
        1.6, 3.1, 6.1, 10.1, 15.1, 23.1, 30.1, 35.0, 40.0, 45.0, 50.0, 55.0]) 

    if x < y[k] or x > y[k + dim + 1]: 
     return 0 

    elif dim != 0: 
     ret = ((x - y[k])/(y[k + dim] - y[k])) * basis2(dim - 1, k, x) + (
      (y[k + dim + 1] - x)/(y[k + dim + 1] - y[k + 1])) * basis2(dim - 1, k + 1, x) 
     return ret 

    else: 
     return 1.0 

w = np.array([20.0, 23.1, 30.0]) 
func = lambda x: basis2(3, 14, x) 
vec = map(func, w) 

func2 = np.vectorize(basis2) 
vec2 = func2(3, 14, w) 

print vec # = [0, 0.0, 0.23335417007039491] 
print vec2 # = [0 0 0] 

ответ

6

Как говорит строка документации:

Тип данных на выходе vectorized определяется путем вызова функции с первым элементом ввода. Этого можно избежать , указав аргумент otypes.

вам нужно добавить otypes аргумент:

func2 = np.vectorize(basis2, otypes="d") 

или изменить return 0 к return 0.0 в basis2().

Смежные вопросы