HighPerformanceMathematics/homotopy-continuation
2
0
Fork 0
mirror of https://github.com/BachoSeven/homotopy-continuation.git synced 2026-10-06 22:55:15 +00:00
Code Issues Packages Projects Releases Wiki Activity

chore: update plots after step reset fix

Browse Source
This commit is contained in:
BachoSeven
2023-08-24 19:06:42 +02:00
parent 6edc93d215
commit 5b47562d12
4 changed files with 4 additions and 14 deletions
Download Patch File Download Diff File Expand all files Collapse all files
plots/solutions1.png
BIN
View File
Binary file not shown.
Side by Side Swipe Overlay

Before

Width:  |  Height:  |  Size: 16 KiB

After

Width:  |  Height:  |  Size: 16 KiB

plots/solutions2.png
BIN
View File
Binary file not shown.
Side by Side Swipe Overlay

Before

Width:  |  Height:  |  Size: 14 KiB

After

Width:  |  Height:  |  Size: 14 KiB

plots/solutions3.png
BIN
View File
Binary file not shown.
Side by Side Swipe Overlay

Before

Width:  |  Height:  |  Size: 16 KiB

After

Width:  |  Height:  |  Size: 16 KiB

solve.jl
+4 -14
View File
@@ -20,7 +20,6 @@ function solve(F, (G, roots) = start_system(F), maxsteps=10000)
# F=homogenize(F) # F=homogenize(F)
H=homotopy(F,G) H=homotopy(F,G)
solutions = [] solutions = []
step_array = []
@time Threads.@threads for r in roots @time Threads.@threads for r in roots
t = 1.0 t = 1.0
@@ -37,10 +36,9 @@ function solve(F, (G, roots) = start_system(F), maxsteps=10000)
steps += 1 steps += 1
end end
push!(solutions, x0) push!(solutions, x0)
push!(step_array, steps)
end end
return (solutions, step_array) return solutions
end end
# Input polynomial system # Input polynomial system
@@ -49,17 +47,9 @@ F = [x*y - 1, x^2 + y^2 - 4]
T = [x*y - 1, x^2 + y^2 - 2] T = [x*y - 1, x^2 + y^2 - 2]
C = [x^3 - y + 5x^2 - 10, 2x^2 - y - 10] C = [x^3 - y + 5x^2 - 10, 2x^2 - y - 10]
(sF, sf) = solve(F) sF = filter(u -> imag(u[1]) < 0.1 && imag(u[2]) < 0.1, solve(F))
(sT, st) = solve(T) sT = filter(u -> imag(u[1]) < 0.1 && imag(u[2]) < 0.1, solve(T))
(sC, sc) = solve(C) sC = filter(u -> imag(u[1]) < 0.1 && imag(u[2]) < 0.1, solve(C))
println(sf)
println(st)
println(sc)
sF = filter(u -> imag(u[1]) < 0.1 && imag(u[2]) < 0.1, sF)
sT = filter(u -> imag(u[1]) < 0.1 && imag(u[2]) < 0.1, sT)
sC = filter(u -> imag(u[1]) < 0.1 && imag(u[2]) < 0.1, sC)
# Plotting the system and the real solutions # Plotting the system and the real solutions
ENV["GKSwstype"]="nul" ENV["GKSwstype"]="nul"