Работа с срезами памяти Cython и вещание - PullRequest
Новая
       7

Работа с срезами памяти Cython и вещание

0 голосов
/ 30 марта 2020

У меня возникли проблемы с расшифровкой этой ошибки Cython: 

---> 83     return RVKDE_loop(S, X, sigmas, n, m)

_cython_magic_e23143d4f8f8ff2526f1f1d1d47eacfb.pyx in _cython_magic_e23143d4f8f8ff2526f1f1d1d47eacfb.RVKDE_loop()

TypeError: unsupported operand type(s) for -: 'float' and '_cython_magic_e23143d4f8f8ff2526f1f1d1d47eacfb._memoryviewslice'

Вот код спагетти Cython, который выдал ошибку: 

%%cython
cimport cython
import numpy as np
from libc.stdio cimport printf

ctypedef fused dfloat:
    cython.float
    cython.double
    int


def RVKDE_loop(dfloat[:] S, dfloat[:] X, dfloat[:] sigmas, int n, int m):
    total = []
    for i in range(n):
        norms = -(abs(X[i]-S)**2)
        exp_denom = (2*(sigmas[i]**2))
        a_denom = ((sigmas[i]*np.sqrt(2*np.pi))**m)
        a = (1/a_denom)
        value = a*np.exp(norms/exp_denom)
        total += [value]
    return (1/n)*np.sum(total,axis=1)

Вот несколько примеров ввода 

S: np.array([-1.3015387, -0.07296862,  0.2387931,   0.38824359,  0.47182825,  0.75062962,
  1.3190391,   1.86540763,  2.62434536,  2.74481176])

X: np.array([-5.,         -3.77777778, -2.55555556, -1.33333333, -0.11111111,  1.11111111,
  2.33333333,  3.55555556,  4.77777778,  6.        ])

sigmas: np.array([2.43158064, 0.59796393, 0.27334147, 0.14160981, 0.14160981, 0.53204325,
 1.06717756, 1.06717756, 0.21537329, 0.21537329])

n: 10
m: 1

Любые советы о том, как я могу go исправить это, были бы великолепны. Обратите внимание, что функция cython делает то, что я хочу, когда это обычный python код; это просто слишком медленно на мой вкус в обычном python.

...