Measures.hh

Go to the documentation of this file.

#ifndef __DIADES__UTILS__MEASURES__HH
#define __DIADES__UTILS__MEASURES__HH

#include<set>

using namespace std;

namespace Diades
{
  namespace Utils
  {

        template<typename InputIterator1, typename InputIterator2>
        int sequenceHamming(InputIterator1 first1, InputIterator1 last1,
                                InputIterator2 first2, InputIterator2 last2)
        {
          int result = 0;
          while((first1 != last1)
                        &&
                        (first2 != last2))
                {
                  if(*first1 != *first2)
                        {
                          ++result;
                        }
                  ++first1;
                  ++first2;
                }
          if((first1 != last1) || ((first2 != last2)))
                {
                  result = -1;
                }
          return result;
        }

        
        /*******************************************************************************/

        template<typename InputIterator1, typename InputIterator2>
        int orderedHamming(InputIterator1 first1, InputIterator1 last1,
                                InputIterator2 first2, InputIterator2 last2)
        {
          int result = 0;
          
          while((first1 != last1)
                        &&
                        (first2 != last2))
                {
                  if(*first1==*first2)
                        {
                          ++first1;
                          ++first2;
                        }
                  else
                        {
                          ++result;
                          if(*first1 < *first2)
                                {
                                  ++first1;
                                }
                          else
                                {
                                  ++first2;
                                }
                        }
                }
          while(first1 != last1)
                {
                  ++first1;
                  ++result;
                }
          while(first2 != last2)
                {
                  ++first2;
                  ++result;
                }
          return result;
        }
        
        /*******************************************************************************/
        

        template<typename InputIterator1, typename InputIterator2>
        int intersectionSize(InputIterator1 first1, InputIterator1 last1,
                                                 InputIterator2 first2, InputIterator2 last2)
        {
          int result = 0;
          while((first1 != last1)
                        &&
                        (first2 != last2))
                {
                  if(*first1==*first2)
                        {
                          ++first1;
                          ++first2;
                          ++result;
                        }
                  else
                        {
                          if(*first1 < *first2)
                                {
                                  ++first1;
                                }
                          else
                                {
                                  ++first2;
                                }
                        }
                }
          return result;
        }


        /*******************************************************************************/

        template<typename InputIterator1, typename InputIterator2>
        int differenceSize(InputIterator1 first1, InputIterator1 last1,
                                                 InputIterator2 first2, InputIterator2 last2)
        {
          int result = 0;
          while((first1 != last1)
                        &&
                        (first2 != last2))
                {
                  if(*first1==*first2)
                        {
                          ++first1;
                          ++first2;
                        }
                  else
                        {
                          if(*first1 < *first2)
                                {
                                  ++first1;
                                  ++result;
                                }
                          else
                                {
                                  ++first2;
                                }
                        }
                }
          while(first1 != last1)
                {
                  ++first1;
                  ++result;
                }
          return result;
        }       
  };

  
};

#endif