Fix union process

Author: juliohm

src/processes/union.jl

@@ -18,12 +18,12 @@ struct UnionSampling end
 
 default_sampling_algorithm(::UnionProcess) = UnionSampling()
 
-function rand_single(p::UnionProcess, ::Geometry, algo::UnionSampling)
-  pp₁ = rand(p.p₁, r)
-  pp₂ = rand(p.p₂, r)
+function rand_single(p::UnionProcess, g::Geometry, algo::UnionSampling)
+  pp₁ = rand(p.p₁, g)
+  pp₂ = rand(p.p₂, g)
 
-  X = coordinates(pp₁)
-  Y = coordinates(pp₂)
+  X = coordinates(pp₁, 1:nelements(pp₁))
+  Y = coordinates(pp₂, 1:nelements(pp₂))
 
   PointSet(hcat(X, Y))
 end

src/thinning.jl

@@ -14,14 +14,14 @@ abstract type AbstractThinning end
 
 Thin spatial point process `p` with thinning method `t`.
 """
-thin(p::PointProcess, t::AbstractThinning) = @error "not implemented"
+function thin end
 
 """
     thin(pp, t)
 
 Thin spatial point pattern `pp` with thinning method `t`.
 """
-thin(pp::PointSet, t::AbstractThinning) = @error "not implemented"
+thin(::PointSet, ::AbstractThinning)
 
 #-----------------
 # IMPLEMENTATIONS

test/processes.jl

@@ -16,4 +16,17 @@
   @testset "Poisson" begin
     # TODO
   end
+
+  @testset "Union" begin
+    b  = Box((0.,0.), (100.,100.))
+    p₁ = BinomialProcess(50)
+    p₂ = BinomialProcess(50)
+    p  = p₁ ∪ p₂ # 100 points
+
+    s = rand(p, b, 2)
+    @test length(s) == 2
+    @test s[1] isa PointSet
+    @test s[2] isa PointSet
+    @test nelements.(s) == [100,100]
+  end
 end