gsTemplateTools.h

Go to the documentation of this file.

#pragma once

#include <gsCore/gsExport.h>
#include <utility>
#include <complex>

namespace gismo
{

namespace util {

#if __cplusplus >= 201103L
//see also http://lists.boost.org/Archives/boost/2009/04/151209.php
// has_move_constructor is not working with MSVC up to VS2019
template <typename T> struct has_move_constructor
{
    typedef char yes[1];
    typedef char no[2];
    
    struct AmbiguousConverter
    {
        operator T&& ();
        operator const T& ();
    };
    template <typename C> static no& test(decltype( new C( AmbiguousConverter{} )));
    template <typename> static yes& test(...);
    enum { value = (sizeof(test<T>(0)) == sizeof(yes)) };
};
#endif

#if __cplusplus >= 201103L || _MSC_VER >= 1600

using std::conditional;
using std::enable_if;
using std::false_type;
using std::integral_constant;
using std::is_base_of;
using std::is_integral;
using std::is_same;
using std::reference_wrapper;
using std::remove_const;
using std::remove_cv;
using std::remove_volatile;
using std::true_type;
using std::make_unsigned;
using std::make_signed;
using std::is_signed;

# define GS_BIND1ST(_op,_arg) std::bind(_op, _arg, std::placeholders::_1)
# define GS_BIND2ND(_op,_arg) std::bind(_op, std::placeholders::_1, _arg)

#else

# define GS_BIND1ST(_op,_arg) std::bind1st(_op,_arg)
# define GS_BIND2ND(_op,_arg) std::bind2nd(_op,_arg)

// template <typename T> struct has_move_constructor { enum { value = 0 }; };

template<bool B, class T, class F> struct conditional { typedef T type; };
template<class T, class F> struct conditional<false, T, F> { typedef F type; };

template<bool B, class T = void> struct enable_if {};
template<class T> struct enable_if<true, T> { typedef T type;};

template<class T, class U> struct is_same { enum { value = 0 }; };
template<class T>          struct is_same<T, T> { enum { value = 1 }; };

template<typename U> struct is_pointer     { static const bool value = false; };
template<typename U> struct is_pointer<U*> { static const bool value = true ; };

template <typename B, typename D> struct Host
{ operator B*() const; operator D*(); };
template <typename B, typename D>
struct is_base_of
{
    typedef char (&yes)[1];
    typedef char (&no)[2];
    template <typename T> static yes check(D*, T);
    static no check(B*, int);
    static const bool value = sizeof(check(Host<B,D>(), int())) == sizeof(yes);
};

template <class T>
class reference_wrapper
{
public:
    typedef T type;
    
    reference_wrapper(const T& ref) : _ptr(&const_cast<T&>(ref)) { }
    reference_wrapper(const reference_wrapper&o) : _ptr(o._ptr)  { }
    
    reference_wrapper& operator=(const reference_wrapper& o)
    { _ptr = o._ptr; return *this; }
    
