Activity for Drew LaMar
-
Drew LaMar posted a comment on discussion Help
Hi, Bruno. Question: why is your free parameter for the limit cycles alpha? I changed it to mu, the same as the equilibrium point curve, and got the attached diagram. Does this look close to what you are expecting? Also, if you add PCargs.StopAtPoints = 'B' to your equilibrium point curve arguments, the curve will stop at the predetermined boundary points (plus and minus 1).
-
Hi, Bruno. I think it's the PCargs.LocBifPoints = ['all'] line. I changed it to PCargs.LocBifPoints = 'all' and it worked.
-
Hi, Ahmet. The best thing to do would be to use the 'initpoint' argument of your input dictionary and feed in a good initial guess to an equilibrium point. It looks like, based on your output, that you are at a limit cycle and not an equilibrium point. If your system is bistable, you can change initial conditions to see if you can't find an equilibrium point. You might have to change parameters to get this to occur. In any case, once you have computed a solution that settles to an equilibrium point,...
-
Actually, just realized this is not PyCont related, so it's probably best if another member of the team chases after this one. Also just noticed that you are creating things symbolically, which I don't think VFGEN deals with. Anyone else from the PyDSTool team wanna try their hand at this?
-
Hi, Ahmet. I'll take a look at your code, but in the meantime, you might want to try VFGEN by Weckesser. It's a markup language for differential equations models. What's nice about it is you only need to write up the model once, and then VFGEN can create XPPAUT and PyDSTOOL model files from the template. Pretty nice. I'll report back if I find where your model is having trouble.
-
Yes, my apologies for being absent on this one - other job commitments kept me away...
-
You can specify the initdirec argument of the continuer class, as a dictionary.
-
Hi, Maurizio. I think the non-reproducibility is due to the random default initialization...
-
Hi, Maurizio. I've seen this a few times - I think it's a bug. Any chance you could...
-
Hi, Emil. Not sure I understand what is going on now. What do you mean by the free...
-
If you are able to continue and find bifurcation points, then the number of points...
-
Hi, Danny. It's been awhile since I've taken a look at the code, but I'm pretty sure...
-
Hi, Maurizio. The numerical algorithm to detect bifurcation points is different than...
-
Got it. Forgot that I can check the build.log. :/ Yep, sinf was already defined in...
-
Hi, everyone. Sorry I’ve been out of the loop lately. I have a specific model that...
-
Hmmm, you are correct that it is a dictionary, although I'm not sure what's going...
-
Yes, essentially - some function of the state variables and parameters that remains...
-
This is a common issue with continuation. In the past, I've gotten this problem when...
-
No problem, glad it works. I think Rob would be best at answering your second question...
-
Awesome, thanks, FLorian. I'll pull it soon. Drew
-
Sounds good. Go ahead and post! On Aug 27, 2014, at 10:47 AM, Maurizio De Pitta'...
-
Sounds good. Please post as an enhancement on github. Drew
-
I don't think so. The most natural thing to do would be to extend the functionality...
-
Here is what I got: >>> PyCont.gensys.set(pdomain={'p1': [0.0, 2.0]}) Traceback (most...
-
Oh, check out PyCont_Logistic.py example to see how pdomain is set in the beginning....
-
Oh, I think I understand what's going on. The pdomain argument is given to the ODE...
-
Hi, Maurizio. Let me look into the PD issue. It sounds like you want to simply double...
-
Hi, Florian. I am planning on diving into the AUTO code soon to implement tracking...
-
Hi, Maurizio. Does the StopAtPoint set to the boundary work when you specify parameter...
-
Hi, Maurizio. First of all, PyCont can't continue NS points yet for limit cycle curves,...
-
Hi, Maurizio. First of all, PyCont can't continue NS points yet for limit cycle curves,...
-
Your "flow" of questions have been invaluable! So thank you for that. I wasn't implying...
-
Hi, Maurizio. I know EXACTLY where this error is coming from. Essentially, to calculate...
-
Hi, Maurizio. I know EXACTLY where this error is coming from. Essentially, to calculate...
-
Hi again, Maurizio. Another good example of fine tuning of plots is in examples/PyCont_ABReaction.py....
-
I'll put in the bug issue. :/ Lots of bugs lately, but I REALLY appreciate you all...
-
Yeah, it currently is not possible to start again from an "MX" point for an LC-C...
-
Oops, sorry, here's the axis handle as well: PyCont.plot.fig1.axes1.axes Drew
-
Hi, Maurizio. You can get the figure handle as follows: PyCont.plot.fig1.fig Dre...
-
Yes, but it would be nice though to be able to just delete the portion of the curve...
-
Hi, Maurizio. I think the update command is the correct way to go about doing what...
-
Hi, Maurizio. I think the current code only can extend curves from the endpoints....
-
That, my friend, would be awesome. :) You are more than welcome to join the party....
-
Hi, Greg. You are correct: that is what the code is supposed to do. I will put in...
-
Hi, Greg. You are correct: that is what the code is supposed to do. I will put in...
-
Hi, Greg. I found the problem: You are naming your curves the same as bifurcation...
-
Hi, Greg. Strange, I was unable to reproduce your problem. I tried all kinds of variations...
-
I don't think we ever had that functionality, actually. This just highlights the...
-
Hi, Greg. I think I understand what is going on. When gam > 1, there are no LP points...
-
Hi, Greg. Any chance you could share your full code? If it's research code you are...
-
Hi Greg. The points that you are seeing beyond the BT point (gam < 0.5) are neutral...
-
Hi Greg. Yeah, just to be sure I took a look at the eigenvalues and sure enough,...
-
Hi, Greg. I've attached the 2-parameter bifurcation diagram I was able to reproduce...
-
Hi, Greg. Here's the 2-parameter bifurcation diagram I was able to reproduce from...
-
Hi, everyone. Sorry I'm late to the conversation. I agree that post-calculation is...
-
Taking a look at your system, you have a lot of symmetry going on. In particular,...
-
Hi, Henrik. The MoorePenrose corrector is much more stable. I changed the corrector...
-
Hi Matt. The negative "D" value does not indicate a break in the continuation, as...
1