libm/math/generic/
fmaximum.rs

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