1// Numeric extensions -*- C++ -*-
2
3// Copyright (C) 2002 Free Software Foundation, Inc.
4//
5// This file is part of the GNU ISO C++ Library.  This library is free
6// software; you can redistribute it and/or modify it under the
7// terms of the GNU General Public License as published by the
8// Free Software Foundation; either version 2, or (at your option)
9// any later version.
10
11// This library is distributed in the hope that it will be useful,
12// but WITHOUT ANY WARRANTY; without even the implied warranty of
13// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14// GNU General Public License for more details.
15
16// You should have received a copy of the GNU General Public License along
17// with this library; see the file COPYING.  If not, write to the Free
18// Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307,
19// USA.
20
21// As a special exception, you may use this file as part of a free software
22// library without restriction.  Specifically, if other files instantiate
23// templates or use macros or inline functions from this file, or you compile
24// this file and link it with other files to produce an executable, this
25// file does not by itself cause the resulting executable to be covered by
26// the GNU General Public License.  This exception does not however
27// invalidate any other reasons why the executable file might be covered by
28// the GNU General Public License.
29
30/*
31 *
32 * Copyright (c) 1994
33 * Hewlett-Packard Company
34 *
35 * Permission to use, copy, modify, distribute and sell this software
36 * and its documentation for any purpose is hereby granted without fee,
37 * provided that the above copyright notice appear in all copies and
38 * that both that copyright notice and this permission notice appear
39 * in supporting documentation.  Hewlett-Packard Company makes no
40 * representations about the suitability of this software for any
41 * purpose.  It is provided "as is" without express or implied warranty.
42 *
43 *
44 * Copyright (c) 1996
45 * Silicon Graphics Computer Systems, Inc.
46 *
47 * Permission to use, copy, modify, distribute and sell this software
48 * and its documentation for any purpose is hereby granted without fee,
49 * provided that the above copyright notice appear in all copies and
50 * that both that copyright notice and this permission notice appear
51 * in supporting documentation.  Silicon Graphics makes no
52 * representations about the suitability of this software for any
53 * purpose.  It is provided "as is" without express or implied warranty.
54 */
55
56/** @file ext/numeric
57 *  This file is a GNU extension to the Standard C++ Library (possibly
58 *  containing extensions from the HP/SGI STL subset).  You should only
59 *  include this header if you are using GCC 3 or later.
60 */
61
62#ifndef _EXT_NUMERIC
63#define _EXT_NUMERIC 1
64
65#pragma GCC system_header
66
67#include <bits/concept_check.h>
68#include <numeric>
69
70#include <ext/functional> // For identity_element
71
72namespace __gnu_cxx
73{
74  // Returns __x ** __n, where __n >= 0.  _Note that "multiplication"
75  // is required to be associative, but not necessarily commutative.
76  template<typename _Tp, typename _Integer, typename _MonoidOperation>
77    _Tp
78    __power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op)
79    {
80      if (__n == 0)
81	return identity_element(__monoid_op);
82      else {
83	while ((__n & 1) == 0) {
84	  __n >>= 1;
85	  __x = __monoid_op(__x, __x);
86	}
87
88	_Tp __result = __x;
89	__n >>= 1;
90	while (__n != 0) {
91	  __x = __monoid_op(__x, __x);
92	  if ((__n & 1) != 0)
93	    __result = __monoid_op(__result, __x);
94	  __n >>= 1;
95	}
96	return __result;
97      }
98    }
99
100  template<typename _Tp, typename _Integer>
101    inline _Tp
102    __power(_Tp __x, _Integer __n)
103    { return __power(__x, __n, std::multiplies<_Tp>()); }
104
105  /**
106   *  This is an SGI extension.
107   *  @ingroup SGIextensions
108   *  @doctodo
109  */
110  // Alias for the internal name __power.  Note that power is an extension,
111  // not part of the C++ standard.
112  template<typename _Tp, typename _Integer, typename _MonoidOperation>
113    inline _Tp
114    power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op)
115    { return __power(__x, __n, __monoid_op); }
116
117  /**
118   *  This is an SGI extension.
119   *  @ingroup SGIextensions
120   *  @doctodo
121  */
122  template<typename _Tp, typename _Integer>
123    inline _Tp
124    power(_Tp __x, _Integer __n)
125    { return __power(__x, __n); }
126
127  /**
128   *  This is an SGI extension.
129   *  @ingroup SGIextensions
130   *  @doctodo
131  */
132  // iota is not part of the C++ standard.  It is an extension.
133  template<typename _ForwardIter, typename _Tp>
134    void
135    iota(_ForwardIter __first, _ForwardIter __last, _Tp __value)
136    {
137      // concept requirements
138      __glibcxx_function_requires(_Mutable_ForwardIteratorConcept<_ForwardIter>)
139      __glibcxx_function_requires(_ConvertibleConcept<_Tp,
140	    typename std::iterator_traits<_ForwardIter>::value_type>)
141
142      while (__first != __last)
143	*__first++ = __value++;
144    }
145} // namespace __gnu_cxx
146
147#endif
148
149