Implement boundstype

Author: OlivierHnt

src/decorations/intervals.jl

@@ -34,22 +34,7 @@ end
 DecoratedInterval(I::DecoratedInterval, dec::DECORATION) = DecoratedInterval(I.interval, dec)
 
 function DecoratedInterval(a::T, b::S, d::DECORATION) where {T<:Real, S<:Real}
-    BoundsType = promote_type(default_bound(), T, S)
-    is_valid_interval(a, b) || return DecoratedInterval(Interval{BoundsType}(a, b), ill)
-    return DecoratedInterval(Interval{BoundsType}(a, b), d)
-end
-function DecoratedInterval(a::T, b::S, d::DECORATION) where {T<:AbstractFloat, S<:Real}
-    BoundsType = promote_type(T, S)
-    is_valid_interval(a, b) || return DecoratedInterval(Interval{BoundsType}(a, b), ill)
-    return DecoratedInterval(Interval{BoundsType}(a, b), d)
-end
-function DecoratedInterval(a::T, b::S, d::DECORATION) where {T<:Real, S<:AbstractFloat}
-    BoundsType = promote_type(T, S)
-    is_valid_interval(a, b) || return DecoratedInterval(Interval{BoundsType}(a, b), ill)
-    return DecoratedInterval(Interval{BoundsType}(a, b), d)
-end
-function DecoratedInterval(a::T, b::S, d::DECORATION) where {T<:AbstractFloat, S<:AbstractFloat}
-    BoundsType = promote_type(T, S)
+    BoundsType = boundstype(T, S)
     is_valid_interval(a, b) || return DecoratedInterval(Interval{BoundsType}(a, b), ill)
     return DecoratedInterval(Interval{BoundsType}(a, b), d)
 end
@@ -67,22 +52,7 @@ end
 DecoratedInterval(I::Interval) = DecoratedInterval{default_bound()}(I)
 
 function DecoratedInterval(a::T, b::S) where {T<:Real, S<:Real}
-    BoundsType = promote_type(default_bound(), T, S)
-    is_valid_interval(a, b) || return DecoratedInterval(Interval{BoundsType}(a, b), ill)
-    return DecoratedInterval(Interval{BoundsType}(a, b))
-end
-function DecoratedInterval(a::T, b::S) where {T<:AbstractFloat, S<:Real}
-    BoundsType = promote_type(T, S)
-    is_valid_interval(a, b) || return DecoratedInterval(Interval{BoundsType}(a, b), ill)
-    return DecoratedInterval(Interval{BoundsType}(a, b))
-end
-function DecoratedInterval(a::T, b::S) where {T<:Real, S<:AbstractFloat}
-    BoundsType = promote_type(T, S)
-    is_valid_interval(a, b) || return DecoratedInterval(Interval{BoundsType}(a, b), ill)
-    return DecoratedInterval(Interval{BoundsType}(a, b))
-end
-function DecoratedInterval(a::T, b::S) where {T<:AbstractFloat, S<:AbstractFloat}
-    BoundsType = promote_type(T, S)
+    BoundsType = boundstype(T, S)
     is_valid_interval(a, b) || return DecoratedInterval(Interval{BoundsType}(a, b), ill)
     return DecoratedInterval(Interval{BoundsType}(a, b))
 end

src/intervals/construction.jl

@@ -34,6 +34,12 @@ Interval{Float64}
 """
 default_bound() = Float64
 
+# Produce the type of the bounds of an interval when not explicitly imposed
+boundstype(::Type{T}, ::Type{S}) where {T<:AbstractFloat, S<:AbstractFloat} = promote_type(T, S)
+boundstype(::Type{T}, ::Type{S}) where {T<:AbstractFloat, S} = promote_type(T, S)
+boundstype(::Type{T}, ::Type{S}) where {T, S<:AbstractFloat} = promote_type(T, S)
+boundstype(::Type{T}, ::Type{S}) where {T, S} = promote_type(default_bound(), T, S)
+
 @inline _normalisezero(a::Real) = ifelse(iszero(a), zero(a), a)
 
 """
@@ -88,11 +94,7 @@ function interval(::Type{T}, a, b) where {T}
     @warn "Invalid input, empty interval is returned"
     return emptyinterval(T)
 end
-interval(a::T, b::S) where {T<:AbstractFloat, S<:AbstractFloat} =
-    interval(promote_type(T, S), a, b)
-interval(a::T, b::S) where {T<:AbstractFloat, S} = interval(promote_type(T, S), a, b)
-interval(a::T, b::S) where {T, S<:AbstractFloat} = interval(promote_type(T, S), a, b)
-interval(a::T, b::S) where {T, S} = interval(promote_type(default_bound(), T, S), a, b)
+interval(a::T, b::S) where {T, S} = interval(boundstype(T, S), a, b)
 
 # Real: `is_valid_interval(a, a) != true`
 interval(::Type{T}, a::Real) where {T} = interval(T, a, a)