Subversion Repositories Kolibri OS

Rev

Go to most recent revision | Blame | Last modification | View Log | RSS feed

  1. // Numeric extensions -*- C++ -*-
  2.  
  3. // Copyright (C) 2002-2013 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 3, 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. // Under Section 7 of GPL version 3, you are granted additional
  17. // permissions described in the GCC Runtime Library Exception, version
  18. // 3.1, as published by the Free Software Foundation.
  19.  
  20. // You should have received a copy of the GNU General Public License and
  21. // a copy of the GCC Runtime Library Exception along with this program;
  22. // see the files COPYING3 and COPYING.RUNTIME respectively.  If not, see
  23. // <http://www.gnu.org/licenses/>.
  24.  
  25. /*
  26.  *
  27.  * Copyright (c) 1994
  28.  * Hewlett-Packard Company
  29.  *
  30.  * Permission to use, copy, modify, distribute and sell this software
  31.  * and its documentation for any purpose is hereby granted without fee,
  32.  * provided that the above copyright notice appear in all copies and
  33.  * that both that copyright notice and this permission notice appear
  34.  * in supporting documentation.  Hewlett-Packard Company makes no
  35.  * representations about the suitability of this software for any
  36.  * purpose.  It is provided "as is" without express or implied warranty.
  37.  *
  38.  *
  39.  * Copyright (c) 1996
  40.  * Silicon Graphics Computer Systems, Inc.
  41.  *
  42.  * Permission to use, copy, modify, distribute and sell this software
  43.  * and its documentation for any purpose is hereby granted without fee,
  44.  * provided that the above copyright notice appear in all copies and
  45.  * that both that copyright notice and this permission notice appear
  46.  * in supporting documentation.  Silicon Graphics makes no
  47.  * representations about the suitability of this software for any
  48.  * purpose.  It is provided "as is" without express or implied warranty.
  49.  */
  50.  
  51. /** @file ext/numeric
  52.  *  This file is a GNU extension to the Standard C++ Library (possibly
  53.  *  containing extensions from the HP/SGI STL subset).
  54.  */
  55.  
  56. #ifndef _EXT_NUMERIC
  57. #define _EXT_NUMERIC 1
  58.  
  59. #pragma GCC system_header
  60.  
  61. #include <bits/concept_check.h>
  62. #include <numeric>
  63.  
  64. #include <ext/functional> // For identity_element
  65.  
  66. namespace __gnu_cxx _GLIBCXX_VISIBILITY(default)
  67. {
  68. _GLIBCXX_BEGIN_NAMESPACE_VERSION
  69.  
  70.   // Returns __x ** __n, where __n >= 0.  _Note that "multiplication"
  71.   // is required to be associative, but not necessarily commutative.
  72.   template<typename _Tp, typename _Integer, typename _MonoidOperation>
  73.     _Tp
  74.     __power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op)
  75.     {
  76.       if (__n == 0)
  77.         return identity_element(__monoid_op);
  78.       else
  79.         {
  80.           while ((__n & 1) == 0)
  81.             {
  82.               __n >>= 1;
  83.               __x = __monoid_op(__x, __x);
  84.             }
  85.  
  86.           _Tp __result = __x;
  87.           __n >>= 1;
  88.           while (__n != 0)
  89.             {
  90.               __x = __monoid_op(__x, __x);
  91.               if ((__n & 1) != 0)
  92.                 __result = __monoid_op(__result, __x);
  93.               __n >>= 1;
  94.             }
  95.           return __result;
  96.         }
  97.     }
  98.  
  99.   template<typename _Tp, typename _Integer>
  100.     inline _Tp
  101.     __power(_Tp __x, _Integer __n)
  102.     { return __power(__x, __n, std::multiplies<_Tp>()); }
  103.  
  104.   /**
  105.    *  This is an SGI extension.
  106.    *  @ingroup SGIextensions
  107.    *  @doctodo
  108.   */
  109.   // Alias for the internal name __power.  Note that power is an extension,
  110.   // not part of the C++ standard.
  111.   template<typename _Tp, typename _Integer, typename _MonoidOperation>
  112.     inline _Tp
  113.     power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op)
  114.     { return __power(__x, __n, __monoid_op); }
  115.  
  116.   /**
  117.    *  This is an SGI extension.
  118.    *  @ingroup SGIextensions
  119.    *  @doctodo
  120.   */
  121.   template<typename _Tp, typename _Integer>
  122.     inline _Tp
  123.     power(_Tp __x, _Integer __n)
  124.     { return __power(__x, __n); }
  125.  
  126. #if __cplusplus >= 201103L
  127.   using std::iota;
  128. #else
  129.   /**
  130.    *  This is an SGI extension.
  131.    *  @ingroup SGIextensions
  132.    *  @doctodo
  133.   */
  134.   // iota is not part of the C++ standard.  It is an extension.
  135.   template<typename _ForwardIter, typename _Tp>
  136.     void
  137.     iota(_ForwardIter __first, _ForwardIter __last, _Tp __value)
  138.     {
  139.       // concept requirements
  140.       __glibcxx_function_requires(_Mutable_ForwardIteratorConcept<_ForwardIter>)
  141.       __glibcxx_function_requires(_ConvertibleConcept<_Tp,
  142.             typename std::iterator_traits<_ForwardIter>::value_type>)
  143.  
  144.       while (__first != __last)
  145.         *__first++ = __value++;
  146.     }
  147. #endif  // C++11
  148.  
  149. _GLIBCXX_END_NAMESPACE_VERSION
  150. } // namespace
  151.  
  152. #endif
  153.