// Copyright (c) 2019-2024 Alexander Medvednikov. All rights reserved.
// Use of this source code is governed by an MIT license
// that can be found in the LICENSE file.
module math

union U32_F32 {
    u u32
    f f32
}

// f32_bits returns the IEEE 754 binary representation of f,
// with the sign bit of f and the result in the same bit position.
// f32_bits(f32_from_bits(x)) == x.
@[inline]
pub fn f32_bits(f f32) u32 {
    $if tinyc {
        return *unsafe { &u32(&f) } // this is faster for tcc, but causes `error: dereferencing type-punned pointer will break strict-aliasing rules` on gcc, with -cstrict
    }
    return unsafe {
        U32_F32{
            f: f
        }.u
    }
}

// f32_from_bits returns the floating-point number corresponding
// to the IEEE 754 binary representation b, with the sign bit of b
// and the result in the same bit position.
// f32_from_bits(f32_bits(x)) == x.
@[inline]
pub fn f32_from_bits(b u32) f32 {
    $if tinyc {
        return *unsafe { &f32(&b) }
    }
    return unsafe {
        U32_F32{
            u: b
        }.f
    }
}

union U64_F64 {
    u u64
    f f64
}

// f64_bits returns the IEEE 754 binary representation of f,
// with the sign bit of f and the result in the same bit position,
// and f64_bits(f64_from_bits(x)) == x.
@[inline]
pub fn f64_bits(f f64) u64 {
    $if tinyc {
        return *unsafe { &u64(&f) }
    }
    return unsafe {
        U64_F64{
            f: f
        }.u
    }
}

// f64_from_bits returns the floating-point number corresponding
// to the IEEE 754 binary representation b, with the sign bit of b
// and the result in the same bit position.
// f64_from_bits(f64_bits(x)) == x.
@[inline]
pub fn f64_from_bits(b u64) f64 {
    $if tinyc {
        return *unsafe { &f64(&b) }
    }
    return unsafe {
        U64_F64{
            u: b
        }.f
    }
}

// with_set_low_word sets low word of `f` to `lo`
@[inline]
pub fn with_set_low_word(f f64, lo u32) f64 {
    mut tmp := f64_bits(f)
    tmp &= 0xffffffff_00000000
    tmp |= u64(lo)
    return f64_from_bits(tmp)
}

// with_set_high_word sets high word of `f` to `lo`
@[inline]
pub fn with_set_high_word(f f64, hi u32) f64 {
    mut tmp := f64_bits(f)
    tmp &= 0x00000000_ffffffff
    tmp |= u64(hi) << 32
    return f64_from_bits(tmp)
}

// get_high_word returns high part of the word of `f`.
@[inline]
pub fn get_high_word(f f64) u32 {
    return u32(f64_bits(f) >> 32)
}