#ifndef STORM_UTILITY_PARAMETRIC_H
#define STORM_UTILITY_PARAMETRIC_H

#include "storm/adapters/RationalFunctionAdapter.h"

#include <map>

namespace storm {
    namespace utility {
        namespace parametric {
            
            /*!
             * Access the type of variables from a given function type
             */
            template<typename FunctionType>
            struct VariableType { typedef void type; };
            
            /*!
             * Acess the type of coefficients from a given function type
             */
            template<typename FunctionType>
            struct CoefficientType { typedef void type; };
            
#ifdef STORM_HAVE_CARL
            template<>
            struct VariableType<storm::RationalFunction> { typedef storm::RationalFunctionVariable type; };
            template<>
            struct CoefficientType<storm::RationalFunction> { typedef storm::RationalFunctionCoefficient type; };
#endif

            template<typename FunctionType> using Valuation = std::map<typename VariableType<FunctionType>::type, typename CoefficientType<FunctionType>::type>;

            /*!
             * Evaluates the given function wrt. the given valuation
             */
            template<typename FunctionType>
            typename CoefficientType<FunctionType>::type evaluate(FunctionType const& function, Valuation<FunctionType> const& valuation);
            
            /*!
             *  Add all variables that occur in the given function to the the given set
             */
            template<typename FunctionType>
            void gatherOccurringVariables(FunctionType const& function, std::set<typename VariableType<FunctionType>::type>& variableSet);
            
            /*!
             *  Checks whether the function is linear (in one parameter)
             */
            template<typename FunctionType>
            bool isLinear(FunctionType const& function);
            
            /*!
             *  Checks whether the function is a multilinear polynomial, i.e., a polynomial which only considers variables with exponent at most 1
             */
            template<typename FunctionType>
            bool isMultiLinearPolynomial(FunctionType const& function);
            
        }
        
    }
}


#endif /* STORM_UTILITY_PARAMETRIC_H */