Remove _extrema_rf optimization for AbstractFloat (follow advise)

Add TODO

Co-Authored-By: Takafumi Arakaki <29282+tkf@users.noreply.github.com>

Author: N5N3

base/mpfr.jl

@@ -706,17 +706,6 @@ function min(x::BigFloat, y::BigFloat)
     return z
 end
 
-# We don't want allocations during `extrema`.
-# TODO: avoid allocations in `min/max`
-function Base._extrema_rf(x::NTuple{2,BigFloat}, y::NTuple{2,BigFloat})
-    (x1, x2), (y1, y2) = x, y
-    isnan(x1) && return x
-    isnan(y1) && return y
-    z1 = x1 < y1 || signbit(x1) > signbit(y1) ? x1 : y1
-    z2 = x2 < y2 || signbit(x2) > signbit(y2) ? y2 : x2
-    z1, z2
-end
-
 function modf(x::BigFloat)
     zint = BigFloat()
     zfloat = BigFloat()

base/reduce.jl

@@ -855,15 +855,8 @@ end
 ExtremaMap(::Type{T}) where {T} = ExtremaMap{Type{T}}(T)
 @inline (f::ExtremaMap)(x) = (y = f.f(x); (y, y))
 
+# TODO: optimize for inputs <: AbstractFloat
 @inline _extrema_rf((min1, max1), (min2, max2)) = (min(min1, min2), max(max1, max2))
-# optimization for AbstractFloat
-function _extrema_rf(x::NTuple{2,T}, y::NTuple{2,T}) where {T<:AbstractFloat}
-    (x1, x2), (y1, y2) = x, y
-    anynan = isnan(x1)|isnan(y1)
-    z1 = ifelse(anynan, x1-y1, ifelse(signbit(x1-y1), x1, y1))
-    z2 = ifelse(anynan, x1-y1, ifelse(signbit(x2-y2), y2, x2))
-    z1, z2
-end
 
 ## findmax, findmin, argmax & argmin