glen/dyna3
1
0
Fork 0
Code Issues 11 Pull Requests 1 Actions Projects Releases 2 Wiki Activity

Test different sign patterns for cosines

Browse Source 
It seems like there are real solutions if and only if the product of the
cosines is positive.
This commit is contained in:
Aaron Fenyes 2024-02-15 16:25:09 -08:00
parent ba365174d3
commit f2000e5731
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
engine-proto/Engine.jl
View File

@ -44,7 +44,7 @@ tangencies = [
Engine.AlignsWithBy{CoeffType}( Engine.AlignsWithBy{CoeffType}(
spheres[n], spheres[n],
spheres[mod1(n+1, length(spheres))], spheres[mod1(n+1, length(spheres))],
CoeffType(-1)^n CoeffType([1, 1, 1][n])
) )
for n in 1:3 for n in 1:3
] ]
@ -90,7 +90,7 @@ vbls = Variable.(symbols(coordring))
system = CompiledSystem(System(eqns_tan_sph, variables = vbls)) system = CompiledSystem(System(eqns_tan_sph, variables = vbls))
norm2 = vec -> real(dot(conj.(vec), vec)) norm2 = vec -> real(dot(conj.(vec), vec))
rng = MersenneTwister(6071) rng = MersenneTwister(6071)
n_planes = 3 n_planes = 36
samples = [] samples = []
for _ in 1:n_planes for _ in 1:n_planes
real_solns = solution.(Engine.Numerical.real_samples(system, freedom, rng = rng)) real_solns = solution.(Engine.Numerical.real_samples(system, freedom, rng = rng))