|| ||           || Doxygen's  

---

HxUtilityMath.h

/*
 *  Copyright (c) 1998, University of Amsterdam, The Netherlands.
 *  All rights reserved.
 *
 *
 *  Author(s):
 *  Silvia D. Olabarriaga  (silvia@wins.uva.nl)
 */

#ifndef HxUtilityMath_h
#define HxUtilityMath_h

#include <cmath>

#define EPS 1.0e-10

inline int
roundd(double a)
{
     return ( a >= 0.0 ? ( (int) floor(a) + 0.5 ) :
        ( (int) ceil(a) - 0.5 ) );
}

//
// square two numbers
//

inline 
int square(int x)
{
   return x*x;
}

inline 
double square(double x)
{
   return x*x;
}

//
// swap two numbers
//

inline void 
swap(int& x, int& y)
{
   int tmp = x;
   x = y;
   y = tmp;
}

inline void 
swap(double& x, double& y)
{
   double tmp = x;
   x = y;
   y = tmp;
}

//
// smaller of two numbers
//

inline int 
minimum(int x, int y)
{
   if (x < y)
      return x;
   else
      return y;
}

inline double 
minimum(double x, double y)
{
   if (x < y)
      return x;
   else
      return y;
}

//
// larger of two numbers
//

inline int 
maximum(int x, int y)
{
   if (x > y)
      return x;
   else
      return y;
}

inline double 
maximum(double x, double y)
{
   if (x > y)
      return x;
   else
      return y;
}

//
// absolute value
//

inline int 
absolute(int x)
{
   if (x < 0)
      x = -x;
   return x;
}

inline double 
absolute(double x)
{
   if (x < 0.0)
      x = -x;
   return x;
}

//
// check if number in interval 
//

inline int 
inRange(double x, double min, double max)
{
   return (x >= min && x <= max);
}

//
// compare doubles with precision
//

inline int 
compare(double x, double y)
{
   if (absolute(x - y) < EPS)
      return 1;
   else
      return 0;
}

//
// clip number to interval
//

inline int 
clip(int val, int min, int max)
{
   if(val < min)
      val = min;
   if(val > max)
      val = max;
   return val;
}

inline double 
clip(double val, double min, double max)
{
   if(val < min)
      val = min;
   if(val > max)
      val = max;
   return val;
}

//
// sign of a number
//

inline int 
sign(double x)
{
   if (x < 0.0)
    return -1;
   else if (x == 0.0)
        return 0;
   else return 1;
}

#endif 

---