1/*
2 *
3 *
4 * Copyright (c) 1994
5 * Hewlett-Packard Company
6 *
7 * Copyright (c) 1996,1997
8 * Silicon Graphics Computer Systems, Inc.
9 *
10 * Copyright (c) 1997
11 * Moscow Center for SPARC Technology
12 *
13 * Copyright (c) 1999
14 * Boris Fomitchev
15 *
16 * This material is provided "as is", with absolutely no warranty expressed
17 * or implied. Any use is at your own risk.
18 *
19 * Permission to use or copy this software for any purpose is hereby granted
20 * without fee, provided the above notices are retained on all copies.
21 * Permission to modify the code and to distribute modified code is granted,
22 * provided the above notices are retained, and a notice that the code was
23 * modified is included with the above copyright notice.
24 *
25 */
26#ifndef _STLP_NUMERIC_C
27#define _STLP_NUMERIC_C
28
29#ifndef _STLP_INTERNAL_NUMERIC_H
30# include <stl/_numeric.h>
31#endif
32
33_STLP_BEGIN_NAMESPACE
34
35_STLP_MOVE_TO_PRIV_NAMESPACE
36
37template <class _InputIterator, class _OutputIterator, class _Tp,
38          class _BinaryOperation>
39_OutputIterator
40__partial_sum(_InputIterator __first, _InputIterator __last,
41              _OutputIterator __result, _Tp*, _BinaryOperation __binary_op) {
42  _STLP_DEBUG_CHECK(__check_range(__first, __last))
43  if (__first == __last) return __result;
44  *__result = *__first;
45
46  _Tp __val = *__first;
47  while (++__first != __last) {
48    __val = __binary_op(__val, *__first);
49    *++__result = __val;
50  }
51  return ++__result;
52}
53
54template <class _InputIterator, class _OutputIterator, class _Tp,
55          class _BinaryOperation>
56_OutputIterator
57__adjacent_difference(_InputIterator __first, _InputIterator __last,
58                      _OutputIterator __result, _Tp*,
59                      _BinaryOperation __binary_op) {
60  _STLP_DEBUG_CHECK(__check_range(__first, __last))
61  if (__first == __last) return __result;
62  *__result = *__first;
63  _Tp __val = *__first;
64  while (++__first != __last) {
65    _Tp __tmp = *__first;
66    *++__result = __binary_op(__tmp, __val);
67    __val = __tmp;
68  }
69  return ++__result;
70}
71
72
73template <class _Tp, class _Integer, class _MonoidOperation>
74_Tp __power(_Tp __x, _Integer __n, _MonoidOperation __opr) {
75  _STLP_MPWFIX_TRY
76  if (__n == 0)
77    return __identity_element(__opr);
78  else {
79    while ((__n & 1) == 0) {
80      __n >>= 1;
81      __x = __opr(__x, __x);
82    }
83    _Tp __result = __x;
84  _STLP_MPWFIX_TRY
85    __n >>= 1;
86    while (__n != 0) {
87      __x = __opr(__x, __x);
88      if ((__n & 1) != 0)
89        __result = __opr(__result, __x);
90      __n >>= 1;
91    }
92    return __result;
93  _STLP_MPWFIX_CATCH
94  }
95  _STLP_MPWFIX_CATCH_ACTION(__x = _Tp())
96}
97
98_STLP_MOVE_TO_STD_NAMESPACE
99
100_STLP_END_NAMESPACE
101
102#endif /*  _STLP_NUMERIC_C */
103
104// Local Variables:
105// mode:C++
106// End:
107