]>
Commit | Line | Data |
---|---|---|
01b0d4b8 PC |
1 | // Numeric extensions -*- C++ -*- |
2 | ||
3 | // Copyright (C) 2002 Free Software Foundation, Inc. | |
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 | ||
56 | #ifndef _EXT_NUMERIC | |
57 | #define _EXT_NUMERIC | |
58 | ||
59 | #pragma GCC system_header | |
60 | #include <bits/concept_check.h> | |
54c1bf78 | 61 | #include <numeric> |
01b0d4b8 PC |
62 | |
63 | #include <ext/functional> // For identity_element | |
64 | ||
65 | namespace __gnu_cxx | |
66 | { | |
67 | // Returns __x ** __n, where __n >= 0. _Note that "multiplication" | |
68 | // is required to be associative, but not necessarily commutative. | |
69 | ||
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); | |
76 | else { | |
77 | while ((__n & 1) == 0) { | |
78 | __n >>= 1; | |
79 | __x = __monoid_op(__x, __x); | |
80 | } | |
81 | ||
82 | _Tp __result = __x; | |
83 | __n >>= 1; | |
84 | while (__n != 0) { | |
85 | __x = __monoid_op(__x, __x); | |
86 | if ((__n & 1) != 0) | |
87 | __result = __monoid_op(__result, __x); | |
88 | __n >>= 1; | |
89 | } | |
90 | return __result; | |
91 | } | |
92 | } | |
93 | ||
94 | template<typename _Tp, typename _Integer> | |
95 | inline _Tp | |
96 | __power(_Tp __x, _Integer __n) | |
97 | { return __power(__x, __n, std::multiplies<_Tp>()); } | |
98 | ||
99 | // Alias for the internal name __power. Note that power is an extension, | |
100 | // not part of the C++ standard. | |
101 | ||
102 | template<typename _Tp, typename _Integer, typename _MonoidOperation> | |
103 | inline _Tp | |
104 | power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op) | |
105 | { return __power(__x, __n, __monoid_op); } | |
106 | ||
107 | template<typename _Tp, typename _Integer> | |
108 | inline _Tp | |
109 | power(_Tp __x, _Integer __n) | |
110 | { return __power(__x, __n); } | |
111 | ||
112 | // iota is not part of the C++ standard. It is an extension. | |
113 | ||
114 | template<typename _ForwardIter, typename _Tp> | |
115 | void | |
116 | iota(_ForwardIter __first, _ForwardIter __last, _Tp __value) | |
117 | { | |
118 | // concept requirements | |
119 | __glibcpp_function_requires(_Mutable_ForwardIteratorConcept<_ForwardIter>) | |
120 | __glibcpp_function_requires(_ConvertibleConcept<_Tp, | |
121 | typename std::iterator_traits<_ForwardIter>::value_type>) | |
122 | ||
123 | while (__first != __last) | |
124 | *__first++ = __value++; | |
125 | } | |
126 | ||
127 | } // namespace __gnu_cxx | |
128 | ||
129 | #endif /* _EXT_NUMERIC */ | |
130 |