Engine prototype #13
@ -21,7 +21,7 @@ end
|
||||
# === example ===
|
||||
|
||||
# the Lorentz form
|
||||
Q = diagm([-1, 1, 1, 1, 1])
|
||||
Q = diagm([1, 1, 1, 1, -1])
|
||||
|
||||
# initialize the partial gram matrix for an arrangement of seven spheres in
|
||||
# which spheres 1 through 5 are mutually tangent, and spheres 3 through 7 are
|
||||
@ -54,11 +54,11 @@ gram[1, 6] = gram[6, 1]
|
||||
# in this initial guess, the mutual tangency condition is satisfied for spheres
|
||||
# 1 through 5
|
||||
guess = sqrt(0.5) * BigFloat[
|
||||
1 1 1 1 2 0.2 0.1;
|
||||
0 0 0 0 -sqrt(6) 0.3 -0.2;
|
||||
1 1 -1 -1 0 -0.1 0.3;
|
||||
1 -1 1 -1 0 -0.5 0.4;
|
||||
1 -1 -1 1 0 0.1 -0.2
|
||||
1 -1 -1 1 0 0.1 -0.2;
|
||||
0 0 0 0 -sqrt(6) 0.3 -0.2;
|
||||
1 1 1 1 2 0.2 0.1;
|
||||
]
|
||||
|
||||
# search parameters
|
||||
|
Loading…
Reference in New Issue
Block a user