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

Finding limit cycles #11

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
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions