math/math2.h

Go to the documentation of this file.

//-------------------------------------------------------------------
// Math2.h
//-------------------------------------------------------------------
#ifndef MATH2_H_
#define MATH2_H_ 1

#include "debug/debug.h"
#include <math.h>

//-definitions-------------------------------------------------------
#ifndef M_PI
#define M_PI    3.14159265359
#endif

#define ROUNDING_ERROR  (1e-5)
#define IsZero(x)       (fabs(x)<ROUNDING_ERROR)

//-function-prototypes-----------------------------------------------
template< class T >    bool Between( const T value, const T left, const T right );
double  Clamp( const double value, const double min, const double max );
double  CosDeg( const double value );
template< class T >    T Max( const T a, const T b );
double  Max( const double a, const double b, const double c );
double  Max( const double a, const double b, const double c, const double d );
double  Max( const double a, const double b, const double c, const double d, const double e );
template< class T >    T Min( const T a, const T b );
double  Min( const double a, const double b, const double c );
double  Min( const double a, const double b, const double c, const double d );
double  MinMag( const double a, const double b );
double  MinMag( const double a, const double b, const double c );
template< class T >    T OddOneOut( const T a, const T b, const T c );
int     Pow( const int a, const int b );
double  Rad2Deg( const double radians );
double  Deg2Rad( const double degrees );
double  RandNorm( const double mean, const double standard_devation );
double  SinDeg( const double value );
double  Sqr( const double a );
int     TwoToThePower( const int power );


template< class T >    
bool Between( const T value, const T left, const T right )
{
    if( value < math::Min( left, right ) ) return false;
    if( value > math::Max( left, right ) ) return false;
    return true;
}

namespace math
{
    template< class T >    
    inline T Max( const T a, const T b )
    {
        if( a > b ) return a;
        return b;
    }

    template< class T >    
    inline T Min( const T a, const T b )
    {
        if( a < b ) return a;
        return b;
    }
}

template< class T >
T OddOneOut( const T a, const T b, const T c )
{
    if( a == b )
    {
        return c;
    }

    if( b == c )
    {
        return a;
    }

    if( c == a )
    {
        return b;
    }
    IJG_Assert( false );
    return a;
}



#endif

---