]>
Commit | Line | Data |
---|---|---|
01b0d4b8 PC |
1 | // Numeric extensions -*- C++ -*- |
2 | ||
d962e073 | 3 | // Copyright (C) 2002, 2004 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 | |
8 | // Free Software Foundation; either version 2, 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 | // You should have received a copy of the GNU General Public License along | |
17 | // with this library; see the file COPYING. If not, write to the Free | |
18 | // Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307, | |
19 | // USA. | |
20 | ||
21 | // As a special exception, you may use this file as part of a free software | |
22 | // library without restriction. Specifically, if other files instantiate | |
23 | // templates or use macros or inline functions from this file, or you compile | |
24 | // this file and link it with other files to produce an executable, this | |
25 | // file does not by itself cause the resulting executable to be covered by | |
26 | // the GNU General Public License. This exception does not however | |
27 | // invalidate any other reasons why the executable file might be covered by | |
28 | // the GNU General Public License. | |
29 | ||
30 | /* | |
31 | * | |
32 | * Copyright (c) 1994 | |
33 | * Hewlett-Packard Company | |
34 | * | |
35 | * Permission to use, copy, modify, distribute and sell this software | |
36 | * and its documentation for any purpose is hereby granted without fee, | |
37 | * provided that the above copyright notice appear in all copies and | |
38 | * that both that copyright notice and this permission notice appear | |
39 | * in supporting documentation. Hewlett-Packard Company makes no | |
40 | * representations about the suitability of this software for any | |
41 | * purpose. It is provided "as is" without express or implied warranty. | |
42 | * | |
43 | * | |
44 | * Copyright (c) 1996 | |
45 | * Silicon Graphics Computer Systems, Inc. | |
46 | * | |
47 | * Permission to use, copy, modify, distribute and sell this software | |
48 | * and its documentation for any purpose is hereby granted without fee, | |
49 | * provided that the above copyright notice appear in all copies and | |
50 | * that both that copyright notice and this permission notice appear | |
51 | * in supporting documentation. Silicon Graphics makes no | |
52 | * representations about the suitability of this software for any | |
53 | * purpose. It is provided "as is" without express or implied warranty. | |
54 | */ | |
55 | ||
ffe94f83 PE |
56 | /** @file ext/numeric |
57 | * This file is a GNU extension to the Standard C++ Library (possibly | |
0aa06b18 | 58 | * containing extensions from the HP/SGI STL subset). |
ffe94f83 PE |
59 | */ |
60 | ||
01b0d4b8 | 61 | #ifndef _EXT_NUMERIC |
b2acb86f | 62 | #define _EXT_NUMERIC 1 |
01b0d4b8 PC |
63 | |
64 | #pragma GCC system_header | |
b2acb86f | 65 | |
01b0d4b8 | 66 | #include <bits/concept_check.h> |
54c1bf78 | 67 | #include <numeric> |
01b0d4b8 PC |
68 | |
69 | #include <ext/functional> // For identity_element | |
70 | ||
71 | namespace __gnu_cxx | |
72 | { | |
73 | // Returns __x ** __n, where __n >= 0. _Note that "multiplication" | |
fa30fe72 | 74 | // is required to be associative, but not necessarily commutative. |
01b0d4b8 PC |
75 | template<typename _Tp, typename _Integer, typename _MonoidOperation> |
76 | _Tp | |
77 | __power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op) | |
78 | { | |
79 | if (__n == 0) | |
80 | return identity_element(__monoid_op); | |
d962e073 PC |
81 | else |
82 | { | |
83 | while ((__n & 1) == 0) | |
84 | { | |
85 | __n >>= 1; | |
86 | __x = __monoid_op(__x, __x); | |
87 | } | |
01b0d4b8 | 88 | |
d962e073 | 89 | _Tp __result = __x; |
01b0d4b8 | 90 | __n >>= 1; |
d962e073 PC |
91 | while (__n != 0) |
92 | { | |
93 | __x = __monoid_op(__x, __x); | |
94 | if ((__n & 1) != 0) | |
95 | __result = __monoid_op(__result, __x); | |
96 | __n >>= 1; | |
97 | } | |
98 | return __result; | |
01b0d4b8 | 99 | } |
01b0d4b8 PC |
100 | } |
101 | ||
102 | template<typename _Tp, typename _Integer> | |
103 | inline _Tp | |
104 | __power(_Tp __x, _Integer __n) | |
105 | { return __power(__x, __n, std::multiplies<_Tp>()); } | |
106 | ||
fc7f0a80 PE |
107 | /** |
108 | * This is an SGI extension. | |
109 | * @ingroup SGIextensions | |
110 | * @doctodo | |
111 | */ | |
b2acb86f BK |
112 | // Alias for the internal name __power. Note that power is an extension, |
113 | // not part of the C++ standard. | |
01b0d4b8 PC |
114 | template<typename _Tp, typename _Integer, typename _MonoidOperation> |
115 | inline _Tp | |
116 | power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op) | |
117 | { return __power(__x, __n, __monoid_op); } | |
118 | ||
fc7f0a80 PE |
119 | /** |
120 | * This is an SGI extension. | |
121 | * @ingroup SGIextensions | |
122 | * @doctodo | |
123 | */ | |
01b0d4b8 PC |
124 | template<typename _Tp, typename _Integer> |
125 | inline _Tp | |
126 | power(_Tp __x, _Integer __n) | |
127 | { return __power(__x, __n); } | |
128 | ||
fc7f0a80 PE |
129 | /** |
130 | * This is an SGI extension. | |
131 | * @ingroup SGIextensions | |
132 | * @doctodo | |
133 | */ | |
b2acb86f | 134 | // iota is not part of the C++ standard. 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 | } | |
01b0d4b8 PC |
147 | } // namespace __gnu_cxx |
148 | ||
b2acb86f | 149 | #endif |
01b0d4b8 | 150 |