//
// Copyright (C) 2024-2026 Intel Corporation
// SPDX-License-Identifier: Apache-2.0
//

#pragma once

#include <cmath>
#include <cstdint>
#include <limits>
#include <ostream>

namespace vpux {
namespace type {
/**
 * @brief Class to represent the f8e4m3 type.
 */
class float8_e5m2 {
public:
    float8_e5m2() = default;
    float8_e5m2(uint32_t sign, uint32_t biased_exponent, uint32_t fraction);
    float8_e5m2(float value);

    template <typename I>
    explicit float8_e5m2(I value): m_value{float8_e5m2{static_cast<float>(value)}.m_value} {
    }

    size_t size() const;
    template <typename T>
    bool operator==(const T& other) const;
    template <typename T>
    bool operator!=(const T& other) const {
        return !(*this == other);
    }

    template <typename T>
    bool operator<(const T& other) const;
    template <typename T>
    bool operator<=(const T& other) const;
    template <typename T>
    bool operator>(const T& other) const;
    template <typename T>
    bool operator>=(const T& other) const;
    template <typename T>
    float8_e5m2 operator+(const T& other) const;
    template <typename T>
    float8_e5m2 operator+=(const T& other);
    template <typename T>
    float8_e5m2 operator-(const T& other) const;
    template <typename T>
    float8_e5m2 operator-=(const T& other);
    template <typename T>
    float8_e5m2 operator*(const T& other) const;
    template <typename T>
    float8_e5m2 operator*=(const T& other);
    template <typename T>
    float8_e5m2 operator/(const T& other) const;
    template <typename T>
    float8_e5m2 operator/=(const T& other);

    operator float() const;

    static constexpr float8_e5m2 from_bits(uint8_t bits) {
        return float8_e5m2(bits, true);
    }
    uint8_t to_bits() const;
    friend std::ostream& operator<<(std::ostream& out, const float8_e5m2& obj) {
        out << static_cast<float>(obj);
        return out;
    }

private:
    constexpr float8_e5m2(uint8_t x, bool): m_value{x} {
    }

    uint8_t m_value;
};

#if defined(_MSC_VER)
#pragma warning(push)
#pragma warning(disable : 4756)
#endif
template <typename T>
bool float8_e5m2::operator==(const T& other) const {
#if defined(__GNUC__)
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wfloat-equal"
#endif
    return (static_cast<float>(*this) == static_cast<float>(other));
#if defined(__GNUC__)
#pragma GCC diagnostic pop
#endif
}

template <typename T>
bool float8_e5m2::operator<(const T& other) const {
    return (static_cast<float>(*this) < static_cast<float>(other));
}

template <typename T>
bool float8_e5m2::operator<=(const T& other) const {
    return (static_cast<float>(*this) <= static_cast<float>(other));
}

template <typename T>
bool float8_e5m2::operator>(const T& other) const {
    return (static_cast<float>(*this) > static_cast<float>(other));
}

template <typename T>
bool float8_e5m2::operator>=(const T& other) const {
    return (static_cast<float>(*this) >= static_cast<float>(other));
}

template <typename T>
float8_e5m2 float8_e5m2::operator+(const T& other) const {
    return {static_cast<float>(*this) + static_cast<float>(other)};
}

template <typename T>
float8_e5m2 float8_e5m2::operator+=(const T& other) {
    return *this = *this + other;
}

template <typename T>
float8_e5m2 float8_e5m2::operator-(const T& other) const {
    return {static_cast<float>(*this) - static_cast<float>(other)};
}

template <typename T>
float8_e5m2 float8_e5m2::operator-=(const T& other) {
    return *this = *this - other;
}

template <typename T>
float8_e5m2 float8_e5m2::operator*(const T& other) const {
    return {static_cast<float>(*this) * static_cast<float>(other)};
}

template <typename T>
float8_e5m2 float8_e5m2::operator*=(const T& other) {
    return *this = *this * other;
}

template <typename T>
float8_e5m2 float8_e5m2::operator/(const T& other) const {
    return {static_cast<float>(*this) / static_cast<float>(other)};
}

template <typename T>
float8_e5m2 float8_e5m2::operator/=(const T& other) {
    return *this = *this / other;
}
#if defined(_MSC_VER)
#pragma warning(pop)
#endif

}  // namespace type
}  // namespace vpux

namespace std {
template <>
class numeric_limits<vpux::type::float8_e5m2> {
public:
    static constexpr bool is_specialized = true;
    static constexpr vpux::type::float8_e5m2 min() noexcept {
        return vpux::type::float8_e5m2::from_bits(0b00000100);  // minimum positive normalized value
    }
    static constexpr vpux::type::float8_e5m2 max() noexcept {
        return vpux::type::float8_e5m2::from_bits(0b01111011);
    }
    static constexpr vpux::type::float8_e5m2 lowest() noexcept {
        return vpux::type::float8_e5m2::from_bits(0b11111011);
    }
    static constexpr int digits = 3;
    static constexpr int digits10 = 0;

    static constexpr bool is_signed = true;
    static constexpr bool is_integer = false;
    static constexpr bool is_exact = false;

    static constexpr int radix = 2;

    static constexpr vpux::type::float8_e5m2 epsilon() noexcept {
        return vpux::type::float8_e5m2::from_bits(0b00000001);
    }
    static constexpr vpux::type::float8_e5m2 round_error() noexcept {
        return vpux::type::float8_e5m2::from_bits(0b00111000);
    }

    static constexpr int min_exponent = -13;
    static constexpr int min_exponent10 = -4;
    static constexpr int max_exponent = 16;
    static constexpr int max_exponent10 = 4;

    static constexpr bool has_infinity = true;
    static constexpr bool has_quiet_NaN = true;
    static constexpr bool has_signaling_NaN = true;

    static constexpr float_denorm_style has_denorm = denorm_present;
    static constexpr bool has_denorm_loss = false;

    static constexpr vpux::type::float8_e5m2 infinity() noexcept {
        return vpux::type::float8_e5m2::from_bits(0b01111100);
    }
    static constexpr vpux::type::float8_e5m2 quiet_NaN() noexcept {
        return vpux::type::float8_e5m2::from_bits(0b01111111);
    }
    static constexpr vpux::type::float8_e5m2 signaling_NaN() noexcept {
        return vpux::type::float8_e5m2::from_bits(0b01111101);
    }
    static constexpr vpux::type::float8_e5m2 denorm_min() noexcept {
        return vpux::type::float8_e5m2::from_bits(0b00000001);  // minimum positive denormalized value
    }
    static constexpr bool is_iec559 = false;
    static constexpr bool is_bounded = true;
    static constexpr bool is_modulo = false;
    static constexpr bool traps = false;
    static constexpr bool tinyness_before = false;
    static constexpr float_round_style round_style = round_to_nearest;
};
}  // namespace std
