]>
Commit | Line | Data |
---|---|---|
01b0d4b8 PC |
1 | // Numeric extensions -*- C++ -*- |
2 | ||
8d9254fc | 3 | // Copyright (C) 2002-2020 Free Software Foundation, Inc. |
01b0d4b8 PC |
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 | |
748086b7 | 8 | // Free Software Foundation; either version 3, or (at your option) |
01b0d4b8 PC |
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 | ||
748086b7 JJ |
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. | |
01b0d4b8 | 19 | |
748086b7 JJ |
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/>. | |
01b0d4b8 PC |
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 | ||
ffe94f83 PE |
51 | /** @file ext/numeric |
52 | * This file is a GNU extension to the Standard C++ Library (possibly | |
0aa06b18 | 53 | * containing extensions from the HP/SGI STL subset). |
ffe94f83 PE |
54 | */ |
55 | ||
01b0d4b8 | 56 | #ifndef _EXT_NUMERIC |
b2acb86f | 57 | #define _EXT_NUMERIC 1 |
01b0d4b8 PC |
58 | |
59 | #pragma GCC system_header | |
b2acb86f | 60 | |
01b0d4b8 | 61 | #include <bits/concept_check.h> |
54c1bf78 | 62 | #include <numeric> |
01b0d4b8 PC |
63 | |
64 | #include <ext/functional> // For identity_element | |
65 | ||
12ffa228 BK |
66 | namespace __gnu_cxx _GLIBCXX_VISIBILITY(default) |
67 | { | |
68 | _GLIBCXX_BEGIN_NAMESPACE_VERSION | |
3cbc7af0 | 69 | |
01b0d4b8 | 70 | // Returns __x ** __n, where __n >= 0. _Note that "multiplication" |
fa30fe72 | 71 | // is required to be associative, but not necessarily commutative. |
01b0d4b8 PC |
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); | |
d962e073 PC |
78 | else |
79 | { | |
80 | while ((__n & 1) == 0) | |
81 | { | |
82 | __n >>= 1; | |
83 | __x = __monoid_op(__x, __x); | |
84 | } | |
01b0d4b8 | 85 | |
d962e073 | 86 | _Tp __result = __x; |
01b0d4b8 | 87 | __n >>= 1; |
d962e073 PC |
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; | |
01b0d4b8 | 96 | } |
01b0d4b8 PC |
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 | ||
fc7f0a80 PE |
104 | /** |
105 | * This is an SGI extension. | |
106 | * @ingroup SGIextensions | |
107 | * @doctodo | |
108 | */ | |
b2acb86f BK |
109 | // Alias for the internal name __power. Note that power is an extension, |
110 | // not part of the C++ standard. | |
01b0d4b8 PC |
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 | ||
fc7f0a80 PE |
116 | /** |
117 | * This is an SGI extension. | |
118 | * @ingroup SGIextensions | |
119 | * @doctodo | |
120 | */ | |
01b0d4b8 PC |
121 | template<typename _Tp, typename _Integer> |
122 | inline _Tp | |
123 | power(_Tp __x, _Integer __n) | |
124 | { return __power(__x, __n); } | |
125 | ||
734f5023 | 126 | #if __cplusplus >= 201103L |
7c16382a JY |
127 | using std::iota; |
128 | #else | |
fc7f0a80 PE |
129 | /** |
130 | * This is an SGI extension. | |
131 | * @ingroup SGIextensions | |
132 | * @doctodo | |
133 | */ | |
3263fb9c | 134 | // iota is not part of the standard until C++11. It is an extension. |
01b0d4b8 | 135 | template<typename _ForwardIter, typename _Tp> |
fa30fe72 | 136 | void |
01b0d4b8 PC |
137 | iota(_ForwardIter __first, _ForwardIter __last, _Tp __value) |
138 | { | |
139 | // concept requirements | |
3d7c150e BK |
140 | __glibcxx_function_requires(_Mutable_ForwardIteratorConcept<_ForwardIter>) |
141 | __glibcxx_function_requires(_ConvertibleConcept<_Tp, | |
01b0d4b8 PC |
142 | typename std::iterator_traits<_ForwardIter>::value_type>) |
143 | ||
144 | while (__first != __last) | |
145 | *__first++ = __value++; | |
146 | } | |
734f5023 | 147 | #endif // C++11 |
3cbc7af0 | 148 | |
12ffa228 BK |
149 | _GLIBCXX_END_NAMESPACE_VERSION |
150 | } // namespace | |
01b0d4b8 | 151 | |
b2acb86f | 152 | #endif |