math_utils.h

Go to the documentation of this file.

#ifndef MATH_UTILS_H
#define MATH_UTILS_H
 
#include "config.h"             // For PRECISION_TOLERANCE and MyErrorHandler
#include "simple_type_traits.h" // For is_floating_point_v and is_same_v
 
namespace math
{
 
    template <typename T>
    constexpr T sqrt(T x) noexcept
    {
        static_assert(is_floating_point_v<T>, "sqrt requires a floating-point type");
        if (x < T(0))
            MyErrorHandler::error("sqrt of negative value");
        return is_same_v<T, float> ? __builtin_sqrtf(x) : __builtin_sqrt(x);
    }
 
    template <typename T>
    constexpr T abs(T x) noexcept
    {
        if constexpr (is_same_v<T, float>)
            return __builtin_fabsf(x);
        else if constexpr (is_same_v<T, double>)
            return __builtin_fabs(x);
        else if constexpr (is_same_v<T, int>)
            return __builtin_abs(x);
        else if constexpr (is_same_v<T, long>)
            return __builtin_labs(x);
        else
            return x < T(0) ? -x : x;
    }
 
    template <typename T>
    constexpr T sin(T x) noexcept
    {
        static_assert(is_floating_point_v<T>, "sin requires a floating-point type");
        return is_same_v<T, float> ? __builtin_sinf(x) : __builtin_sin(x);
    }
 
    template <typename T>
    constexpr T cos(T x) noexcept
    {
        static_assert(is_floating_point_v<T>, "cos requires a floating-point type");
        return is_same_v<T, float> ? __builtin_cosf(x) : __builtin_cos(x);
    }
 
    template <typename T>
    constexpr T acos(T x) noexcept
    {
        static_assert(is_floating_point_v<T>, "acos requires a floating-point type");
        return is_same_v<T, float> ? __builtin_acosf(x) : __builtin_acos(x);
    }
 
} // namespace math
 
#endif // MATH_UTILS_H