add union of HyperRectangles

Author: maltezfaria

src/Interpolation/element.jl

@@ -75,6 +75,8 @@ range_dimension(el::AbstractElement{R,T}) where {R,T} = length(T)
 
 range_dimension(t::Type{<:AbstractElement{R,T}}) where {R,T} = length(T)
 
+center(el::AbstractElement{D}) where {D} = el(center(D()))
+
 """
     struct LagrangeElement{D,Np,T} <: AbstractElement{D,T}
 

src/Interpolation/hyperrectangle.jl

@@ -152,6 +152,12 @@ function Base.intersect(r1::HyperRectangle, r2::HyperRectangle)
     return HyperRectangle(lb, ub)
 end
 
+function Base.union(r1::HyperRectangle, r2::HyperRectangle)
+    lb = min.(low_corner(r1), low_corner(r2))
+    ub = max.(high_corner(r1), high_corner(r2))
+    return HyperRectangle(lb, ub)
+end
+
 Base.isempty(r::HyperRectangle) = any(low_corner(r) .> high_corner(r))
 
 ######