add rng argument to octorand function

JuliaGeometry · Jan 24, 2023 · ba18e61 · ba18e61
1 parent ae650c6
commit ba18e61
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/octonion.jl b/src/octonion.jl
@@ -130,7 +130,7 @@ function Base.sqrt(o::Octonion)
   exp(0.5 * log(o))
 end
 
-octorand() = octo(randn(), randn(), randn(), randn(), randn(), randn(), randn(), randn())
+octorand(rng::AbstractRNG = Random.GLOBAL_RNG) = octo(randn(rng), randn(rng), randn(rng), randn(rng), randn(rng), randn(rng), randn(rng), randn(rng))
 
 function Base.rand(rng::AbstractRNG, ::Random.SamplerType{Octonion{T}}) where {T<:Real}
   Octonion{T}(rand(rng, T), rand(rng, T), rand(rng, T), rand(rng, T),