libm/math/generic/
fminimum_num.rs

1/* SPDX-License-Identifier: MIT OR Apache-2.0 */
2//! IEEE 754-2019 `minimum`.
3//!
4//! Per the spec, returns:
5//! - `x` if `x < y`
6//! - `y` if `y < x`
7//! - Non-NaN if one operand is NaN
8//! - Logic following +0.0 > -0.0
9//! - Either `x` or `y` if `x == y` and the signs are the same
10//! - qNaN if either operand is a NaN
11//!
12//! Excluded from our implementation is sNaN handling.
13
14use crate::support::Float;
15
16#[inline]
17pub fn fminimum_num<F: Float>(x: F, y: F) -> F {
18    let res =
19        if y.is_nan() || x < y || (x.to_bits() == F::NEG_ZERO.to_bits() && y.is_sign_positive()) {
20            x
21        } else {
22            y
23        };
24
25    // Canonicalize
26    res * F::ONE
27}