include/c++/stl/stl/_numeric.c

   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
  37 template <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
  54 template <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
  73 template <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: