inline real_concept modf(real_concept a, real_concept* ipart)
{
   long double ip;
   long double result = std::modf(a.value(), &ip);
   *ipart = ip;
   return result;
}

inline real_concept modf(real_concept a, real_concept* ipart)
{
   real_concept_base_type ip;
   real_concept_base_type result = std::modf(a.value(), &ip);
   *ipart = ip;
   return result;
}

inline real_concept log10(real_concept a)
{ return std::log10(a.value()); }

inline real_concept ldexp(real_concept a, int expon)
{ return std::ldexp(a.value(), expon); }

inline real_concept frexp(real_concept a, int* expon)
{ return std::frexp(a.value(), expon); }

inline real_concept pow(real_concept a, int b)
{ return std::pow(a.value(), static_cast<real_concept_base_type>(b)); }

inline real_concept sqrt(real_concept a)
{ return std::sqrt(a.value()); }

inline real_concept fmod(real_concept a, real_concept b)
{ return fmodl(a.value(), b.value()); }

inline real_concept cosh(real_concept a)
{ return std::cosh(a.value()); }

inline real_concept atan2(real_concept a, real_concept b)
{ return std::atan2(a.value(), b.value()); }

inline real_concept ceil(real_concept a)
{ return std::ceil(a.value()); }

inline real_concept atan(real_concept a)
{ return std::atan(a.value()); }

inline bool operator >= (const real_concept& a, const real_concept& b)
{ return a.value() >= b.value(); }

inline real_concept pow(real_concept a, int b)
{ return std::pow(a.value(), static_cast<long double>(b)); }

inline real_concept exp(real_concept a)
{ return std::exp(a.value()); }

inline real_concept abs(real_concept a)
{ return std::abs(a.value()); }

inline real_concept sinh(real_concept a)
{ return std::sinh(a.value()); }

inline real_concept floor(real_concept a)
{ return std::floor(a.value()); }

inline real_concept tanh(real_concept a)
{ return std::tanh(a.value()); }

inline bool operator <= (const T& a, const real_concept& b)
{
   return a <= b.value();
}