math/math2.cpp

Go to the documentation of this file.

//-------------------------------------------------------------------
// Math2.cpp
//-------------------------------------------------------------------
#include "math2.h"
#include <math.h>
#include <stdlib.h>

//-------------------------------------------------------------------
double Clamp( const double value, const double min, const double max )
{
    if( value < min ) return min;
    if( value > max ) return max;
    return value;
}

//-------------------------------------------------------------------
double CosDeg( const double value )
{
        return( cos( value / 180 * M_PI ) );
};
//-------------------------------------------------------------------
double SinDeg( const double value ) 
{
        return( sin( value / 180 * M_PI ) );
}
//-------------------------------------------------------------------
double  Rad2Deg( const double radians )
{
        return radians / M_PI * 180.0 ;
};
//-------------------------------------------------------------------
double  Deg2Rad( const double degrees )
{
        return (degrees / 180.0) * M_PI;
};
//-------------------------------------------------------------------
double Max( const double a, const double b )
{
        //IMPROVE: make these templates
        if( a > b )
        {
                return a ;
        }
        return b ;
};
//-------------------------------------------------------------------
double Min( const double a, const double b )
{
        //IMPROVE: make these templates
        if( a < b )
        {
                return a ;
        }
        return b ; 
};
//-------------------------------------------------------------------
double Max( const double a, const double b, const double c )
{
        return( Max( a, Max( b, c ) ) ) ;
};
//-------------------------------------------------------------------
double  Max( const double a, const double b, const double c, const double d )
{
        return Max( Max( a, b ), Max( c, d ) );
}
//-------------------------------------------------------------------
double  Max( const double a, const double b, const double c, const double d, const double e )
{
        return Max( Max( a, b, c ), Max( d, e ) );
}
//-------------------------------------------------------------------
double Min( const double a, const double b, const double c )
{
        return ( Min( a, Min( b, c ) ) ) ;
};
//-------------------------------------------------------------------
double  Min( const double a, const double b, const double c, const double d )
{
        return Min( Min( a, b ), Min( c, d ) );
}
//-------------------------------------------------------------------
int     Min( const int a, const int b )
{
    if( a < b )
    {
        return a;
    }
    else
    {
        return b;
    }
}
//-------------------------------------------------------------------
unsigned int Min( const unsigned int a, const unsigned int b )
{
    if( a < b )
    {
        return a;
    }
    else
    {
        return b;
    }
}
//-------------------------------------------------------------------
double MinMag( const double a, const double b )
{
        //IMPROVE: make these templates
        if( fabs(a) < fabs(b) )
        {
                return a ;
        }
        return b ; 
}
//-------------------------------------------------------------------
double MinMag( const double a, const double b, const double c )
{
        return ( MinMag( a, MinMag( b, c ) ) );
}
//-------------------------------------------------------------------
int Pow( const int a, const int b )
{
    IJG_Assert( b >= 0 );
    int returnMe = 1;
    int i;
    for( i = 0; i < b; ++i )
    {
        returnMe *= a;
    }
    return returnMe;
}
//-------------------------------------------------------------------
double RandNorm( const double mean, const double standard_devation )
{
        // Normal Distribution, based on the crude model in the WinNT 3.5 manual:
        //   Normal = Mean + (-7 + Sum(14 Random Numbers [0..1])) * Standard Deviation

        const int numbers = 7;

        double sum = 0;

        for (int i = 0; i < (2*numbers); i++ )
        {
                sum += double( rand() ) / RAND_MAX ; //this is a random number between [0..1]
        }

        return ( mean + (-numbers + sum )*standard_devation );
}
//-------------------------------------------------------------------
double Sqr( const double a ) 
{
        return a*a ;
};
//-------------------------------------------------------------------
int TwoToThePower( const int power )
{
    IJG_Assert( power >= 0 );
    IJG_Assert( power <= 31 );
    int returnMe = 1;
    returnMe <<= power;
    return returnMe;
}
//-------------------------------------------------------------------

---