1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16

use core::f64;

#[inline]
pub fn fdim(x: f64, y: f64) -> f64 {
    if x.is_nan() {
        x
    } else if y.is_nan() {
        y
    } else {
        if x > y {
            x - y
        } else {
            0.0
        }
    }
}