1 // Numeric extensions -*- C++ -*-
3 // Copyright (C) 2002-2022 Free Software Foundation, Inc.
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 3, or (at your option)
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.
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.
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/>.
28 * Hewlett-Packard Company
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.
40 * Silicon Graphics Computer Systems, Inc.
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.
52 * This file is a GNU extension to the Standard C++ Library (possibly
53 * containing extensions from the HP/SGI STL subset).
57 #define _EXT_NUMERIC 1
59 #pragma GCC system_header
61 #include <bits/concept_check.h>
64 #include <ext/functional> // For identity_element
66 namespace __gnu_cxx _GLIBCXX_VISIBILITY(default)
68 _GLIBCXX_BEGIN_NAMESPACE_VERSION
70 // Returns __x ** __n, where __n >= 0. _Note that "multiplication"
71 // is required to be associative, but not necessarily commutative.
72 template<typename _Tp, typename _Integer, typename _MonoidOperation>
74 __power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op)
77 return identity_element(__monoid_op);
80 while ((__n & 1) == 0)
83 __x = __monoid_op(__x, __x);
90 __x = __monoid_op(__x, __x);
92 __result = __monoid_op(__result, __x);
99 template<typename _Tp, typename _Integer>
101 __power(_Tp __x, _Integer __n)
102 { return __power(__x, __n, std::multiplies<_Tp>()); }
105 * This is an SGI extension.
106 * @ingroup SGIextensions
109 // Alias for the internal name __power. Note that power is an extension,
110 // not part of the C++ standard.
111 template<typename _Tp, typename _Integer, typename _MonoidOperation>
113 power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op)
114 { return __power(__x, __n, __monoid_op); }
117 * This is an SGI extension.
118 * @ingroup SGIextensions
121 template<typename _Tp, typename _Integer>
123 power(_Tp __x, _Integer __n)
124 { return __power(__x, __n); }
126 #if __cplusplus >= 201103L
130 * This is an SGI extension.
131 * @ingroup SGIextensions
134 // iota is not part of the standard until C++11. It is an extension.
135 template<typename _ForwardIter, typename _Tp>
137 iota(_ForwardIter __first, _ForwardIter __last, _Tp __value)
139 // concept requirements
140 __glibcxx_function_requires(_Mutable_ForwardIteratorConcept<_ForwardIter>)
141 __glibcxx_function_requires(_ConvertibleConcept<_Tp,
142 typename std::iterator_traits<_ForwardIter>::value_type>)
144 while (__first != __last)
145 *__first++ = __value++;
149 _GLIBCXX_END_NAMESPACE_VERSION