Add missing constructor, fix typo

src/specialitensornetworks.jl

@@ -70,14 +70,26 @@ Build an ITensor network on a graph specified by the inds network s.
 Bond_dim is given by link_space and entries are randomized.
 The random distribution is based on the input argument `distribution`.
 """
-function random_tensornetwork(distribution::Distribution, s::IndsNetwork; kwargs...)
+function random_tensornetwork(
+  rng::AbstractRNG, distribution::Distribution, s::IndsNetwork; kwargs...
+)
   return ITensorNetwork(s; kwargs...) do v
-    return inds -> itensor(rand(distribution, dim.(inds)...), inds)
+    return inds -> itensor(rand(rng, distribution, dim.(inds)...), inds)
   end
 end
 
+function random_tensornetwork(distribution::Distribution, s::IndsNetwork; kwargs...)
+  return random_tensornetwork(Random.default_rng(), distribution, s; kwargs...)
+end
+
+@traitfn function random_tensornetwork(
+  rng::AbstractRNG, distribution::Distribution, g::::IsUnderlyingGraph; kwargs...
+)
+  return random_tensornetwork(rng, distribution, IndsNetwork(g); kwargs...)
+end
+
 @traitfn function random_tensornetwork(
   distribution::Distribution, g::::IsUnderlyingGraph; kwargs...
 )
-  return random_tensornetwork(distribution, IndsNetwork(g); kwargs...)
+  return random_tensornetwork(Random.default_rng(), distribution, g; kwargs...)
 end

test/test_itensornetwork.jl

@@ -349,7 +349,7 @@ const elts = (Float32, Float64, Complex{Float32}, Complex{Float64})
     dims = (2, 2)
     g = named_grid(dims)
     s = siteinds("S=1/2", g)
-    rng = StableRN(1234)
+    rng = StableRNG(1234)
     ψ = random_tensornetwork(rng, s; link_space=2)
     @test scalartype(ψ) == Float64
     ϕ = NDTensors.convert_scalartype(new_eltype, ψ)