math2.h

Go to the documentation of this file.

#pragma once

/*+
________________________________________________________________________

 (C) dGB Beheer B.V.; (LICENSE) http://opendtect.org/OpendTect_license.txt
 Author:        A.H. Bril
 Date:          Aug 2005
________________________________________________________________________

-*/

#include "basicmod.h"
#include "gendefs.h"
#include "odcomplex.h"

/* Functions with some extra facilities added to math.h

   Each function has a float and a double version. This may be very
   old-fashioned, but who cares.

   At the end, there is also an all-integer x to the power y.

 */


namespace Math
{

mGlobal(Basic) unsigned int SetBits( unsigned int curflags,
                                     unsigned int mask, bool yn );

mGlobal(Basic) bool AreBitsSet( unsigned int curflags,
                                unsigned int mask, bool all=true );

mGlobal(Basic) bool IsNormalNumber(float);

mGlobal(Basic) float PowerOf(float,float);

mGlobal(Basic) float ASin(float);

mGlobal(Basic) float ACos(float);

mGlobal(Basic) float Log(float);

mGlobal(Basic) float Log10(float);

mGlobal(Basic) float Sqrt(float);

mGlobal(Basic) float BesselI0(float);

mGlobal(Basic) float_complex Sqrt(const float_complex&);

mGlobal(Basic) float Exp(float);

mGlobal(Basic) float toDB(float);

mGlobal(Basic) inline unsigned int Abs( unsigned int i )    { return i; }
mGlobal(Basic) inline od_uint64 Abs( od_uint64 i )          { return i; }
mGlobal(Basic) unsigned int Abs(int i);
mGlobal(Basic) od_uint64 Abs(od_int64 i);
mGlobal(Basic) float Abs(float);
mGlobal(Basic) double Abs(double);
mGlobal(Basic) float Abs(float_complex);

mGlobal(Basic) float Floor(float);
mGlobal(Basic) double Floor(double);
mGlobal(Basic) float Ceil(float);
mGlobal(Basic) double Ceil(double);

mGlobal(Basic)  float Sin(float);
mGlobal(Basic)  double Sin(double);
mGlobal(Basic)  float Cos(float);
mGlobal(Basic)  double Cos(double);
mGlobal(Basic)  float Sinc(float);
mGlobal(Basic)  double Sinc(double);
mGlobal(Basic)  float Atan2(float y, float x);
mGlobal(Basic)  double Atan2(double y, double x);
mGlobal(Basic)  float Atan2(float_complex timeval);

inline float    toDegrees(float r)      { return (mIsUdf(r) ? r : r*mRad2DegF);}
inline double   toDegrees(double r)     { return (mIsUdf(r) ? r : r*mRad2DegD);}
inline float    toRadians(float d)      { return (mIsUdf(d) ? d : d*mDeg2RadF);}
inline double   toRadians(double d)     { return (mIsUdf(d) ? d : d*mDeg2RadD);}

inline float    degFromNorth( float azimuth )
                {
                    const float deg = 90 - mRad2DegF * azimuth;
                    return deg < 0 ? deg + 360 : deg;
                }

mGlobal(Basic) int LCMOf(int,int); 
mGlobal(Basic) int HCFOf(int,int); 
mGlobal(Basic) bool IsNormalNumber(double);
mGlobal(Basic) double PowerOf(double,double);
mGlobal(Basic) double ASin(double);
mGlobal(Basic) double ACos(double);
mGlobal(Basic) double Log(double);
mGlobal(Basic) double Log10(double);
mGlobal(Basic) double Exp(double);
mGlobal(Basic) double Sqrt(double);
mGlobal(Basic) double toDB(double);
mGlobal(Basic) double BesselI0(double);

template <class iT,class iPOW> inline
iT IntPowerOf( iT i, iPOW p )
{
    iT ret = 1;
    while ( p )
    {
        if ( p > 0 )
            { ret *= i; p--; }
        else
            { ret /= i; p++; }
    }
    return ret;
}

mGlobal(Basic) double IntPowerOf(double,int);
mGlobal(Basic) float IntPowerOf(float,int);


} // namespace Math