]>
Commit | Line | Data |
---|---|---|
01b0d4b8 PC |
1 | // Numeric extensions -*- C++ -*- |
2 | ||
748086b7 | 3 | // Copyright (C) 2002, 2004, 2005, 2009 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 | ||
3cbc7af0 BK |
66 | _GLIBCXX_BEGIN_NAMESPACE(__gnu_cxx) |
67 | ||
01b0d4b8 | 68 | // Returns __x ** __n, where __n >= 0. _Note that "multiplication" |
fa30fe72 | 69 | // is required to be associative, but not necessarily commutative. |
01b0d4b8 PC |
70 | template<typename _Tp, typename _Integer, typename _MonoidOperation> |
71 | _Tp | |
72 | __power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op) | |
73 | { | |
74 | if (__n == 0) | |
75 | return identity_element(__monoid_op); | |
d962e073 PC |
76 | else |
77 | { | |
78 | while ((__n & 1) == 0) | |
79 | { | |
80 | __n >>= 1; | |
81 | __x = __monoid_op(__x, __x); | |
82 | } | |
01b0d4b8 | 83 | |
d962e073 | 84 | _Tp __result = __x; |
01b0d4b8 | 85 | __n >>= 1; |
d962e073 PC |
86 | while (__n != 0) |
87 | { | |
88 | __x = __monoid_op(__x, __x); | |
89 | if ((__n & 1) != 0) | |
90 | __result = __monoid_op(__result, __x); | |
91 | __n >>= 1; | |
92 | } | |
93 | return __result; | |
01b0d4b8 | 94 | } |
01b0d4b8 PC |
95 | } |
96 | ||
97 | template<typename _Tp, typename _Integer> | |
98 | inline _Tp | |
99 | __power(_Tp __x, _Integer __n) | |
100 | { return __power(__x, __n, std::multiplies<_Tp>()); } | |
101 | ||
fc7f0a80 PE |
102 | /** |
103 | * This is an SGI extension. | |
104 | * @ingroup SGIextensions | |
105 | * @doctodo | |
106 | */ | |
b2acb86f BK |
107 | // Alias for the internal name __power. Note that power is an extension, |
108 | // not part of the C++ standard. | |
01b0d4b8 PC |
109 | template<typename _Tp, typename _Integer, typename _MonoidOperation> |
110 | inline _Tp | |
111 | power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op) | |
112 | { return __power(__x, __n, __monoid_op); } | |
113 | ||
fc7f0a80 PE |
114 | /** |
115 | * This is an SGI extension. | |
116 | * @ingroup SGIextensions | |
117 | * @doctodo | |
118 | */ | |
01b0d4b8 PC |
119 | template<typename _Tp, typename _Integer> |
120 | inline _Tp | |
121 | power(_Tp __x, _Integer __n) | |
122 | { return __power(__x, __n); } | |
123 | ||
fc7f0a80 PE |
124 | /** |
125 | * This is an SGI extension. | |
126 | * @ingroup SGIextensions | |
127 | * @doctodo | |
128 | */ | |
b2acb86f | 129 | // iota is not part of the C++ standard. It is an extension. |
01b0d4b8 | 130 | template<typename _ForwardIter, typename _Tp> |
fa30fe72 | 131 | void |
01b0d4b8 PC |
132 | iota(_ForwardIter __first, _ForwardIter __last, _Tp __value) |
133 | { | |
134 | // concept requirements | |
3d7c150e BK |
135 | __glibcxx_function_requires(_Mutable_ForwardIteratorConcept<_ForwardIter>) |
136 | __glibcxx_function_requires(_ConvertibleConcept<_Tp, | |
01b0d4b8 PC |
137 | typename std::iterator_traits<_ForwardIter>::value_type>) |
138 | ||
139 | while (__first != __last) | |
140 | *__first++ = __value++; | |
141 | } | |
3cbc7af0 BK |
142 | |
143 | _GLIBCXX_END_NAMESPACE | |
01b0d4b8 | 144 | |
b2acb86f | 145 | #endif |
01b0d4b8 | 146 |