summaryrefslogtreecommitdiffstats
path: root/libc-bottom-half/sources/math/fmin-fmax.c
blob: 92937984eb13d94a46339380a52ffe0cce473f58 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
// Wasm's `min` and `max` operators implement the IEEE 754-2019
// `minimum` and `maximum` operations, meaning that given a choice
// between NaN and a number, they return NaN. This differs from
// the C standard library's `fmin` and `fmax` functions, which
// return the number. However, we can still use wasm's builtins
// by handling the NaN cases explicitly, and it still turns out
// to be faster than doing the whole operation in
// target-independent C. And, it's smaller.

#include <math.h>

float fminf(float x, float y) {
    if (isnan(x)) return y;
    if (isnan(y)) return x;
    return __builtin_wasm_min_f32(x, y);
}

float fmaxf(float x, float y) {
    if (isnan(x)) return y;
    if (isnan(y)) return x;
    return __builtin_wasm_max_f32(x, y);
}

double fmin(double x, double y) {
    if (isnan(x)) return y;
    if (isnan(y)) return x;
    return __builtin_wasm_min_f64(x, y);
}

double fmax(double x, double y) {
    if (isnan(x)) return y;
    if (isnan(y)) return x;
    return __builtin_wasm_max_f64(x, y);
}