Vectornaut/dyna3
1
0
Fork 0
forked from StudioInfinity/dyna3
Code Pull requests Activity

Make results reproducible

Browse source
This commit is contained in:
Aaron Fenyes 2024-02-15 16:00:46 -08:00
parent 8d8bc9162c
commit ae5db0f9ea
2 changed files with 7 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

6
engine-proto/Engine.jl
View file

@ -89,11 +89,11 @@ vbls = Variable.(symbols(coordring))
# test a random witness set
system = CompiledSystem(System(small_eqns_tan_sph, variables = vbls))
norm2 = vec -> real(dot(conj.(vec), vec))
Random.seed!(6071)
n_planes = 36
rng = MersenneTwister(6071)
n_planes = 3
samples = []
for _ in 1:n_planes
real_solns = solution.(Engine.Numerical.real_samples(system, freedom))
real_solns = solution.(Engine.Numerical.real_samples(system, freedom, rng = rng))
for soln in real_solns
if all(norm2(soln - samp) > 1e-4*length(gens(coordring)) for samp in samples)
push!(samples, soln)