eigen_compat.hpp

Go to the documentation of this file.

#pragma once

 
#ifdef ENABLE_EIGEN_SUPPORT
#include <Eigen/Dense>
 
// Foward Declarations
namespace sigma {
template<typename T>
class Uncertain;
template<typename T>
class Interval;
} // namespace sigma

#define EIGEN_NUMTRAITS(float_type)                                          \
                   \
    template<>                                                               \
    struct NumTraits<sigma::Uncertain<float_type>> : NumTraits<float_type> { \
                                            \
        using Uncertain = sigma::Uncertain<float_type>;                      \
                                   \
        using Real = Uncertain;                                              \
                            \
        using NonInteger = Uncertain;                                        \
                                \
        using Literal = Uncertain;                                           \
                                 \
        using Nested = Uncertain;                                            \
        enum {                                                               \
            IsComplex             = 0,                                       \
            IsInteger             = 0,                                       \
            IsSigned              = 1,                                       \
            RequireInitialization = 1,                                       \
            ReadCost              = 1,                                       \
            AddCost               = 3,                                       \
            MulCost               = 3                                        \
        };                                                                   \
    }

#define EIGEN_NUMTRAITS_INTERVAL(float_type)                                \
                   \
    template<>                                                              \
    struct NumTraits<sigma::Interval<float_type>> : NumTraits<float_type> { \
                                            \
        using Interval = sigma::Interval<float_type>;                       \
                                  \
        using Real = Interval;                                              \
                           \
        using NonInteger = Interval;                                        \
                               \
        using Literal = Interval;                                           \
                                \
        using Nested = Interval;                                            \
        enum {                                                              \
            IsComplex             = 0,                                      \
            IsInteger             = 0,                                      \
            IsSigned              = 1,                                      \
            RequireInitialization = 1,                                      \
            ReadCost              = 1,                                      \
            AddCost               = 3,                                      \
            MulCost               = 3                                       \
        };                                                                  \
    }

namespace Eigen {
 
EIGEN_NUMTRAITS(float);
EIGEN_NUMTRAITS(double);
 
EIGEN_NUMTRAITS_INTERVAL(float);
EIGEN_NUMTRAITS_INTERVAL(double);
 
} // namespace Eigen
 
#undef EIGEN_NUMTRAITS
#undef EIGEN_NUMTRAITS_INTERVAL
#endif // ENABLE_EIGEN_SUPPORT