Given a directory that has a Bertini input file that has ParameterHomotopy set to 2, a start file, for startparameters, this function uses parameter homotopies to perform a monodromy homotopy.
i1 : f="x^3+x*y^2+y+y^3+x-2"; h="a1*x+a2*y-1"; |
i3 : writeStartFile(storeBM2Files,{{1,0}},NameStartFile=>"startSave");
|
i4 : writeParameterFile(storeBM2Files,{1,0},NameParameterFile=>"start_parameters");
|
i5 : makeB'InputFile(storeBM2Files,
B'Configs=>{{PARAMETERHOMOTOPY,2},{MPTYPE,2}},AffVariableGroup=>{x,y},ParameterGroup=>{a1,a2}, B'Polynomials=>{f,h} )
o5 = /var/folders/j_/gx42s4z576z_vj47_ym0j5xm0000gn/T/M2-74671-0/0/input
o5 : File
|
i6 : b'PHMonodromyCollect(storeBM2Files,
NameStartFile=>"startSave",
NameSolutionsFile=>"simple_raw_solutions",
NumberOfLoops=>10,NumSolBound=>3,
MonodromyStartParameters=>{1,0} )
/var/folders/j_/gx42s4z576z_vj47_ym0j5xm0000gn/T/M2-74671-0/0/
(Monodromy loop number, 1)
~/bertini
~/bertini
~/bertini
(Monodromy loop number, 2)
~/bertini
~/bertini
~/bertini
(Monodromy loop number, 3)
~/bertini
~/bertini
~/bertini
(Monodromy loop number, 4)
~/bertini
~/bertini
~/bertini
(Monodromy loop number, 5)
~/bertini
~/bertini
~/bertini
3
o6 = {{1, -.5-.866025*ii}, {1, -.5+.866025*ii}, {1,
------------------------------------------------------------------------
-8.6739e-17-9.48317e-16*ii}}
o6 : List
|
i7 : f="x^3+x*y^2+y+y^3+x-2"; h="a1*x+a2*y-1"; |
i9 : makeB'InputFile(storeBM2Files,
B'Configs=>{{PARAMETERHOMOTOPY,2},{MPTYPE,2}},AffVariableGroup=>{x,y},ParameterGroup=>{a1,a2}, B'Polynomials=>{f,h} )
o9 = /var/folders/j_/gx42s4z576z_vj47_ym0j5xm0000gn/T/M2-74671-0/0/input
o9 : File
|
i10 : b'PHMonodromyCollect(storeBM2Files,
StorageFolder=>"StoreFiles",
MonodromyStartPoints=>{{1,0}},
NumberOfLoops=>10,NumSolBound=>3,
MonodromyStartParameters=>{1,0} )
/var/folders/j_/gx42s4z576z_vj47_ym0j5xm0000gn/T/M2-74671-0/0/
(Monodromy loop number, 1)
~/bertini
~/bertini
~/bertini
(Monodromy loop number, 2)
~/bertini
~/bertini
~/bertini
(Monodromy loop number, 3)
~/bertini
~/bertini
~/bertini
3
o10 = {{1, -.5-.866025*ii}, {1, -.5+.866025*ii}, {1,
-----------------------------------------------------------------------
7.11849e-16-3.45223e-16*ii}}
o10 : List
|