тройная интеграция по массивам python - PullRequest
Новая
       129

тройная интеграция по массивам python

0 голосов
/ 22 октября 2018

У меня есть следующий тройной интеграл:

интеграл

данные: t 1D массив (размер 124), pr 1D массив (размер 10),Lat 1D массив (241) и LON 1D массив (размер 480) V 4D массив (124, 10, 241, 480)

Я хочу интегрировать v по t, pr и lon в каждом лат. Итак, я использовалследующий код: 

def M(T, lam, P, V, phi):

    return integrate.tplquad(V*R*np.cos(phi), 0, T,
                         lambda T: 0, lambda T: lam,
                         lambda T,lam: 0, lambda T,lam: P)[0]

for i in range(lat.shape[0]):
    a = M(t, lon, pr, v[:, :, i, :], lat[i])

Но я получил ошибку:

ValueError: Значение истинности массива с более чем одним элементом неоднозначно.Используйте a.any () или a.all ()

То, как работает tplquad, до сих пор сбивает меня с толку.Любая помощь в выполнении интеграла с tplquad или других функций?

...