/*

 *

 *

 * Copyright (c) 1994

 * Hewlett-Packard Company

 *

 * Copyright (c) 1996,1997

 * Silicon Graphics Computer Systems, Inc.

 *

 * Copyright (c) 1997

 * Moscow Center for SPARC Technology

 *

 * Copyright (c) 1999 

 * Boris Fomitchev

 *

 * This material is provided "as is", with absolutely no warranty expressed

 * or implied. Any use is at your own risk.

 *

 * Permission to use or copy this software for any purpose is hereby granted 

 * without fee, provided the above notices are retained on all copies.

 * Permission to modify the code and to distribute modified code is granted,

 * provided the above notices are retained, and a notice that the code was

 * modified is included with the above copyright notice.

 *

 */

#ifndef _STLP_NUMERIC_C

#define _STLP_NUMERIC_C

#ifndef _STLP_INTERNAL_NUMERIC_H

# include <stl/_numeric.h>

#endif

_STLP_BEGIN_NAMESPACE

template <class _InputIterator, class _OutputIterator, class _Tp,

          class _BinaryOperation>

_OutputIterator 

__partial_sum(_InputIterator __first, _InputIterator __last, 

              _OutputIterator __result, _Tp*, _BinaryOperation __binary_op)

{

  _STLP_DEBUG_CHECK(__check_range(__first, __last))

  if (__first == __last) return __result;

  *__result = *__first;

  _Tp __val = *__first;

  while (++__first != __last) {

    __val = __binary_op(__val, *__first);

    *++__result = __val;

  }

  return ++__result;

}

template <class _InputIterator, class _OutputIterator, class _Tp, 

          class _BinaryOperation>

_OutputIterator

__adjacent_difference(_InputIterator __first, _InputIterator __last, 

                      _OutputIterator __result, _Tp*,

                      _BinaryOperation __binary_op) {

  _STLP_DEBUG_CHECK(__check_range(__first, __last))

  if (__first == __last) return __result;

  *__result = *__first;

  _Tp __val = *__first;

  while (++__first != __last) {

    _Tp __tmp = *__first;

    *++__result = __binary_op(__tmp, __val);

    __val = __tmp;

  }

  return ++__result;

}

template <class _Tp, class _Integer, class _MonoidOperation>

_Tp __power(_Tp __x, _Integer __n, _MonoidOperation __opr) 

{

  _STLP_MPWFIX_TRY

  if (__n == 0)

    return __identity_element(__opr);

  else {

    while ((__n & 1) == 0) {

      __n >>= 1;

      __x = __opr(__x, __x);

    }

    _Tp __result = __x;

	_STLP_MPWFIX_TRY

    __n >>= 1;

    while (__n != 0) {

      __x = __opr(__x, __x);

      if ((__n & 1) != 0)

        __result = __opr(__result, __x);

      __n >>= 1;

    }

    return __result;

	_STLP_MPWFIX_CATCH

  }

  _STLP_MPWFIX_CATCH_ACTION(__x = _Tp())

}

_STLP_END_NAMESPACE

#endif /*  _STLP_NUMERIC_C */

// Local Variables:

// mode:C++

// End: