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

   1 // -*- C++ -*- C math library.
   2
   3 // Copyright (C) 2000 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, 59 Temple Place - Suite 330, Boston, MA 02111-1307,
  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 _CPP_BITS_CMATH_TCC
  33 #define _CPP_BITS_CMATH_TCC 1
  34
  35 namespace std {
  36   export template<typename _Tp>
  37     _Tp
  38     __cmath_power(_Tp __x, unsigned int __n)
  39     {
  40       _Tp __y = __n % 2 ? __x : 1;
  41
  42       while (__n >>= 1)
  43         {
  44           __x = __x * __x;
  45           if (__n % 2)
  46             __y = __y * __x;
  47         }
  48
  49       return __y;
  50     }
  51 }
  52
  53 #endif