H3DMath.h

//    Copyright 2004, SenseGraphics AB
//
//    This file is part of H3D API.
//
//    H3D API is free software; you can redistribute it and/or modify
//    it under the terms of the GNU General Public License as published by
//    the Free Software Foundation; either version 2 of the License, or
//    (at your option) any later version.
//
//    H3D API is distributed in the hope that it will be useful,
//    but WITHOUT ANY WARRANTY; without even the implied warranty of
//    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//    GNU General Public License for more details.
//
//    You should have received a copy of the GNU General Public License
//    along with H3D API; if not, write to the Free Software
//    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
//
//    A commercial license is also available. Please contact us at 
//    www.sensegraphics.com for more information.
//
//
//

#ifndef __H3DMATH_H__
#define __H3DMATH_H__

#include "Exception.h"
#include <math.h>

namespace H3D {

  namespace Constants {
    static const double pi = 3.1415926535897932384626433832795;
  }

  inline bool isPowerOfTwo( unsigned int i ) {
    return (((i - 1) & i) == 0);
  }

  inline unsigned int nextPowerOfTwo( unsigned int i ) {
    i--;
    i |= i >> 16;
    i |= i >> 8;
    i |= i >> 4;
    i |= i >> 2;
    i |= i >> 1;
    i++;
    return i;
  }

  template< class F >
  inline F H3DAbs( F f ) {
    return fabs( f );
  }

  template< class F >
  inline F H3DSqrt( F f ) {
    return sqrt( f );
  }

  template< class F >
  inline F H3DExp( F f ) {
    return exp( f );
  }

  template< class F >
  inline F H3DLog( F f ) {
    return log( f );
  }

  template< class F >
  inline F H3DLog10( F f ) {
    return log10( f );
  }

  template< class F, class T >
  inline F H3DPow( F f, T t ) {
    return pow( f, t );
  }

  inline double H3DPow( int f, int t ) {
    return pow( (double)f, (double)t );
  }
  
  template< class A >
  inline A H3DMax( const A &a, const A &b  ) {
    return max( a, b );
  }
    
  template< class A >
  inline A H3DMin( const A &a, const A &b  ) {
    return min( a, b );
  }

  template< class F >
  inline F H3DCos( F d ) {
    return cos( d );
  }

  template< class F >
  inline F H3DCosh( F d ) {
    return cosh( d );
  }

  template< class F >
  inline F H3DAcos( F d ) {
    return acos( d );
  }

  template< class F >
  inline  F H3DSin( F d ) {
    return sin( d );
  }

  template< class F >
  inline F H3DSinh( F d ) {
    return sinh( d );
  }

  template< class F >
  inline F H3DAsin( F d ) {
    return asin( d );
  }

  template< class F >
  inline F H3DTan( F d ) {
    return tan( d );
  }

  template< class F >
  inline F H3DTanh( F d ) {
    return tanh( d );
  }

  template< class F >
  inline F H3DAtan( F d ) {
    return atan( d );
  }

  template< class F >
  inline F H3DAtan2( F y, F x ) {
    return atan2( y, x );
  }

  template< class F >
  inline F H3DFloor( F f ) {
    return floor( f );
  }

  template< class F >
  inline F H3DCeil( F f ) {
    return ceil( f );
  }
}

#endif

---