197403Sobrien// Numeric extensions -*- C++ -*- 297403Sobrien 3169691Skan// Copyright (C) 2002, 2004, 2005 Free Software Foundation, Inc. 497403Sobrien// 597403Sobrien// This file is part of the GNU ISO C++ Library. This library is free 697403Sobrien// software; you can redistribute it and/or modify it under the 797403Sobrien// terms of the GNU General Public License as published by the 897403Sobrien// Free Software Foundation; either version 2, or (at your option) 997403Sobrien// any later version. 1097403Sobrien 1197403Sobrien// This library is distributed in the hope that it will be useful, 1297403Sobrien// but WITHOUT ANY WARRANTY; without even the implied warranty of 1397403Sobrien// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 1497403Sobrien// GNU General Public License for more details. 1597403Sobrien 1697403Sobrien// You should have received a copy of the GNU General Public License along 1797403Sobrien// with this library; see the file COPYING. If not, write to the Free 18169691Skan// Software Foundation, 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, 1997403Sobrien// USA. 2097403Sobrien 2197403Sobrien// As a special exception, you may use this file as part of a free software 2297403Sobrien// library without restriction. Specifically, if other files instantiate 2397403Sobrien// templates or use macros or inline functions from this file, or you compile 2497403Sobrien// this file and link it with other files to produce an executable, this 2597403Sobrien// file does not by itself cause the resulting executable to be covered by 2697403Sobrien// the GNU General Public License. This exception does not however 2797403Sobrien// invalidate any other reasons why the executable file might be covered by 2897403Sobrien// the GNU General Public License. 2997403Sobrien 3097403Sobrien/* 3197403Sobrien * 3297403Sobrien * Copyright (c) 1994 3397403Sobrien * Hewlett-Packard Company 3497403Sobrien * 3597403Sobrien * Permission to use, copy, modify, distribute and sell this software 3697403Sobrien * and its documentation for any purpose is hereby granted without fee, 3797403Sobrien * provided that the above copyright notice appear in all copies and 3897403Sobrien * that both that copyright notice and this permission notice appear 3997403Sobrien * in supporting documentation. Hewlett-Packard Company makes no 4097403Sobrien * representations about the suitability of this software for any 4197403Sobrien * purpose. It is provided "as is" without express or implied warranty. 4297403Sobrien * 4397403Sobrien * 4497403Sobrien * Copyright (c) 1996 4597403Sobrien * Silicon Graphics Computer Systems, Inc. 4697403Sobrien * 4797403Sobrien * Permission to use, copy, modify, distribute and sell this software 4897403Sobrien * and its documentation for any purpose is hereby granted without fee, 4997403Sobrien * provided that the above copyright notice appear in all copies and 5097403Sobrien * that both that copyright notice and this permission notice appear 5197403Sobrien * in supporting documentation. Silicon Graphics makes no 5297403Sobrien * representations about the suitability of this software for any 5397403Sobrien * purpose. It is provided "as is" without express or implied warranty. 5497403Sobrien */ 5597403Sobrien 5697403Sobrien/** @file ext/numeric 5797403Sobrien * This file is a GNU extension to the Standard C++ Library (possibly 58169691Skan * containing extensions from the HP/SGI STL subset). 5997403Sobrien */ 6097403Sobrien 6197403Sobrien#ifndef _EXT_NUMERIC 62132720Skan#define _EXT_NUMERIC 1 6397403Sobrien 6497403Sobrien#pragma GCC system_header 65132720Skan 6697403Sobrien#include <bits/concept_check.h> 6797403Sobrien#include <numeric> 6897403Sobrien 6997403Sobrien#include <ext/functional> // For identity_element 7097403Sobrien 71169691Skan_GLIBCXX_BEGIN_NAMESPACE(__gnu_cxx) 72169691Skan 7397403Sobrien // Returns __x ** __n, where __n >= 0. _Note that "multiplication" 7497403Sobrien // is required to be associative, but not necessarily commutative. 7597403Sobrien template<typename _Tp, typename _Integer, typename _MonoidOperation> 7697403Sobrien _Tp 7797403Sobrien __power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op) 7897403Sobrien { 7997403Sobrien if (__n == 0) 8097403Sobrien return identity_element(__monoid_op); 81169691Skan else 82169691Skan { 83169691Skan while ((__n & 1) == 0) 84169691Skan { 85169691Skan __n >>= 1; 86169691Skan __x = __monoid_op(__x, __x); 87169691Skan } 8897403Sobrien 89169691Skan _Tp __result = __x; 9097403Sobrien __n >>= 1; 91169691Skan while (__n != 0) 92169691Skan { 93169691Skan __x = __monoid_op(__x, __x); 94169691Skan if ((__n & 1) != 0) 95169691Skan __result = __monoid_op(__result, __x); 96169691Skan __n >>= 1; 97169691Skan } 98169691Skan return __result; 9997403Sobrien } 10097403Sobrien } 10197403Sobrien 10297403Sobrien template<typename _Tp, typename _Integer> 10397403Sobrien inline _Tp 10497403Sobrien __power(_Tp __x, _Integer __n) 10597403Sobrien { return __power(__x, __n, std::multiplies<_Tp>()); } 10697403Sobrien 107102782Skan /** 108102782Skan * This is an SGI extension. 109102782Skan * @ingroup SGIextensions 110102782Skan * @doctodo 111102782Skan */ 112132720Skan // Alias for the internal name __power. Note that power is an extension, 113132720Skan // not part of the C++ standard. 11497403Sobrien template<typename _Tp, typename _Integer, typename _MonoidOperation> 11597403Sobrien inline _Tp 11697403Sobrien power(_Tp __x, _Integer __n, _MonoidOperation __monoid_op) 11797403Sobrien { return __power(__x, __n, __monoid_op); } 11897403Sobrien 119102782Skan /** 120102782Skan * This is an SGI extension. 121102782Skan * @ingroup SGIextensions 122102782Skan * @doctodo 123102782Skan */ 12497403Sobrien template<typename _Tp, typename _Integer> 12597403Sobrien inline _Tp 12697403Sobrien power(_Tp __x, _Integer __n) 12797403Sobrien { return __power(__x, __n); } 12897403Sobrien 129102782Skan /** 130102782Skan * This is an SGI extension. 131102782Skan * @ingroup SGIextensions 132102782Skan * @doctodo 133102782Skan */ 134132720Skan // iota is not part of the C++ standard. It is an extension. 13597403Sobrien template<typename _ForwardIter, typename _Tp> 136132720Skan void 13797403Sobrien iota(_ForwardIter __first, _ForwardIter __last, _Tp __value) 13897403Sobrien { 13997403Sobrien // concept requirements 140132720Skan __glibcxx_function_requires(_Mutable_ForwardIteratorConcept<_ForwardIter>) 141132720Skan __glibcxx_function_requires(_ConvertibleConcept<_Tp, 14297403Sobrien typename std::iterator_traits<_ForwardIter>::value_type>) 14397403Sobrien 14497403Sobrien while (__first != __last) 14597403Sobrien *__first++ = __value++; 14697403Sobrien } 14797403Sobrien 148169691Skan_GLIBCXX_END_NAMESPACE 149169691Skan 150132720Skan#endif 15197403Sobrien 152