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} ) <<<<<<< HEAD o5 = /var/folders/j_/gx42s4z576z_vj47_ym0j5xm0000gn/T/M2-7067-0/0/input ======= o5 = /var/folders/46/9b86vqxj4hjcngvy7kd7sb140000gn/T/M2-22983-0/0/input >>>>>>> origin/master o5 : File |
i6 : bertiniImageMonodromyCollect(storeBM2Files, NameStartFile=>"startSave", StorageFolder=>"StoreFiles", AffVariableGroup=>{x,y}, NumberOfLoops=>10,NumSolBound=>3, MonodromyStartParameters=>{1,0} ) Starting monodromy loop number: 1. Current fiber size: 1 Starting monodromy loop number: 2. Current fiber size: 2 Starting monodromy loop number: 3. Current fiber size: 2 Starting monodromy loop number: 4. Current fiber size: 2 Starting monodromy loop number: 5. Current fiber size: 3 Number of loops: 5. NumSolBound has been reached. o6 = {{1, 0}, {1, -.5+.866025*ii}, {1, -.5-.866025*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} ) <<<<<<< HEAD o9 = /var/folders/j_/gx42s4z576z_vj47_ym0j5xm0000gn/T/M2-7067-0/0/input ======= o9 = /var/folders/46/9b86vqxj4hjcngvy7kd7sb140000gn/T/M2-22983-0/0/input >>>>>>> origin/master o9 : File |
i10 : bertiniImageMonodromyCollect(storeBM2Files, AffVariableGroup=>{x,y}, StorageFolder=>"StoreFiles", MonodromyStartPoints=>{{1,0}}, NumberOfLoops=>10,NumSolBound=>3) Starting monodromy loop number: 1. Current fiber size: 2 Starting monodromy loop number: 2. Current fiber size: 3 Number of loops: 2. NumSolBound has been reached. o10 = {{1, 0}, {1, -.5+.866025*ii}, {1, -.5-.866025*ii}} o10 : List |
i11 : f="x^3+x*y^2+y+y^3+x-2"; h="a1*x+a2*y-1"; |
i13 : makeB'InputFile(storeBM2Files, B'Configs=>{{PARAMETERHOMOTOPY,2},{MPTYPE,2}},AffVariableGroup=>{x,y},ParameterGroup=>{a1,a2}, B'Polynomials=>{f,h} ) <<<<<<< HEAD o13 = /var/folders/j_/gx42s4z576z_vj47_ym0j5xm0000gn/T/M2-7067-0/0/input ======= o13 = /var/folders/46/9b86vqxj4hjcngvy7kd7sb140000gn/T/M2-22983-0/0/input >>>>>>> origin/master o13 : File |
i14 : s1=bertiniImageMonodromyCollect(storeBM2Files, AffVariableGroup=>{x,y}, MonodromyStartPoints=>{{1,0}}, MonodromyStartParameters=>{1,0}, NumberOfLoops=>10,NumSolBound=>3, MonodromyStartParameters=>{1,0}, ImageCoordinates=>({x})) Starting monodromy loop number: 1. Current fiber size: 1 Starting monodromy loop number: 2. Current fiber size: 1 Starting monodromy loop number: 3. Current fiber size: 1 Starting monodromy loop number: 4. Current fiber size: 1 Starting monodromy loop number: 5. Current fiber size: 1 Starting monodromy loop number: 6. Current fiber size: 1 Starting monodromy loop number: 7. Current fiber size: 1 Starting monodromy loop number: 8. Current fiber size: 1 Starting monodromy loop number: 9. Current fiber size: 1 Starting monodromy loop number: 10. Current fiber size: 1 NumberOfLoops has been reached. o14 = {{1, 0}} o14 : List |
i15 : #s1==1 o15 = true |
i16 : s2=bertiniImageMonodromyCollect(storeBM2Files, MonodromyStartPoints=>{{1,0}}, MonodromyStartParameters=>{1,0}, NumberOfLoops=>10,NumSolBound=>3, MonodromyStartParameters=>{1,0}, AffVariableGroup=>{x,y}, ImageCoordinates=>{"x^3+x*y^2+y+y^3+x-2"}) Starting monodromy loop number: 1. Current fiber size: 1 Starting monodromy loop number: 2. Current fiber size: 1 Starting monodromy loop number: 3. Current fiber size: 1 Starting monodromy loop number: 4. Current fiber size: 1 Starting monodromy loop number: 5. Current fiber size: 1 Starting monodromy loop number: 6. Current fiber size: 1 Starting monodromy loop number: 7. Current fiber size: 1 Starting monodromy loop number: 8. Current fiber size: 1 Starting monodromy loop number: 9. Current fiber size: 1 Starting monodromy loop number: 10. Current fiber size: 1 NumberOfLoops has been reached. o16 = {{1, 0}} o16 : List |
i17 : #s2==1 o17 = true |