Threshold.hpp

//  Threshold.hpp
//  
//  Author : Jeromy Tompkins
//  Date   : 8/14/2013

#ifndef THRESHOLD_H
#define THRESHOLD_H

#include "Trace.hpp"
#include "AlgoIterator.hpp"

namespace TrAnal {



    template<class T>
        TrIterator<T> Threshold(const TrIterator<T>& begin,
                const TrIterator<T>& end, 
                double thresh)
        {
            TrIterator<T> it=begin;
            while ( it<end && *it<thresh ) ++it;              

            return it;
        }


    template<class T>
        TrIterator<T> Threshold(const TrRange<T>& range, double thresh) { return Threshold<T>(range.begin(), range.end(), thresh);}


    template<class T>
        TrIterator<T> Threshold(AlgoIterator<T>& begin,
                const TrIterator<T>& end, 
                double thresh)
        {
            // Default return is end
            TrIterator<T> ret_it = end;

            // change the name of begin
            AlgoIterator<T>& it = begin;

            while ( it<end && *it<thresh ) {
                ++it;              
            }

            // Store result if threshold exceeded
            if ( *it >= thresh) ret_it = it.min_extent(); 

            return ret_it;
        }

} // end namespace
#endif