Komplexe Typkonvertierungsfehler in Numpy vs c Complex.h Library ⇐ Python

[Guest]

Ich habe eine große Python -Methode, die eine Reihe von NDarray -Mengen von dType = np.complex128 berechnet. Ich habe mit Cython mit der Erhöhung der Geschwindigkeit experimentiert. Bei der Umwandlung einer Python -Methode in Cython wird der Überblick über die Verwendung von Numpy -Berechnungen im Vergleich zu C -Bibliotheken, um dies zu tun. Bei der Erforschung der letzteren ist dies eine verwässerte Cython-Methode zur Berechnung des komplexen Exponentials eines Eingangsnumpy-Arrays, x: < /p>
import numpy as np
cimport numpy as cnp
cnp.import_array()

cdef extern from "complex.h":
double complex cexp(double complex z)

DTYPEC = np.complex128
ctypedef cnp.complex128_t DTYPEC_t

cpdef cnp.ndarray[DTYPEC_t, ndim=1] stupid_exp(cnp.ndarray[DTYPEC_t, ndim=1] x):

cdef int num_pts, j
cdef double complex tmp

num_pts = len(x)
cdef cnp.ndarray[DTYPEC_t, ndim=1] result = np.zeros(num_pts, dtype=DTYPEC)

for j in range(num_pts):
result[j] = cexp(x[j])

return result
< /code>
Auf meinem Windows 10 -Gerät, kompiliert mit Visual Studio 2022, erstellt das Ergebnis [j] = CEXP (x [j]) eine Fehlerpaar: < /p>
error C2440: 'function': cannot convert from '__pyx_t_double_complex' to '_Dcomplex'
error C2440: '=': cannot convert from '_Dcomplex' to '__pyx_t_double_complex'
< /code>
Dies besagt, dass Numpy Complex128 und Cs Doppelkomplex nicht als der gleiche Datentyp angesehen werden. Ein LLM schlug eine explizite Typumwandlung vor: < /p>
tmp = x[j]
tmp = cexp(tmp)
result[j] = tmp
< /code>
, aber dies behebt die Fehler nicht. Irgendwelche Gedanken darüber, wie ich diesen einfachen Code in Cython richtig kompilieren kann?