-
-
Notifications
You must be signed in to change notification settings - Fork 9
Finding limit cycles #11
Description
Hi everyone,
I am relatively new to Julia and I am experiencing an error, that I cannot resolve using the PyDSTool wrapper. I do not know, whether my code is faulty or if something inside my PyDSTool installation is acting up.
Problem description: I have a set of nonlinear differential equations and I want to find stable, periodic solutions, "limit cycles" I suppose. I have two bifurcation parameters, f_N and ome_N. Following the instructions at http://docs.juliadiffeq.org/latest/analysis/bifurcation.html:
u0 = zeros(18);
tspan = [0;30]
p = [0,0.3*2*pi*952.26]
dsargs = build_ode(EOM_nd_9p,u0,tspan,p)
ode = ds[:Generator][:Vode_ODEsystem](dsargs)
ode[:set](pars = Dict("f_N"=>0.0))
ode[:set](pars = Dict("ome_N"=>0.3*2*pi*952.26))
PC = ds[:ContClass](ode)
Everything works fine up to this point. Now I want to calculate a limit cycle curve:
bif = bifurcation_curve(PC,"LC-C",["f_N","ome_N"],
max_num_points=450,
max_stepsize=2,min_stepsize=1e-5,
stepsize=2e-2,loc_bif_points="all",
save_eigen=true,name="LC1",
print_info=true,calc_stab=true)
And this error is thrown:
Error opening file Shell9{getfield(Main, Symbol("##3#7")),getfield(Main, Symbol("##4#8")),getfield(Main, Symbol("##5#9")),Nothing,Nothing,getfield(Main, Symbol("##6#10")),Expr,Expr}_vf.c for writing
ERROR: PyError ($(Expr(:escape, :(ccall(#= C:\Users\Martin.julia\packages\PyCall\rUul9\src\pyfncall.jl:44 =# @pysym(:PyObject_Call), PyPtr, (PyPtr, PyPtr, PyPtr), o, pyargsptr, kw))))) <class 'OSError'>
OSError(OSError(22, 'Invalid argument'),)
File "C:\Users\Martin.julia\packages\Conda\hsaaN\deps\usr\lib\site-packages\PyDSTool\PyCont\ContClass.py", line 212, in newCurve
self.loadAutoMod()
File "C:\Users\Martin.julia\packages\Conda\hsaaN\deps\usr\lib\site-packages\PyDSTool\PyCont\ContClass.py", line 423, in loadAutoMod
self.makeAutoLibSource()
File "C:\Users\Martin.julia\packages\Conda\hsaaN\deps\usr\lib\site-packages\PyDSTool\PyCont\ContClass.py", line 611, in makeAutoLibSource
raise IOError(e)
I already did some investigating, but I did not find any solution. Maybe some module imported from the AUTO97 software is not working.
Thanks in advance!
Martin
System of differential equations attached.
EOM_nd_9p.txt