mingw/lib/gcc/mingw32/4.3.3/include/c++/bits/cmath.tcc

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