precision.hpp

 
 
#ifndef mx_precision_hpp
#define mx_precision_hpp
 
#include <type_traits>
 
#ifdef MX_INCLUDE_BOOST
#include <boost/math/tools/precision.hpp>
#endif
 
namespace mx
{
namespace math
{
namespace func
{
 
/// Get the sqrt(epsilon) where epsilon is machine precision
/** Wrapper for boost function
 *
 * \todo this isn't, but should be constexpr -- need constr sqrt.
 *
 * \ingroup gen_math_precision
 */
template <typename T>
T root_epsilon()
{
#ifdef MX_INCLUDE_BOOST
    return boost::math::tools::root_epsilon<T>();
#else
    static_assert(
        std::is_fundamental<T>::value || !std::is_fundamental<T>::value,
        "root_epsilon<T> not specialized for type T, and MX_INCLUDE_BOOST is not defined, so I can't just use boost." );
    return 0;
#endif
}
 
template <>
float root_epsilon<float>();
 
template <>
double root_epsilon<double>();
 
} // namespace func
} // namespace math
} // namespace mx
 
#endif // mx_precision_hpp