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
 RCS:           $Id$
________________________________________________________________________
 
-*/
 
#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 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);
mGlobal(Basic) int NrSignificantDecimals(double);
 
 
} // namespace Math