Add InhibitionProcess

Author: juliohm

src/PointPatterns.jl

@@ -19,6 +19,7 @@ export
   PointProcess,
   BinomialProcess,
   PoissonProcess,
+  InhibitionProcess,
   ClusterProcess,
   UnionProcess,
   ishomogeneous,

src/processes.jl

@@ -45,5 +45,6 @@ function randsingle end
 
 include("processes/binomial.jl")
 include("processes/poisson.jl")
+include("processes/inhibition.jl")
 include("processes/cluster.jl")
 include("processes/union.jl")

src/processes/inhibition.jl

@@ -0,0 +1,14 @@
+# ------------------------------------------------------------------
+# Licensed under the MIT License. See LICENSE in the project root.
+# ------------------------------------------------------------------
+
+"""
+    InhibitionProcess(r)
+
+An inhibition point process with radius `r`.
+"""
+struct InhibitionProcess{R<:Real} <: PointProcess
+  r::R
+end
+
+randsingle(rng::Random.AbstractRNG, p::InhibitionProcess, g) = PointSet(sample(rng, g, MinDistanceSampling(p.r)))

test/processes.jl

@@ -20,7 +20,8 @@
   binom = BinomialProcess(100)
   poisson1 = PoissonProcess(100.0)
   poisson2 = PoissonProcess(λ)
-  procs = [binom, poisson1, poisson2]
+  inhibit  = InhibitionProcess(0.1)
+  procs = [binom, poisson1, poisson2, inhibit]
 
   @testset "Basic" begin
     for p in procs, g in geoms
@@ -54,6 +55,12 @@
     @test isnothing(rand(PoissonProcess(100.0), pp))
   end
 
+  @testset "Inhibition" begin
+  end
+
+  @testset "Cluster" begin
+  end
+
   @testset "Union" begin
     b = Box((0.0, 0.0), (100.0, 100.0))
     p₁ = BinomialProcess(50)