math.h

Go to the documentation of this file.

#pragma once
#include "interop/config.h"
#include <limits>

#include <cmath>
#if defined(HAVE_NO_STD_ISNAN)
#   if defined(HAVE_ISNAN)
#       include <math.h>
    #elif defined(HAVE___ISNAN)
#       include <math.h>
    #elif defined(HAVE_FLOAT_H_ISNAN)
#       include <float.h>
    #endif
    namespace std
    {
        template<typename T>
        bool isnan(const T val)
        {
#   if defined(HAVE_ISNAN)
            return ::isnan(val);
#   elif defined(HAVE___ISNAN)
            return __isnan(val);
#   elif defined(HAVE_FLOAT_H_ISNAN)
            return _isnan(val);
#   else
            return val != val;
#endif
        }
        template<typename T>
        bool isinf(const T val)
        {
#   if defined(HAVE_ISNAN)
            return ::isinf(val);
#   elif defined(HAVE___ISNAN)
            return __isinf(val);
#   elif defined(HAVE_FLOAT_H_ISNAN)
            return _isinf(val);
#   else
            return !isnan(val) && isnan(val-val);
#endif
        }
    }
#endif

namespace illumina { namespace interop { namespace util
{
    template<typename T>
    float float_cast(const T val, const T sentinel = std::numeric_limits<T>::max())
    {
        if (val == sentinel)
            return std::numeric_limits<float>::quiet_NaN();
        return static_cast<float>(val);
    }
}}}