libalgebra-demo/libalgebra/ConstPower.h

/* *************************************************************

Copyright 2010 Terry Lyons, Stephen Buckley, Djalil Chafai, 
Greg Gyurkó and Arend Janssen. 

Distributed under the terms of the GNU General Public License, 
Version 3. (See accompanying file License.txt)

************************************************************* */




#pragma once
#ifndef ConstPower_h__
#define ConstPower_h__


template <unsigned arg, unsigned exp>
struct ConstPower
{
        static  enum
        {
                ans = ConstPower < arg, exp /2 >::ans * ConstPower < arg, exp /2 >::ans * ((exp%2 == 0) ? 1 : arg)
};
};

template < unsigned arg>
struct ConstPower < arg, 0>
{
static  enum
{
ans = 1
};
};


template < unsigned arg, unsigned exp >
struct TestConstPower
{
static  enum
{
intermediate = ConstPower < arg, exp>::ans * TestConstPower<arg, exp - 1 >::intermediate,
ans = (intermediate == ConstPower < arg, (ConstPower < exp, 2 >::ans + exp)/2>::ans) && TestConstPower<arg - 1,
exp >::ans
};
};

template < unsigned arg >
struct TestConstPower < arg, 0 >
{
static  enum
{
intermediate = ConstPower < arg, 0 >::ans,
ans = (intermediate == 1)
};
};

template < unsigned exp >
struct TestConstPower < 1, exp >
{
static  enum
{
intermediate = ConstPower < 1, exp >::ans,
ans = (intermediate == 1)
};
};

template <>
struct TestConstPower < 1, 0 >
{
static  enum
{
intermediate = ConstPower < 1, 0 >::ans,
ans = (intermediate == 1)
};
};


#endif // ConstPower_h__