    operator T& () const { return *_ptr; }
    T& get() const { return *_ptr; }

private:
    T* _ptr;
};

template<class T, T v>
    struct integral_constant {
    static const T value = v;
    typedef T value_type;
    typedef integral_constant type;
    value_type operator()() const { return value; }
};

typedef integral_constant<bool, true>  true_type;
typedef integral_constant<bool, false> false_type;

template<typename> struct is_integral_base             : false_type {};

template<> struct is_integral_base<bool>               : true_type {};
template<> struct is_integral_base<char>               : true_type {};
template<> struct is_integral_base<signed char>        : true_type {};
template<> struct is_integral_base<unsigned char>      : true_type {};
template<> struct is_integral_base<wchar_t>            : true_type {};
template<> struct is_integral_base<short>              : true_type {};
template<> struct is_integral_base<int>                : true_type {};
template<> struct is_integral_base<long>               : true_type {};
template<> struct is_integral_base<long long>          : true_type {};
template<> struct is_integral_base<unsigned short>     : true_type {};
template<> struct is_integral_base<unsigned int>       : true_type {};
template<> struct is_integral_base<unsigned long>      : true_type {};
template<> struct is_integral_base<unsigned long long> : true_type {};

template< class T > struct remove_const                  { typedef T type; };
template< class T > struct remove_const<const T>         { typedef T type; };

template< class T > struct remove_volatile               { typedef T type; };
template< class T > struct remove_volatile<volatile T>   { typedef T type; };

template< class T >
struct remove_cv { typedef typename remove_volatile<typename remove_const<T>::type>::type type; };

template<typename T> struct is_integral: is_integral_base<typename remove_cv<T>::type> {};

template<class T>
struct make_unsigned;
#define GISMO_MAKE_UNSIGNED(signed_type)     \
template<>                                   \
struct make_unsigned<signed signed_type> {   \
    typedef unsigned signed_type type;       \
};                                           \
template<>                                   \
struct make_unsigned<unsigned signed_type> { \
    typedef unsigned signed_type type;       \
};
template<>
struct make_unsigned<char> {
    typedef unsigned char type;
};
GISMO_MAKE_UNSIGNED(char)
GISMO_MAKE_UNSIGNED(short)
GISMO_MAKE_UNSIGNED(int)
GISMO_MAKE_UNSIGNED(long)
GISMO_MAKE_UNSIGNED(long long)
#undef GISMO_MAKE_UNSIGNED

template<class T>
struct make_signed;
#define GISMO_MAKE_SIGNED(unsigned_type)     \
template<>                                   \
struct make_signed<signed unsigned_type> {   \
    typedef signed unsigned_type type;       \
};                                           \
template<>                                   \
struct make_signed<unsigned unsigned_type> { \
    typedef signed unsigned_type type;       \
};
template<>
struct make_signed<char> {
    typedef signed char type;
};
GISMO_MAKE_SIGNED(char)
GISMO_MAKE_SIGNED(short)
GISMO_MAKE_SIGNED(int)
GISMO_MAKE_SIGNED(long)
GISMO_MAKE_SIGNED(long long)
#undef GISMO_MAKE_SIGNED

template<class T>
struct is_signed;
#define GISMO_IS_SIGNED(type)        \
template<>                           \
struct is_signed<signed type> {      \
    static const bool value = true;  \
};                                   \
template <>                          \
struct is_signed<unsigned type> {    \
    static const bool value = false; \
};
GISMO_IS_SIGNED(char)
GISMO_IS_SIGNED(short)
GISMO_IS_SIGNED(int)
GISMO_IS_SIGNED(long)
GISMO_IS_SIGNED(long long)
#undef GISMO_IS_SIGNED


#endif

template<typename T> struct remove_pointer {typedef T type;};
template<typename T> struct remove_pointer<T*> {typedef typename remove_pointer<T>::type type;};

template <class T> struct is_complex : public false_type {};
template <class T> struct is_complex<const T > : public is_complex<T>{};
template <class T> struct is_complex<volatile const T > : public is_complex<T>{};
template <class T> struct is_complex<volatile T > : public is_complex<T>{};
template <class T> struct is_complex<std::complex<T> > : public true_type{};

template <class T>
typename make_unsigned<T>::type to_unsigned(T t) {
    return t;
}
template <class T>
typename make_signed<T>::type to_signed(T t) {
    return t;
}

template<class T1, class T2>
bool less(T1 t1, T2 t2)
{
    typedef typename util::make_signed<T1>::type signedT1;
    typedef typename util::make_signed<T2>::type signedT2;
    typedef typename util::make_unsigned<T1>::type unsignedT1;
    typedef typename util::make_unsigned<T2>::type unsignedT2;

    if (is_signed<T1>::value == is_signed<T2>::value) // all is_signed are optimized out at compile time
    {
        if (is_signed<T1>::value) // both signed, cast to signedTx
            return (static_cast<signedT1>(t1) < static_cast<signedT2>(t2));
        else                      // both unsigned, cast to unsignedTx
            return (static_cast<unsignedT1>(t1) < static_cast<unsignedT2>(t2));
    }
    if (is_signed<T1>::value && !is_signed<T2>::value)
    {
        if (t1 < 0)
            return true;
        return (static_cast<unsignedT1>(t1) < static_cast<unsignedT2>(t2));
    }
    if (!is_signed<T1>::value && is_signed<T2>::value)
    {
        if (t2 < 0)
            return false;
        return (static_cast<unsignedT1>(t1) < static_cast<unsignedT2>(t2));
    }
}

template<class T1, class T2>
bool less_equal(T1 t1, T2 t2)
{
    typedef typename util::make_signed<T1>::type signedT1;
    typedef typename util::make_signed<T2>::type signedT2;
    typedef typename util::make_unsigned<T1>::type unsignedT1;
    typedef typename util::make_unsigned<T2>::type unsignedT2;

    if (is_signed<T1>::value == is_signed<T2>::value) // all is_signed are optimized out at compile time
    {
        if (is_signed<T1>::value) // both signed, cast to signedTx
            return (static_cast<signedT1>(t1) <= static_cast<signedT2>(t2));
        else                      // both unsigned, cast to unsignedTx
            return (static_cast<unsignedT1>(t1) <= static_cast<unsignedT2>(t2));
    }
    if (is_signed<T1>::value && !is_signed<T2>::value)
    {
        if (t1 < 0)
            return true;
        return (static_cast<unsignedT1>(t1) <= static_cast<unsignedT2>(t2));
    }
    if (!is_signed<T1>::value && is_signed<T2>::value)
    {
        if (t2 < 0)
            return false;
        return (static_cast<unsignedT1>(t1) <= static_cast<unsignedT2>(t2));
    }
}

template<class T1, class T2>
bool greater(T1 t1, T2 t2)
{
    return less(t2, t1);
}

template<class T1, class T2>
bool greater_equal(T1 t1, T2 t2)
{
    return less_equal(t2, t1);
}

template<class T1, class T2>
bool equal(T1 t1, T2 t2)
{
    typedef typename util::make_signed<T1>::type signedT1;
    typedef typename util::make_signed<T2>::type signedT2;
    typedef typename util::make_unsigned<T1>::type unsignedT1;
    typedef typename util::make_unsigned<T2>::type unsignedT2;

    if (is_signed<T1>::value == is_signed<T2>::value)
    {
        if (is_signed<T1>::value)
            return static_cast<signedT1>(t1) == static_cast<signedT2>(t2);
        else
            return static_cast<unsignedT1>(t1) == static_cast<unsignedT2>(t2);
    }
    if (is_signed<T1>::value && !is_signed<T2>::value)
    {
        if (t1 < 0)
            return false;
        return (static_cast<unsignedT1>(t1) == static_cast<unsignedT2>(t2));
    }
    if (!is_signed<T1>::value && is_signed<T2>::value)
    {
        if (t2 < 0)
            return false;
        return (static_cast<unsignedT1>(t1) == static_cast<unsignedT2>(t2));
    }
}

} // end namespace util

} // end namespace gismo