time_index Function

private function time_index(self, t, index) result(status)

Find the interval of the time list containing t, by bisection (times increasing or decreasing). Returns a TIME_ code, or the number of iterations (> 0) when found, with times(index) <= t < times(index+1) for increasing times (>= and > for decreasing). See the TIME_ constants for index in the other cases.

Type Bound

Trajectory

Arguments

Type IntentOptional Attributes Name
class(Trajectory), intent(in) :: self
real(kind=wp), intent(in) :: t

the time
integer, intent(out) :: index

the interval's first index (see the TIME_* constants)

Return Value integer


Calls

proc~~time_index~~CallsGraph proc~time_index Trajectory%time_index proc~time_distance Trajectory%time_distance proc~time_index->proc~time_distance

Called by

proc~~time_index~~CalledByGraph proc~time_index Trajectory%time_index proc~point_at Trajectory%point_at proc~point_at->proc~time_index proc~state_at TrajectoryFollower%state_at proc~state_at->proc~time_index proc~followed_pose followed_pose proc~followed_pose->proc~state_at proc~world_matrix_at ReferenceFrame%world_matrix_at proc~followed_pose->proc~world_matrix_at proc~trajectory_path_position TrajectoryPath%trajectory_path_position proc~trajectory_path_position->proc~point_at proc~update~2 TrajectoryFollower%update proc~update~2->proc~state_at proc~local_matrix_at ReferenceFrame%local_matrix_at proc~local_matrix_at->proc~followed_pose proc~update ReferenceFrame%update proc~update->proc~followed_pose proc~update->proc~update~2 proc~drawable_update DrawableTrajectory%drawable_update proc~drawable_update->proc~update proc~fm_update FrameManager%fm_update proc~fm_update->proc~update proc~sphere_update Sphere%sphere_update proc~sphere_update->proc~update proc~terrain_update TerrainBody%terrain_update proc~terrain_update->proc~update proc~world_matrix_at->proc~local_matrix_at proc~world_matrix_at->proc~world_matrix_at proc~data_to_shown data_to_shown proc~data_to_shown->proc~world_matrix_at proc~to_shown to_shown proc~to_shown->proc~data_to_shown proc~art_point_between art_point_between proc~art_point_between->proc~to_shown proc~to_shown_at to_shown_at proc~to_shown_at->proc~to_shown proc~art_point art_point proc~art_point->proc~to_shown_at proc~art_points art_points proc~art_points->proc~to_shown_at proc~intermediates_by_time intermediates_by_time proc~intermediates_by_time->proc~art_point_between proc~marker_point marker_point proc~marker_point->proc~to_shown_at proc~marker_point->proc~art_point proc~update_ribbon_tail RibbonArtist%update_ribbon_tail proc~update_ribbon_tail->proc~art_point_between proc~update_ribbon_tail->proc~art_point proc~update_tail CurveArtist%update_tail proc~update_tail->proc~art_point_between proc~update_tail->proc~art_point proc~write_segments write_segments proc~write_segments->proc~to_shown_at proc~add_intermediates add_intermediates proc~add_intermediates->proc~intermediates_by_time proc~intermediates_by_data intermediates_by_data proc~add_intermediates->proc~intermediates_by_data proc~intermediates_by_distance intermediates_by_distance proc~add_intermediates->proc~intermediates_by_distance proc~curve_pick curve_pick proc~curve_pick->proc~art_point proc~intermediates_by_data->proc~art_point proc~intermediates_by_distance->proc~art_point proc~marker_update MarkerArtist%marker_update proc~marker_update->proc~marker_point proc~marker_update->proc~add_intermediates proc~update_bounds update_bounds proc~marker_update->proc~update_bounds proc~segment_update SegmentArtist%segment_update proc~segment_update->proc~write_segments proc~update_bounds->proc~art_points proc~write_curve CurveArtist%write_curve proc~write_curve->proc~art_points proc~write_ribbon RibbonArtist%write_ribbon proc~write_ribbon->proc~art_points proc~drawable_pick DrawableTrajectory%drawable_pick proc~drawable_pick->proc~curve_pick

Source Code

    integer function time_index(self, t, index) result(status)
        class(Trajectory), intent(in) :: self
        real(wp), intent(in) :: t                       !! the time
        integer, intent(out) :: index                   !! the interval's first index (see the TIME_* constants)
        integer :: low, high, mid, iter, n
        real(wp) :: dir, dist
        n = self%n_times
        if (n == 0) then
            index = 0
            status = TIME_OUT_OF_RANGE
            return
        end if
        dir = merge(1.0_wp, -1.0_wp, self%t(1) <= self%t(n))
        dist = self%time_distance(t)
        if (dist > 0.0_wp) then
            index = merge(0, n + 1, dir * t < dir * self%t(1))
            status = TIME_OUT_OF_RANGE
            return
        else if (dist == 0.0_wp) then
            index = merge(1, n, t == self%t(1))
            status = TIME_AT_BOUND
            return
        end if
        low = 1
        high = n
        do iter = 1, 64
            mid = low + (high - low) / 2
            if (dir * t < dir * self%t(mid)) then
                high = mid
            else
                low = mid
            end if
            if (high - low <= 1) then
                index = low
                status = iter
                return
            end if
        end do
        index = low
        status = TIME_NOT_FOUND
    end function time_index