00001 // Numeric extensions -*- C++ -*- 00002 00003 // Copyright (C) 2002 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, 59 Temple Place - Suite 330, Boston, MA 02111-1307, 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 00056 /** @file ext/numeric 00057 * This file is a GNU extension to the Standard C++ Library (possibly 00058 * containing extensions from the HP/SGI STL subset). You should only 00059 * include this header if you are using GCC 3 or later. 00060 */ 00061 00062 #ifndef _EXT_NUMERIC 00063 #define _EXT_NUMERIC 00064 00065 #pragma GCC system_header 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 00076 template<typename _Tp, typename _Integer, typename _MonoidOperation> 00077 _Tp 00078 __power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op) 00079 { 00080 if (__n == 0) 00081 return identity_element(__monoid_op); 00082 else { 00083 while ((__n & 1) == 0) { 00084 __n >>= 1; 00085 __x = __monoid_op(__x, __x); 00086 } 00087 00088 _Tp __result = __x; 00089 __n >>= 1; 00090 while (__n != 0) { 00091 __x = __monoid_op(__x, __x); 00092 if ((__n & 1) != 0) 00093 __result = __monoid_op(__result, __x); 00094 __n >>= 1; 00095 } 00096 return __result; 00097 } 00098 } 00099 00100 template<typename _Tp, typename _Integer> 00101 inline _Tp 00102 __power(_Tp __x, _Integer __n) 00103 { return __power(__x, __n, std::multiplies<_Tp>()); } 00104 00105 // Alias for the internal name __power. Note that power is an extension, 00106 // not part of the C++ standard. 00107 00108 /** 00109 * This is an SGI extension. 00110 * @ingroup SGIextensions 00111 * @doctodo 00112 */ 00113 template<typename _Tp, typename _Integer, typename _MonoidOperation> 00114 inline _Tp 00115 power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op) 00116 { return __power(__x, __n, __monoid_op); } 00117 00118 /** 00119 * This is an SGI extension. 00120 * @ingroup SGIextensions 00121 * @doctodo 00122 */ 00123 template<typename _Tp, typename _Integer> 00124 inline _Tp 00125 power(_Tp __x, _Integer __n) 00126 { return __power(__x, __n); } 00127 00128 // iota is not part of the C++ standard. It is an extension. 00129 00130 /** 00131 * This is an SGI extension. 00132 * @ingroup SGIextensions 00133 * @doctodo 00134 */ 00135 template<typename _ForwardIter, typename _Tp> 00136 void 00137 iota(_ForwardIter __first, _ForwardIter __last, _Tp __value) 00138 { 00139 // concept requirements 00140 __glibcpp_function_requires(_Mutable_ForwardIteratorConcept<_ForwardIter>) 00141 __glibcpp_function_requires(_ConvertibleConcept<_Tp, 00142 typename std::iterator_traits<_ForwardIter>::value_type>) 00143 00144 while (__first != __last) 00145 *__first++ = __value++; 00146 } 00147 00148 } // namespace __gnu_cxx 00149 00150 #endif /* _EXT_NUMERIC */ 00151