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