In This Package:

Power.h

Go to the documentation of this file.

// ============================================================================
#ifndef LHCBMATH_POWER_H 
#define LHCBMATH_POWER_H 1
// ============================================================================
// Include files
// ============================================================================
// STD & STL 
// ============================================================================
#include <cmath>
// ============================================================================
// ============================================================================
namespace Gaudi 
{  
  namespace Math 
  {
    // ========================================================================
    template<typename TYPE>
    inline TYPE pow ( TYPE __x , unsigned long __n )
    {
      TYPE __y = __n % 2 ? __x : 1;
      
      while ( __n >>= 1 )
      {
        __x = __x * __x;
        if ( __n % 2) { __y = __y * __x; }
      } 
      
      return __y ;
    }
    // ========================================================================
  } // end of namespace Math
} // end of namespace Gaudi
// ============================================================================
// The END 
// ============================================================================
#endif // LHCBMATH_POWER_H
// ============================================================================

---