Comments (3)
Thank you for spotting this.
Apart from that, there are many problems with your code. The first fundamental one is that you want to compute periodic orbits of non-autonomous systems. This is not possible in BifurcationKit.jl. Hence, you need to apply special trick to achieve this. You change the problem into:
function F!(f, x, p, t=0)
@unpack a, b1, b2, d, λ, ω = p
f[1] = x[2]
f[2] = - d * x[2] - a * x[1] - b1 * (x[1]^2) - b2 * (x[1]^3) + λ * x[3]
f[3] = x[3]+ω*x[4]-x[3] * (x[3]^2 + x[4]^2)
f[4] = -ω*x[3]+. x[4]-x[4] * (x[3]^2 + x[4]^2)
return f
end
and for example z0=[1., 0.001, 1, 0]
- you should use vectors instead of
z0 = [1.; 0.001]
- you should let DE compute the jacobian for you
- you should use a stiff ODE solver like
Rodas4P2
- for small ODE, you should use a direct solver
ls = DefaultLS()
- you should use a jacobian for the periodic orbits tailored for ODE. Hence, add the following options to continuation
jacobianPO = :autodiffDense,
callbackN = BK.cbMaxNorm(100.),
from bifurcationkit.jl.
Dear @rveltz ,
Many thanks for pointing out all these problems. I have just recently dived into the Julia universe, and a lot of this is new to me. Any point is deeply appreciated. Furthermore, I am impressed by this module. I really hope to use it in my studies, which are focused on free/forced vibrations of mechanical structures with large displacements and deformations. Perhaps one day I could be capable enough to contribute somehow.
Cheers!
from bifurcationkit.jl.
should be fixed on master
from bifurcationkit.jl.
Related Issues (20)
- Continuation of LPC from Bautin/GH point HOT 2
- Out of memory for 1D periodic non-linear system HOT 8
- Plotting with Makie HOT 3
- Getting print statements despite `verbosity=0` HOT 6
- BifurcationKit not working with ArrayPartitions ? HOT 1
- getperiod vs. getPeriod HOT 3
- _isinplace causes incorrect dispatch when both out of place and in place versions of a function is defined HOT 9
- Add `BifurcationProblem` dispatch for ModelingToolkit problems (possibly in an extension?) HOT 8
- BifurcationKit incompatible with latest ModelingToolkit version HOT 8
- record from solution title looks werid in docs HOT 2
- Incorrect calls to `EventSpecialPoint` HOT 2
- Problems with event detection HOT 5
- using `ContinuationPar` directly as input to `bifurcationdiagram` (rather than `(args...) -> opts_br`) HOT 2
- Custom BorderedLinearSolver interface requirements for PALC are missing HOT 4
- Generated function error HOT 2
- record_from_solution shoud have acces to the entire set of parameters HOT 2
- BK question HOT 8
- Problem with plotting a periodic orbit HOT 2
- (potentially long-term) Complete integration with ModelingTookit's nonlinear problem HOT 9
- Issue with passing the `record_from_solution` to the mimally augmented codim 2 continuation HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from bifurcationkit.jl.