adevs
|
#include <adevs_bisection_event_locator.h>
Public Member Functions | |
bisection_event_locator (ode_system< X > *sys, double err_tol) | |
bool | find_events (bool *events, const double *qstart, double *qend, ode_solver< X > *solver, double &h) |
This is a state event locator that uses bisection to find an event in time. The specific algorithm used here is guaranteed to return a timestep that places the state variables on the opposite side of the event surface. That is, the sign of the function is guaranteed to change.
adevs::bisection_event_locator::bisection_event_locator | ( | ode_system< X > * | sys, |
double | err_tol | ||
) |
The locator will attempt to pinpoint events within err_tol of zero on the opposite side of the event surface. An event occurs at time t' if |z(t')| < err_tol and z(t)z(t') < 0, where t is the time at the start of the search.
References adevs::ode_system::numEvents().
bool adevs::bisection_event_locator::find_events | ( | bool * | events, |
const double * | qstart, | ||
double * | qend, | ||
ode_solver< X > * | solver, | ||
double & | h | ||
) | [virtual] |
Find the first state event in the interval [0,h] starting from state qstart. The method returns true if an event is found, setting the events flags to true if the corresponding z entry in the state_event_func above triggered the event. The value of h is overwritten with the event time, and the state of the model at that time is copied to qend.
Implements adevs::event_locator< X >.
References adevs::ode_solver::advance().