Implement rand()

Project.toml

@@ -6,6 +6,7 @@ version = "0.5.12"
 CompositeTypes = "b152e2b5-7a66-4b01-a709-34e65c35f657"
 IntervalSets = "8197267c-284f-5f27-9208-e0e47529a953"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
+Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 

src/DomainSets.jl

@@ -3,6 +3,8 @@ module DomainSets
 using StaticArrays
 using LinearAlgebra, Statistics
 import LinearAlgebra: cross, ×, pinv
+import Random
+using Random: AbstractRNG
 
 using IntervalSets
 using CompositeTypes, CompositeTypes.Display, CompositeTypes.Indexing
@@ -216,6 +218,7 @@ include("domains/indicator.jl")
 include("domains/boundingbox.jl")
 
 include("applications/coordinates.jl")
+include("applications/random.jl")
 include("applications/rotation.jl")
 
 end # module

src/applications/random.jl

@@ -0,0 +1,18 @@
+Random.gentype(::Type{<:Domain{T}}) where T = T
+
+# Random.gentype(::Type{<:ProductDomain{T}}) where T = T
+Base.rand(rng::AbstractRNG, s::Random.SamplerTrivial{<:ProductDomain}) = map(i->(rand(rng, i)), factors(s[]))
+
+# Random.gentype(::Type{DerivedDomain{T}}) where T = Random.gentype(D)
+Base.rand(rng::AbstractRNG, s::Random.SamplerTrivial{<:SimpleLazyDomain}) = rand(rng, superdomain(s[]))
+
+# Ball
+# Mapped
+
+
+# Rectangle
+# Wrapped
+# Product Domain
+
+# # Maybe
+# Sphere

test/runtests.jl

@@ -1,13 +1,13 @@
-using Test, LinearAlgebra, StaticArrays
+using Test, LinearAlgebra, StaticArrays, Random
 
 using DomainSets
 using CompositeTypes.Indexing
 
-include("test_common.jl")
+# include("test_common.jl")
 include("test_maps.jl")
-include("test_generic_domain.jl")
-include("test_specific_domains.jl")
-include("test_canonical.jl")
-include("test_setoperations.jl")
+# include("test_generic_domain.jl")
+# include("test_specific_domains.jl")
+# include("test_canonical.jl")
+# include("test_setoperations.jl")
 include("test_applications.jl")
-include("test_readme.jl")
+# include("test_readme.jl")

test/test_applications.jl

@@ -48,11 +48,32 @@ function test_cart_polar_map(T)
     @test inverse(m2) == m1
 end
 
+function test_rand(T)
+    r = Rectangle(T[-1, 2], T[3, 4])
+    @test @inferred(Random.gentype(r)) == Vector{T}
+    @test typeof(rand(r)) == Random.gentype(r)
+    @test @inferred(rand(r)) in r
 
+    r = Rectangle(SA[T(-1), T(2)], SA[T(3), T(4)])
+    @test @inferred(Random.gentype(r)) == SVector{2, T}
+    @test typeof(rand(r)) == Random.gentype(r)
+    @test @inferred(rand(r)) in r
+
+    hybrid_product = ProductDomain(["a", "b"], 1.0..2.0)
+    @test @inferred(Random.gentype(hybrid_product)) == Tuple{String, Float64}
+    @test typeof(rand(hybrid_product)) == Random.gentype(hybrid_product)
+    @test @inferred(rand(hybrid_product)) in hybrid_product
+
+    # b = Ball(2.0, SA[1.0, 2.0])
+    # @test @inferred(Random.gentype(r)) == SVector{2, T}
+    # @test typeof(rand(r)) == Random.gentype(r)
+    # @test @inferred(rand(r)) in r
+end
 
 function test_applications(T)
     test_rotation_map(T)
     test_cart_polar_map(T)
+    test_rand(T)
 end
 
 @testset "applications" begin