00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057 #ifndef _STL_NUMERIC_H
00058 #define _STL_NUMERIC_H 1
00059
00060 #include <bits/concept_check.h>
00061 #include <debug/debug.h>
00062 #include <bits/move.h>
00063
00064 #ifdef __GXX_EXPERIMENTAL_CXX0X__
00065
00066 namespace std _GLIBCXX_VISIBILITY(default)
00067 {
00068 _GLIBCXX_BEGIN_NAMESPACE_VERSION
00069
00070
00071
00072
00073
00074
00075
00076
00077
00078
00079
00080
00081 template<typename _ForwardIterator, typename _Tp>
00082 void
00083 iota(_ForwardIterator __first, _ForwardIterator __last, _Tp __value)
00084 {
00085
00086 __glibcxx_function_requires(_Mutable_ForwardIteratorConcept<
00087 _ForwardIterator>)
00088 __glibcxx_function_requires(_ConvertibleConcept<_Tp,
00089 typename iterator_traits<_ForwardIterator>::value_type>)
00090 __glibcxx_requires_valid_range(__first, __last);
00091
00092 for (; __first != __last; ++__first)
00093 {
00094 *__first = __value;
00095 ++__value;
00096 }
00097 }
00098
00099 _GLIBCXX_END_NAMESPACE_VERSION
00100 }
00101
00102 #endif
00103
00104 namespace std _GLIBCXX_VISIBILITY(default)
00105 {
00106 _GLIBCXX_BEGIN_NAMESPACE_ALGO
00107
00108
00109
00110
00111
00112
00113
00114
00115
00116
00117
00118
00119 template<typename _InputIterator, typename _Tp>
00120 inline _Tp
00121 accumulate(_InputIterator __first, _InputIterator __last, _Tp __init)
00122 {
00123
00124 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator>)
00125 __glibcxx_requires_valid_range(__first, __last);
00126
00127 for (; __first != __last; ++__first)
00128 __init = __init + *__first;
00129 return __init;
00130 }
00131
00132
00133
00134
00135
00136
00137
00138
00139
00140
00141
00142
00143
00144
00145 template<typename _InputIterator, typename _Tp, typename _BinaryOperation>
00146 inline _Tp
00147 accumulate(_InputIterator __first, _InputIterator __last, _Tp __init,
00148 _BinaryOperation __binary_op)
00149 {
00150
00151 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator>)
00152 __glibcxx_requires_valid_range(__first, __last);
00153
00154 for (; __first != __last; ++__first)
00155 __init = __binary_op(__init, *__first);
00156 return __init;
00157 }
00158
00159
00160
00161
00162
00163
00164
00165
00166
00167
00168
00169
00170
00171
00172
00173 template<typename _InputIterator1, typename _InputIterator2, typename _Tp>
00174 inline _Tp
00175 inner_product(_InputIterator1 __first1, _InputIterator1 __last1,
00176 _InputIterator2 __first2, _Tp __init)
00177 {
00178
00179 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator1>)
00180 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator2>)
00181 __glibcxx_requires_valid_range(__first1, __last1);
00182
00183 for (; __first1 != __last1; ++__first1, ++__first2)
00184 __init = __init + (*__first1 * *__first2);
00185 return __init;
00186 }
00187
00188
00189
00190
00191
00192
00193
00194
00195
00196
00197
00198
00199
00200
00201
00202
00203
00204 template<typename _InputIterator1, typename _InputIterator2, typename _Tp,
00205 typename _BinaryOperation1, typename _BinaryOperation2>
00206 inline _Tp
00207 inner_product(_InputIterator1 __first1, _InputIterator1 __last1,
00208 _InputIterator2 __first2, _Tp __init,
00209 _BinaryOperation1 __binary_op1,
00210 _BinaryOperation2 __binary_op2)
00211 {
00212
00213 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator1>)
00214 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator2>)
00215 __glibcxx_requires_valid_range(__first1, __last1);
00216
00217 for (; __first1 != __last1; ++__first1, ++__first2)
00218 __init = __binary_op1(__init, __binary_op2(*__first1, *__first2));
00219 return __init;
00220 }
00221
00222
00223
00224
00225
00226
00227
00228
00229
00230
00231
00232
00233
00234
00235
00236 template<typename _InputIterator, typename _OutputIterator>
00237 _OutputIterator
00238 partial_sum(_InputIterator __first, _InputIterator __last,
00239 _OutputIterator __result)
00240 {
00241 typedef typename iterator_traits<_InputIterator>::value_type _ValueType;
00242
00243
00244 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator>)
00245 __glibcxx_function_requires(_OutputIteratorConcept<_OutputIterator,
00246 _ValueType>)
00247 __glibcxx_requires_valid_range(__first, __last);
00248
00249 if (__first == __last)
00250 return __result;
00251 _ValueType __value = *__first;
00252 *__result = __value;
00253 while (++__first != __last)
00254 {
00255 __value = __value + *__first;
00256 *++__result = __value;
00257 }
00258 return ++__result;
00259 }
00260
00261
00262
00263
00264
00265
00266
00267
00268
00269
00270
00271
00272
00273
00274
00275 template<typename _InputIterator, typename _OutputIterator,
00276 typename _BinaryOperation>
00277 _OutputIterator
00278 partial_sum(_InputIterator __first, _InputIterator __last,
00279 _OutputIterator __result, _BinaryOperation __binary_op)
00280 {
00281 typedef typename iterator_traits<_InputIterator>::value_type _ValueType;
00282
00283
00284 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator>)
00285 __glibcxx_function_requires(_OutputIteratorConcept<_OutputIterator,
00286 _ValueType>)
00287 __glibcxx_requires_valid_range(__first, __last);
00288
00289 if (__first == __last)
00290 return __result;
00291 _ValueType __value = *__first;
00292 *__result = __value;
00293 while (++__first != __last)
00294 {
00295 __value = __binary_op(__value, *__first);
00296 *++__result = __value;
00297 }
00298 return ++__result;
00299 }
00300
00301
00302
00303
00304
00305
00306
00307
00308
00309
00310
00311
00312
00313
00314
00315 template<typename _InputIterator, typename _OutputIterator>
00316 _OutputIterator
00317 adjacent_difference(_InputIterator __first,
00318 _InputIterator __last, _OutputIterator __result)
00319 {
00320 typedef typename iterator_traits<_InputIterator>::value_type _ValueType;
00321
00322
00323 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator>)
00324 __glibcxx_function_requires(_OutputIteratorConcept<_OutputIterator,
00325 _ValueType>)
00326 __glibcxx_requires_valid_range(__first, __last);
00327
00328 if (__first == __last)
00329 return __result;
00330 _ValueType __value = *__first;
00331 *__result = __value;
00332 while (++__first != __last)
00333 {
00334 _ValueType __tmp = *__first;
00335 *++__result = __tmp - __value;
00336 __value = _GLIBCXX_MOVE(__tmp);
00337 }
00338 return ++__result;
00339 }
00340
00341
00342
00343
00344
00345
00346
00347
00348
00349
00350
00351
00352
00353
00354
00355
00356 template<typename _InputIterator, typename _OutputIterator,
00357 typename _BinaryOperation>
00358 _OutputIterator
00359 adjacent_difference(_InputIterator __first, _InputIterator __last,
00360 _OutputIterator __result, _BinaryOperation __binary_op)
00361 {
00362 typedef typename iterator_traits<_InputIterator>::value_type _ValueType;
00363
00364
00365 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator>)
00366 __glibcxx_function_requires(_OutputIteratorConcept<_OutputIterator,
00367 _ValueType>)
00368 __glibcxx_requires_valid_range(__first, __last);
00369
00370 if (__first == __last)
00371 return __result;
00372 _ValueType __value = *__first;
00373 *__result = __value;
00374 while (++__first != __last)
00375 {
00376 _ValueType __tmp = *__first;
00377 *++__result = __binary_op(__tmp, __value);
00378 __value = _GLIBCXX_MOVE(__tmp);
00379 }
00380 return ++__result;
00381 }
00382
00383 _GLIBCXX_END_NAMESPACE_ALGO
00384 }
00385
00386 #endif