math.hpp

Go to the documentation of this file.

//          Copyright Dominic Koepke 2019 - 2023.
// Distributed under the Boost Software License, Version 1.0.
//    (See accompanying file LICENSE_1_0.txt or copy at
//          https://www.boost.org/LICENSE_1_0.txt)
 
#ifndef SL_UTILITY_MATH_HPP
#define SL_UTILITY_MATH_HPP
 
#pragma once
 
#include <cmath>
#include <concepts>
 
#if __cpp_lib_constexpr_cmath >= 202202L
#define SL_CONSTEXPR_MATH constexpr
#else
#define SL_CONSTEXPR_MATH
#endif
 
namespace sl::math
{
    template <std::floating_point T>
    struct remquo_result
    {
        T remainder;
        int quo;
    };
 
    template <class T>
    remquo_result(T, int) -> remquo_result<T>;
 
    constexpr auto remquo = [](const auto x, const auto y) SL_CONSTEXPR_MATH
            requires requires { std::remquo(x, x, nullptr); }
    {
        int quo{};
        const auto remainder{std::remquo(x, y, &quo)};
        return remquo_result{remainder, quo};
    };
 
    template <std::floating_point T>
    struct frexp_result
    {
        T fraction;
        int exp;
    };
 
    template <class T>
    frexp_result(T, int) -> frexp_result<T>;
 
    constexpr auto frexp = [](const auto value) SL_CONSTEXPR_MATH
            requires requires { std::frexp(value, nullptr); }
    {
        int exp{};
        const auto fraction{std::frexp(value, &exp)};
        return frexp_result{fraction, exp};
    };
 
    template <std::floating_point T>
    struct modf_result
    {
        T integral;
        T fraction;
    };
 
    constexpr auto modf = []<std::floating_point T>(const T value) SL_CONSTEXPR_MATH
    {
        T integral{};
        T fraction{std::modf(value, &integral)};
        return modf_result<T>{integral, fraction};
    };
 
}
 
#endif