utilities.hpp

// This file is part of the imaging2 class library.
//
// University of Innsbruck, Infmath Imaging, 2009.
// http://infmath.uibk.ac.at
//
// All rights reserved.


#ifndef CORE_UTILITIES
#define CORE_UTILITIES


#include <iostream>
#include <sstream>
#include <math.h>
#include <core/imaging2.hpp>

#include <boost/lexical_cast.hpp>


namespace imaging
{
  const double PI = 3.1415926;
  
  const double SQUARE_ROOT_2 = 1.414213562373095049;
  
  const double SQUARE_ROOT_3 = 1.73205080757;

  template <class num_t>
  inline num_t max(const num_t a, const num_t b) { return a > b ? a : b; }

  template <class num_t>
  inline num_t min(const num_t a, const num_t b) { return a < b ? a : b; }

  template <class num_t>
  inline num_t square(num_t a) { return a * a; }

  template <class num_t>
  inline num_t power(num_t a, size_t b)
  { 
    if(b == 0)
      return num_t(1.0);
    else
      return power(a, b - 1) * a;
  }
  
  template <class num_t>
  inline num_t abs(num_t a)
  {
    return a < 0 ? -a : a;
  }
  
  template <>
  inline float abs(float a)
  {
    return fabs(a);
  }
  
  template <>
  inline double abs(double a)
  {
    return fabs(a);
  }
  
  
  template <class num_t>
  inline num_t sign(num_t a)
  {
    if(a == num_t(0))
      return num_t(0);
      
    return a < 0 ? num_t(-1) : num_t(1);
  }
  
  
  float_t delta(size_t base, size_t t);
  
  
  float_t clockwise_difference(float_t angle_1, float_t angle_2);
  
  
  float_t counter_clockwise_difference(float_t angle_1, float_t angle_2);
}

#endif

---