1#[cfg_attr(all(test, assert_no_panic), no_panic::no_panic)]
2pub fn fmax(x: f64, y: f64) -> f64 {
3 // IEEE754 says: maxNum(x, y) is the canonicalized number y if x < y, x if y < x, the
4 // canonicalized number if one operand is a number and the other a quiet NaN. Otherwise it
5 // is either x or y, canonicalized (this means results might differ among implementations).
6 // When either x or y is a signalingNaN, then the result is according to 6.2.
7 //
8 // Since we do not support sNaN in Rust yet, we do not need to handle them.
9 // FIXME(nagisa): due to https://bugs.llvm.org/show_bug.cgi?id=33303 we canonicalize by
10 // multiplying by 1.0. Should switch to the `canonicalize` when it works.
11 (if x.is_nan() || x < y { y } else { x }) * 1.0
12}
13