fix translation map of vcatdomain

JuliaApproximation · Mar 3, 2024 · 343812a · 343812a · daanhb · Mar 3, 2024
1 parent 8742d3a
commit 343812a
Show file tree

Hide file tree

Showing 3 changed files with 7 additions and 3 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,6 +1,6 @@
 name = "DomainSets"
 uuid = "5b8099bc-c8ec-5219-889f-1d9e522a28bf"
-version = "0.7.7"
+version = "0.7.8"
 
 [deps]
 CompositeTypes = "b152e2b5-7a66-4b01-a709-34e65c35f657"

diff --git a/src/generic/productdomain.jl b/src/generic/productdomain.jl
@@ -121,8 +121,9 @@ function map_domain(linmap::GenericLinearMap{SVector{N,S},A}, domain::ProductDom
 end
 
 function map_domain(transmap::Translation{SVector{N,S}}, domain::ProductDomain{SVector{N,T}}) where {N,S,T}
-	vec = unsafe_vector(transmap)
-	ProductDomain{SVector{N,promote_type(S,T)}}(map( (d,v) -> d .+ v, components(domain), vec))
+    vec = unsafe_vector(transmap)
+    ProductDomain{SVector{N,promote_type(S,T)}}(
+            map( (d,v) -> d .+ v, components(domain), tointernalpoint(domain, vec)))
 end
 
 

diff --git a/test/test_domain_product.jl b/test/test_domain_product.jl
@@ -53,6 +53,9 @@ function test_product_domains()
         @test SA[0.5,0.5,0.8] ∈ d3
         @test SA[-1.1,0.3,0.1] ∉ d3
         @test choice(d3) ∈ d3
+
+        v = SA[0.1,0.2,0.3]
+        @test map_domain(Translation(v), productdomain(UnitDisk(), UnitInterval())) isa VcatDomain
     end
     @testset "mixed intervals" begin
         d = (0..1) × (0.0..1)