Give spheres positive radii in examples

This changes the meaning of `indep_val` in the overlapping pyramids
example, so we adjust `indep_val` to get a nice-looking construction.
This commit is contained in:
Aaron Fenyes 2024-07-17 17:22:33 -07:00
parent 2038103d80
commit 4728959ae0
3 changed files with 8 additions and 14 deletions

View File

@ -213,8 +213,8 @@ elements = begin
1 1 -1 -1 0 1 1 -1 -1 0
1 -1 1 -1 0 1 -1 1 -1 0
1 -1 -1 1 0 1 -1 -1 1 0
-0.5 -0.5 -0.5 -0.5 -a-1 0.5 0.5 0.5 0.5 a+1
0.5 0.5 0.5 0.5 -a+1 -0.5 -0.5 -0.5 -0.5 a-1
] ]
end end

View File

@ -23,13 +23,7 @@ end
gram = sparse(J, K, values) gram = sparse(J, K, values)
# set the independent variable # set the independent variable
# indep_val = 2//5
# using gram[6, 2] or gram[7, 1] as the independent variable seems to stall
# convergence, even if its value comes from a known solution, like
#
# gram[6, 2] = 0.9936131705272925
#
indep_val = -9//5
gram[6, 1] = BigFloat(indep_val) gram[6, 1] = BigFloat(indep_val)
gram[1, 6] = gram[6, 1] gram[1, 6] = gram[6, 1]
@ -43,8 +37,8 @@ guess = begin
1 1 -1 -1 0 1 1 -1 -1 0
1 -1 1 -1 0 1 -1 1 -1 0
1 -1 -1 1 0 1 -1 -1 1 0
-0.5 -0.5 -0.5 -0.5 -a-1 0.5 0.5 0.5 0.5 a+1
0.5 0.5 0.5 0.5 -a+1 -0.5 -0.5 -0.5 -0.5 a-1
] + 0.2*Engine.rand_on_shell(fill(BigFloat(-1), 5)), ] + 0.2*Engine.rand_on_shell(fill(BigFloat(-1), 5)),
Engine.rand_on_shell(fill(BigFloat(-1), 2)) Engine.rand_on_shell(fill(BigFloat(-1), 2))
) )

View File

@ -28,7 +28,7 @@ for j in 1:6
push!(values, -1/BigFloat(3)) push!(values, -1/BigFloat(3))
filled = true filled = true
else else
push!(values, -1) push!(values, 1)
filled = true filled = true
end end
if filled if filled
@ -46,8 +46,8 @@ guess = hcat(
1 1 -1 -1 0 1 1 -1 -1 0
1 -1 1 -1 0 1 -1 1 -1 0
1 -1 -1 1 0 1 -1 -1 1 0
0 0 0 0 -1.5 0 0 0 0 1.5
1 1 1 1 -0.5 1 1 1 1 0.5
] + 0.2*Engine.rand_on_shell(fill(BigFloat(-1), 5)), ] + 0.2*Engine.rand_on_shell(fill(BigFloat(-1), 5)),
BigFloat[0, 0, 0, 0, 1] BigFloat[0, 0, 0, 0, 1]
) )