00001 // Numeric extensions -*- C++ -*- 00002 00003 // Copyright (C) 2002, 2004 Free Software Foundation, Inc. 00004 // 00005 // This file is part of the GNU ISO C++ Library. This library is free 00006 // software; you can redistribute it and/or modify it under the 00007 // terms of the GNU General Public License as published by the 00008 // Free Software Foundation; either version 2, or (at your option) 00009 // any later version. 00010 00011 // This library is distributed in the hope that it will be useful, 00012 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00013 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00014 // GNU General Public License for more details. 00015 00016 // You should have received a copy of the GNU General Public License along 00017 // with this library; see the file COPYING. If not, write to the Free 00018 // Software Foundation, 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, 00019 // USA. 00020 00021 // As a special exception, you may use this file as part of a free software 00022 // library without restriction. Specifically, if other files instantiate 00023 // templates or use macros or inline functions from this file, or you compile 00024 // this file and link it with other files to produce an executable, this 00025 // file does not by itself cause the resulting executable to be covered by 00026 // the GNU General Public License. This exception does not however 00027 // invalidate any other reasons why the executable file might be covered by 00028 // the GNU General Public License. 00029 00030 /* 00031 * 00032 * Copyright (c) 1994 00033 * Hewlett-Packard Company 00034 * 00035 * Permission to use, copy, modify, distribute and sell this software 00036 * and its documentation for any purpose is hereby granted without fee, 00037 * provided that the above copyright notice appear in all copies and 00038 * that both that copyright notice and this permission notice appear 00039 * in supporting documentation. Hewlett-Packard Company makes no 00040 * representations about the suitability of this software for any 00041 * purpose. It is provided "as is" without express or implied warranty. 00042 * 00043 * 00044 * Copyright (c) 1996 00045 * Silicon Graphics Computer Systems, Inc. 00046 * 00047 * Permission to use, copy, modify, distribute and sell this software 00048 * and its documentation for any purpose is hereby granted without fee, 00049 * provided that the above copyright notice appear in all copies and 00050 * that both that copyright notice and this permission notice appear 00051 * in supporting documentation. Silicon Graphics makes no 00052 * representations about the suitability of this software for any 00053 * purpose. It is provided "as is" without express or implied warranty. 00054 */ 00055 00061 #ifndef _EXT_NUMERIC 00062 #define _EXT_NUMERIC 1 00063 00064 #pragma GCC system_header 00065 00066 #include <bits/concept_check.h> 00067 #include <numeric> 00068 00069 #include <ext/functional> // For identity_element 00070 00071 namespace __gnu_cxx 00072 { 00073 // Returns __x ** __n, where __n >= 0. _Note that "multiplication" 00074 // is required to be associative, but not necessarily commutative. 00075 template<typename _Tp, typename _Integer, typename _MonoidOperation> 00076 _Tp 00077 __power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op) 00078 { 00079 if (__n == 0) 00080 return identity_element(__monoid_op); 00081 else 00082 { 00083 while ((__n & 1) == 0) 00084 { 00085 __n >>= 1; 00086 __x = __monoid_op(__x, __x); 00087 } 00088 00089 _Tp __result = __x; 00090 __n >>= 1; 00091 while (__n != 0) 00092 { 00093 __x = __monoid_op(__x, __x); 00094 if ((__n & 1) != 0) 00095 __result = __monoid_op(__result, __x); 00096 __n >>= 1; 00097 } 00098 return __result; 00099 } 00100 } 00101 00102 template<typename _Tp, typename _Integer> 00103 inline _Tp 00104 __power(_Tp __x, _Integer __n) 00105 { return __power(__x, __n, std::multiplies<_Tp>()); } 00106 00112 // Alias for the internal name __power. Note that power is an extension, 00113 // not part of the C++ standard. 00114 template<typename _Tp, typename _Integer, typename _MonoidOperation> 00115 inline _Tp 00116 power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op) 00117 { return __power(__x, __n, __monoid_op); } 00118 00124 template<typename _Tp, typename _Integer> 00125 inline _Tp 00126 power(_Tp __x, _Integer __n) 00127 { return __power(__x, __n); } 00128 00134 // iota is not part of the C++ standard. It is an extension. 00135 template<typename _ForwardIter, typename _Tp> 00136 void 00137 iota(_ForwardIter __first, _ForwardIter __last, _Tp __value) 00138 { 00139 // concept requirements 00140 __glibcxx_function_requires(_Mutable_ForwardIteratorConcept<_ForwardIter>) 00141 __glibcxx_function_requires(_ConvertibleConcept<_Tp, 00142 typename std::iterator_traits<_ForwardIter>::value_type>) 00143 00144 while (__first != __last) 00145 *__first++ = __value++; 00146 } 00147 } // namespace __gnu_cxx 00148 00149 #endif 00150