Skip to content
This repository was archived by the owner on Jun 24, 2022. It is now read-only.
This repository was archived by the owner on Jun 24, 2022. It is now read-only.

Finding limit cycles #11

@PotatoDynamics

Description

@PotatoDynamics

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

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions