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
00063 #ifdef __GXX_EXPERIMENTAL_CXX0X__
00064
00065 _GLIBCXX_BEGIN_NAMESPACE(std)
00066
00067
00068
00069
00070
00071
00072
00073
00074
00075
00076
00077
00078 template<typename _ForwardIterator, typename _Tp>
00079 void
00080 iota(_ForwardIterator __first, _ForwardIterator __last, _Tp __value)
00081 {
00082
00083 __glibcxx_function_requires(_Mutable_ForwardIteratorConcept<
00084 _ForwardIterator>)
00085 __glibcxx_function_requires(_ConvertibleConcept<_Tp,
00086 typename iterator_traits<_ForwardIterator>::value_type>)
00087 __glibcxx_requires_valid_range(__first, __last);
00088
00089 for (; __first != __last; ++__first)
00090 {
00091 *__first = __value;
00092 ++__value;
00093 }
00094 }
00095
00096 _GLIBCXX_END_NAMESPACE
00097
00098 #endif
00099
00100 _GLIBCXX_BEGIN_NESTED_NAMESPACE(std, _GLIBCXX_STD_P)
00101
00102
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113 template<typename _InputIterator, typename _Tp>
00114 inline _Tp
00115 accumulate(_InputIterator __first, _InputIterator __last, _Tp __init)
00116 {
00117
00118 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator>)
00119 __glibcxx_requires_valid_range(__first, __last);
00120
00121 for (; __first != __last; ++__first)
00122 __init = __init + *__first;
00123 return __init;
00124 }
00125
00126
00127
00128
00129
00130
00131
00132
00133
00134
00135
00136
00137
00138
00139 template<typename _InputIterator, typename _Tp, typename _BinaryOperation>
00140 inline _Tp
00141 accumulate(_InputIterator __first, _InputIterator __last, _Tp __init,
00142 _BinaryOperation __binary_op)
00143 {
00144
00145 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator>)
00146 __glibcxx_requires_valid_range(__first, __last);
00147
00148 for (; __first != __last; ++__first)
00149 __init = __binary_op(__init, *__first);
00150 return __init;
00151 }
00152
00153
00154
00155
00156
00157
00158
00159
00160
00161
00162
00163
00164
00165
00166
00167 template<typename _InputIterator1, typename _InputIterator2, typename _Tp>
00168 inline _Tp
00169 inner_product(_InputIterator1 __first1, _InputIterator1 __last1,
00170 _InputIterator2 __first2, _Tp __init)
00171 {
00172
00173 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator1>)
00174 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator2>)
00175 __glibcxx_requires_valid_range(__first1, __last1);
00176
00177 for (; __first1 != __last1; ++__first1, ++__first2)
00178 __init = __init + (*__first1 * *__first2);
00179 return __init;
00180 }
00181
00182
00183
00184
00185
00186
00187
00188
00189
00190
00191
00192
00193
00194
00195
00196
00197
00198 template<typename _InputIterator1, typename _InputIterator2, typename _Tp,
00199 typename _BinaryOperation1, typename _BinaryOperation2>
00200 inline _Tp
00201 inner_product(_InputIterator1 __first1, _InputIterator1 __last1,
00202 _InputIterator2 __first2, _Tp __init,
00203 _BinaryOperation1 __binary_op1,
00204 _BinaryOperation2 __binary_op2)
00205 {
00206
00207 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator1>)
00208 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator2>)
00209 __glibcxx_requires_valid_range(__first1, __last1);
00210
00211 for (; __first1 != __last1; ++__first1, ++__first2)
00212 __init = __binary_op1(__init, __binary_op2(*__first1, *__first2));
00213 return __init;
00214 }
00215
00216
00217
00218
00219
00220
00221
00222
00223
00224
00225
00226
00227
00228
00229
00230 template<typename _InputIterator, typename _OutputIterator>
00231 _OutputIterator
00232 partial_sum(_InputIterator __first, _InputIterator __last,
00233 _OutputIterator __result)
00234 {
00235 typedef typename iterator_traits<_InputIterator>::value_type _ValueType;
00236
00237
00238 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator>)
00239 __glibcxx_function_requires(_OutputIteratorConcept<_OutputIterator,
00240 _ValueType>)
00241 __glibcxx_requires_valid_range(__first, __last);
00242
00243 if (__first == __last)
00244 return __result;
00245 _ValueType __value = *__first;
00246 *__result = __value;
00247 while (++__first != __last)
00248 {
00249 __value = __value + *__first;
00250 *++__result = __value;
00251 }
00252 return ++__result;
00253 }
00254
00255
00256
00257
00258
00259
00260
00261
00262
00263
00264
00265
00266
00267
00268
00269 template<typename _InputIterator, typename _OutputIterator,
00270 typename _BinaryOperation>
00271 _OutputIterator
00272 partial_sum(_InputIterator __first, _InputIterator __last,
00273 _OutputIterator __result, _BinaryOperation __binary_op)
00274 {
00275 typedef typename iterator_traits<_InputIterator>::value_type _ValueType;
00276
00277
00278 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator>)
00279 __glibcxx_function_requires(_OutputIteratorConcept<_OutputIterator,
00280 _ValueType>)
00281 __glibcxx_requires_valid_range(__first, __last);
00282
00283 if (__first == __last)
00284 return __result;
00285 _ValueType __value = *__first;
00286 *__result = __value;
00287 while (++__first != __last)
00288 {
00289 __value = __binary_op(__value, *__first);
00290 *++__result = __value;
00291 }
00292 return ++__result;
00293 }
00294
00295
00296
00297
00298
00299
00300
00301
00302
00303
00304
00305
00306 template<typename _InputIterator, typename _OutputIterator>
00307 _OutputIterator
00308 adjacent_difference(_InputIterator __first,
00309 _InputIterator __last, _OutputIterator __result)
00310 {
00311 typedef typename iterator_traits<_InputIterator>::value_type _ValueType;
00312
00313
00314 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator>)
00315 __glibcxx_function_requires(_OutputIteratorConcept<_OutputIterator,
00316 _ValueType>)
00317 __glibcxx_requires_valid_range(__first, __last);
00318
00319 if (__first == __last)
00320 return __result;
00321 _ValueType __value = *__first;
00322 *__result = __value;
00323 while (++__first != __last)
00324 {
00325 _ValueType __tmp = *__first;
00326 *++__result = __tmp - __value;
00327 __value = __tmp;
00328 }
00329 return ++__result;
00330 }
00331
00332
00333
00334
00335
00336
00337
00338
00339
00340
00341
00342
00343
00344 template<typename _InputIterator, typename _OutputIterator,
00345 typename _BinaryOperation>
00346 _OutputIterator
00347 adjacent_difference(_InputIterator __first, _InputIterator __last,
00348 _OutputIterator __result, _BinaryOperation __binary_op)
00349 {
00350 typedef typename iterator_traits<_InputIterator>::value_type _ValueType;
00351
00352
00353 __glibcxx_function_requires(_InputIteratorConcept<_InputIterator>)
00354 __glibcxx_function_requires(_OutputIteratorConcept<_OutputIterator,
00355 _ValueType>)
00356 __glibcxx_requires_valid_range(__first, __last);
00357
00358 if (__first == __last)
00359 return __result;
00360 _ValueType __value = *__first;
00361 *__result = __value;
00362 while (++__first != __last)
00363 {
00364 _ValueType __tmp = *__first;
00365 *++__result = __binary_op(__tmp, __value);
00366 __value = __tmp;
00367 }
00368 return ++__result;
00369 }
00370
00371 _GLIBCXX_END_NESTED_NAMESPACE
00372
00373 #endif