libstdc++-v3/include/c_std/cmath.tcc

   1 // -*- C++ -*- C math library.
   2
   3 // Copyright (C) 2000, 2003, 2004 Free Software Foundation, Inc.
   4 //
   5 // This file is part of the GNU ISO C++ Library.  This library is free
   6 // software; you can redistribute it and/or modify it under the
   7 // terms of the GNU General Public License as published by the
   8 // Free Software Foundation; either version 2, or (at your option)
   9 // any later version.
  10
  11 // This library is distributed in the hope that it will be useful,
  12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  14 // GNU General Public License for more details.
  15
  16 // You should have received a copy of the GNU General Public License along
  17 // with this library; see the file COPYING.  If not, write to the Free
  18 // Software Foundation, 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301,
  19 // USA.
  20
  21 // As a special exception, you may use this file as part of a free software
  22 // library without restriction.  Specifically, if other files instantiate
  23 // templates or use macros or inline functions from this file, or you compile
  24 // this file and link it with other files to produce an executable, this
  25 // file does not by itself cause the resulting executable to be covered by
  26 // the GNU General Public License.  This exception does not however
  27 // invalidate any other reasons why the executable file might be covered by
  28 // the GNU General Public License.
  29
  30 // This file was written by Gabriel Dos Reis <gdr@codesourcery.com>
  31
  32 #ifndef _GLIBCXX_CMATH_TCC
  33 #define _GLIBCXX_CMATH_TCC 1
  34
  35 namespace std
  36 {
  37   template<typename _Tp>
  38     inline _Tp
  39     __cmath_power(_Tp __x, unsigned int __n)
  40     {
  41       _Tp __y = __n % 2 ? __x : 1;
  42
  43       while (__n >>= 1)
  44         {
  45           __x = __x * __x;
  46           if (__n % 2)
  47             __y = __y * __x;
  48         }
  49
  50       return __y;
  51     }
  52 }
  53
  54 #endif