Stopping simulation early on event
8 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
Abed Alnaif
il 25 Nov 2019
Commentato: Abed Alnaif
il 3 Dic 2019
Hello,
Is there a good way to stop a simulation early if an event occurs; i.e., stop the simulation as soon as the event is triggered? Also, I'd like to do this in a way that causes the objective function value to be infinite when using 'sbiofit'.
Thank you,
Abed
0 Commenti
Risposta accettata
Arthur Goldsipe
il 26 Nov 2019
Hi Abed,
Interesting question. The only way an event can stop a simulation early is if the event errors. sbiofit would handle such a siutation by returning Inf for the objective function value. Since you mentioned fitting, I generally recommend setting MaximumWallClock during fitting to ensure that things don't get stuck if a fit encounters a bad parameter value that results in a really long running simulation.
Ok, back to a simulation that errors. The only tricky thing with that is creating an event that is compatible with SimBiology's acceleration (that is, with the subset of the MATLAB language that can be used for code generation). Here's an example of modifying the Lotka-Volterra model to end early due to an event:
First, save the following function to myfunction.m:
function time = myfunction(time)
if time >= 0
error("Ending simulation at time " + time);
end
Then, execute the following:
sbioloadproject('lotka');
threshold = m1.addparameter('threshold',1e6);
m1.addevent('y1 > threshold', 'z = myfunction(time)');
sbioaccelerate(m1); % Ensure this model can be accelerated
simdata = sbiosimulate(m1); % Runs to completion
threshold.Value = 1e3;
simdata = sbiosimulate(m1); % Simulation errors
Più risposte (0)
Community
Più risposte nel SimBiology Community
Vedere anche
Categorie
Scopri di più su Scan Parameter Ranges in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!