neurophysik / jitcdde

Just-in-time compilation for delay differential equations
Other
56 stars 14 forks source link

Facing issues while applying the jitcdde for a 0.2 sec time delay #9

Closed Naliniprasad closed 6 years ago

Naliniprasad commented 6 years ago

I continue my attempt to apply Jitcdde. Now I encounter a issue. While solving a DDE system with three equations and one constant delay. With a set of parameter values where time delay is 0.2 then the system is showing some math domain error. Can you please give me some suggestion regarding such problem. If you provide me your mail id i will send you the code. Here I am providing the error message.

File "", line 1, in runfile('C:/Users/user/Downloads/jitcode-master/jitcode-master/examples/untitled1.py', wdir='C:/Users/user/Downloads/jitcode-master/jitcode-master/examples')

File "C:\ProgramData\Anaconda3\lib\site-packages\spyder\utils\site\sitecustomize.py", line 710, in runfile execfile(filename, namespace)

File "C:\ProgramData\Anaconda3\lib\site-packages\spyder\utils\site\sitecustomize.py", line 101, in execfile exec(compile(f.read(), filename, 'exec'), namespace)

File "C:/Users/user/Downloads/jitcode-master/jitcode-master/examples/untitled1.py", line 22, in data.append( DDE.integrate(time) )

File "C:\ProgramData\Anaconda3\lib\site-packages\jitcdde_jitcdde.py", line 789, in integrate self.DDE.get_next_step(self.dt)

File "C:\ProgramData\Anaconda3\lib\site-packages\jitcdde_python_core.py", line 225, in get_next_step k_3 = self.eval_f(self.t + 0.75delta_t, self.y + 0.75delta_t*k_2)

File "C:\ProgramData\Anaconda3\lib\site-packages\jitcdde_python_core.py", line 218, in eval_f return self.f(t, y, *self.parameters)

File "C:\ProgramData\Anaconda3\lib\site-packages\jitcdde_python_core.py", line 177, in self.f = lambda args: np.array(F(args)).flatten()

File "", line 1, in

ValueError: math domain error

Wrzlprmft commented 6 years ago

I would need your code to give any meaningful help, preferably an minimal example. Usually there should be no problem with posting such a reduced example here, as it has little relation to your actual code. Still, if you do not want to do this, you can find my e-mail address here. Either way, please also tell me which version of JiTCODE you are using.