You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
For many practical pharmacometric models the integrate_ode function is pretty limiting as it only solves the ODE system at given times but cannot modify the system at those times. A more generic function would require events, which pair an event time and some used-defined function that modifies the state of the ODE at that time, for example
This new function then concatenates the vector of observation times and the
event times, removes duplicates, and orders the unique times. The solver is
run to each time, with the resulting state saved if an observation if necessary
and then input to the event functor if necessary.
Thoughts? This would require introducing a type with a method which would
be new to the language. Maybe just having an add_event function with an event_table and function as arguments?
The text was updated successfully, but these errors were encountered:
For many practical pharmacometric models the
integrate_ode
function is pretty limiting as it only solves the ODE system at given times but cannot modify the system at those times. A more generic function would require events, which pair an event time and some used-defined function that modifies the state of the ODE at that time, for exampleThese events would then be collected into a dynamically-sized vector of events,
which would be passed into a new integrate function,
This new function then concatenates the vector of observation times and the
event times, removes duplicates, and orders the unique times. The solver is
run to each time, with the resulting state saved if an observation if necessary
and then input to the event functor if necessary.
Thoughts? This would require introducing a type with a method which would
be new to the language. Maybe just having an
add_event
function with anevent_table
and function as arguments?The text was updated successfully, but these errors were encountered: