| ... | ... | @@ -909,6 +909,35 @@ test "fmin, fminf" { |
| 909 | 909 | } |
| 910 | 910 | } |
| 911 | 911 | |
| 912 | fn generic_fmax(comptime T: type, x: T, y: T) T { |
| 913 | if (isNan(x)) |
| 914 | return y; |
| 915 | if (isNan(y)) |
| 916 | return x; |
| 917 | return if (x < y) y else x; |
| 918 | } |
| 919 | |
| 920 | export fn fmaxf(x: f32, y: f32) callconv(.C) f32 { |
| 921 | return generic_fmax(f32, x, y); |
| 922 | } |
| 923 | |
| 924 | export fn fmax(x: f64, y: f64) callconv(.C) f64 { |
| 925 | return generic_fmax(f64, x, y); |
| 926 | } |
| 927 | |
| 928 | test "fmax, fmaxf" { |
| 929 | inline for ([_]type{ f32, f64 }) |T| { |
| 930 | const nan_val = math.nan(T); |
| 931 | |
| 932 | std.testing.expect(isNan(generic_fmax(T, nan_val, nan_val))); |
| 933 | std.testing.expectEqual(@as(T, 1.0), generic_fmax(T, nan_val, 1.0)); |
| 934 | std.testing.expectEqual(@as(T, 1.0), generic_fmax(T, 1.0, nan_val)); |
| 935 | |
| 936 | std.testing.expectEqual(@as(T, 10.0), generic_fmax(T, 1.0, 10.0)); |
| 937 | std.testing.expectEqual(@as(T, 1.0), generic_fmax(T, 1.0, -1.0)); |
| 938 | } |
| 939 | } |
| 940 | |
| 912 | 941 | // NOTE: The original code is full of implicit signed -> unsigned assumptions and u32 wraparound |
| 913 | 942 | // behaviour. Most intermediate i32 values are changed to u32 where appropriate but there are |
| 914 | 943 | // potentially some edge cases remaining that are not handled in the same way. |
LemonBoy <thatlemon@gmail.com> 2021-04-20 15:44:21+02:00
lib