]> git.ipfire.org Git - thirdparty/gcc.git/blame - libstdc++-v3/include/ext/numeric
doxygroups.cc: Add std::tr1 namespace.
[thirdparty/gcc.git] / libstdc++-v3 / include / ext / numeric
CommitLineData
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
71namespace __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