Fix mixed sign case

As it turns out, our `mod` is not the same as C's `fmod`; that one is `rem`.
JuliaLang · Sep 25, 2021 · 9d4df8f · 9d4df8f
1 parent 0234844
commit 9d4df8f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/base/math.jl b/base/math.jl
@@ -893,7 +893,7 @@ function rem(x::Float64, p::Float64, ::RoundingMode{:Nearest})
         return NaN
     end
     if hp <= 0x7fdfffff
-        x = mod(x, p + p)  # now x < 2p
+        x = rem(x, p + p)  # now x < 2p
     end
     ((hx - hp) | (lx - lp)) == 0 && return 0.0
     x = abs(x)
@@ -932,7 +932,7 @@ function rem(x::Float32, p::Float32, ::RoundingMode{:Nearest})
         return NaN32
     end
     if hp <= 0x7effffff
-        x = mod(x, p + p)  # now x < 2p
+        x = rem(x, p + p)  # now x < 2p
     end
     hx - hp == 0 && return 0.0f0
     x = abs(x)