Невозможно привести данные массива из dtype('complex128') в dtype('float64') в соответствии с правилом "safe". ⇐ Python

[Anonymous]

Я пытаюсь решить задачу оптимизации, как показано ниже. Но каждый раз я получаю сообщение об ошибке. Невозможно привести данные массива из dtype('complex128') в dtype('float64') в соответствии с правилом "safe".
Может ли кто-нибудь помочь мне увидеть, что не так в коде?

def func(vec):

linspec = -(kx**2)+((1.-nu)*(kx**4))
lin = linspec*np.fft.fft(vec)
nlin = np.zeros_like(lin)
nlinre = vec*vec
nlinspec = np.fft.fft(nlinre)
nlin = (0.5*1j*kx*nlinspec)
sol = lin+nlin
rhs = np.zeros_like(sol, dtype='complex')
sol -= rhs
sol = np.fft.ifft(sol).real

return sol

def kssol(u0):

u1 = np.ones((2*Mx,), dtype='complex')
#u1 = 100.*u0
u = scipy.optimize.fsolve(func, u1)

return u


Подробнее здесь: https://stackoverflow.com/questions/491 ... cording-to