finitePRC(model,
ref_traj_period,
evname,
pertcoord,
pertsize=0.05,
settle=5,
verbose=False,
skip=1,
do_pert=<function _default_pert at 0x95002b0>,
keep_trajs=False,
stop_at_t=inf,
force_T=nan)
| source code
|
Return a Pointset with dependent variable 'D_phase', measured from 0 to 1,
where D_phase > 0 is an advance.
Pass a Generator or Model instance for model.
Pass a Trajectory or Pointset for the ref_traj_period argument.
Pass the event name in the model that indicates the periodicity.
Use skip > 1 to sub-sample the points computed along the trajectory at
the skip rate.
Use a do_pert function to do any non-standard perturbation, e.g. if there
are domain boundary conditions that need special treatment. This function
takes four or five arguments (model, ic, pertcoord, pertsize, perttime=None)
and returns the new point ic (not just ic[pertcoord]).
Use settle=0 to perform no forward integration before the time window in
which the perturbation will be applied, or a fraction < 1 to ensure an
integration past the event point (e.g. for non-cycles).
Use stop_at_t to calculate a partial PRC, from perturbation time 0 to this
value.
Use force_T to force the period to be whatever value you like.
Note: Depending on your model, there may be regions of the PRC that are
offset by a constant amount to the rest of the PRC. This is a "wart" that
needs improvement.
|