| 1 | /* SPDX-License-Identifier: MIT OR Apache-2.0 */ |
| 2 | //! IEEE 754-2019 `maximumNumber`. |
|---|
| 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 | |
|---|
| 14 | use crate::support::Float; |
|---|
| 15 | |
|---|
| 16 | #[inline ] |
|---|
| 17 | pub fn fmaximum_num<F: Float>(x: F, y: F) -> F { |
|---|
| 18 | let res: F = if x.is_nan() || x < y || (x.biteq(F::NEG_ZERO) && y.is_sign_positive()) { |
|---|
| 19 | y |
|---|
| 20 | } else { |
|---|
| 21 | x |
|---|
| 22 | }; |
|---|
| 23 | |
|---|
| 24 | // Canonicalize |
|---|
| 25 | res * F::ONE |
|---|
| 26 | } |
|---|
| 27 | |
|---|
