AngleComputer.ih

// IMPLEMENTATION of inline methods.

#include <cstdlib>



// Implementation of inline methods                                          //



inline
float 
DGtal::AngleComputer::cast( float i )
{
  while ( i < 0.0f ) i += (float) (M_PI*2.0);
  while ( i > (double)(M_PI*2.0) ) i -= (float)(M_PI*2.0);
  return i;
}

    
inline
bool
DGtal::AngleComputer::less( float i, float j )
{
  float d = j - i;
  if ( d > 0.0f )
    return d < M_PI;
  else
    return (-d) >= M_PI;
}


inline
float 
DGtal::AngleComputer::posDiff( float j, float i )
{
  return ( i <= j ) ? j - i : j + (float)(M_PI*2.0) - i;
}

inline
float 
DGtal::AngleComputer::deviation( float j, float i )
{
  return less( i, j ) ? posDiff( j, i ) : -posDiff( i, j );
}


inline
float 
DGtal::AngleComputer::min( float i, float j )
{
  return less( i, j ) ? i : j;
}

inline
float 
DGtal::AngleComputer::max( float i, float j )
{
  return less( i, j ) ? j : i;
}



inline
double 
DGtal::AngleComputer::cast( double i )
{
  while ( i < 0.0 ) i += (float)(M_PI*2.0);
  while ( i > (float)(M_PI*2.0) ) i -= (float)(M_PI*2.0);
  return i;
}

    
inline
bool
DGtal::AngleComputer::less( double i, double j )
{
  double d = j - i;
  if ( d > 0.0 )
    return d < M_PI;
  else
    return (-d) >= M_PI;
}


inline
double 
DGtal::AngleComputer::posDiff( double j, double i )
{
  return ( i <= j ) ? j - i : j + (float)(M_PI*2.0) - i;
}

inline
double 
DGtal::AngleComputer::deviation( double j, double i )
{
  return less( i, j ) ? posDiff( j, i ) : -posDiff( i, j );
}


inline
double 
DGtal::AngleComputer::min( double i, double j )
{
  return less( i, j ) ? i : j;
}

inline
double 
DGtal::AngleComputer::max( double i, double j )
{
  return less( i, j ) ? j : i;
}




//                                                                           //