I have system Lienard with fixed parametres. I use alpha 0.20 and w = 0.7315 with t = 500000, when using default algorithm result is uncorrect. I use Feagin12 with abstol 20 and have problem with Instability, if i take abstol 25 i have another problem with Larger maxiters is needed. How can i fixed this? I have article with correctl result, but i need more t for my task with extreme events. I tryed anothers algorithms and this problems is saved or i hads uncorrect results.
When i use t = 100000 or 200000 i don’t have this problems, but i need more t
Have you tried increasing
In article use RK4 with step 0.01 but it is give uncorrect result in Julia. I also have result in python with method RK45 and rtol = 13, atol = 13, with it i have correct result
But what does this mean though? Did you use the RK4 solver in Julia also with non-adaptive time stepping, and set a fixed time step to 0.01 and it didn’t work? And you also used the same initial conditions? And what does “didn’t work” mean? Can you be as precise as possible?
I find this hard to believe, because the DiffEq solvers are very well tested. So it simply sounds unlikely that an RK4 implementation in another code would give “correct” results, and the one in Julia wouldn’t. Actually, if anything, I’d trust the Julia version more
I have same fixed initial condition in all research. I need get this results, that is extreme events (events > red line, dashed cyan line) but i have different result.
RK4 with fixed time step 0.01 and with t 500000 have warning “Larger maxiters is needed”
This seems completely correct to me though? If first image is paper and second image is Julia, this is completely correct.
You can’t expect to have exactly the same time evolution, even if you use exactly the same initial condition and solver algorithm. You have an extremely chaotic system in your hands. Even the fact that you simulate this system in a different computer will lead to different time evolution.
What you should be looking at instead is whether you get the same dynamical characteristics. For example, are the Lyapunov exponents the same in both cases? The fractal dimension of the chaotic attractors? If the system is multistable, do both cases show the same number of attractors?
The resolution of this warning is trivial: increase
maxiters. So we can stop discussing about this case for now and focus on whether any substantial differences exist.
This two images from articles. This image is Julia
Spectre lyapunov exponents the same in both case, attractor one
Increase maxiters before 10 000 000 helped.Thank you. Before this i increae below this value.
Did you set
adaptive=false? If you don’t do
adaptive=false then it is using adaptivity as documented.
RK4 - The canonical Runge-Kutta Order 4 method. Uses a defect control for adaptive stepping using maximum error over the whole interval.
No, don’t use. I now tried with him. Thank you.