Difference between revisions of "ApCoCoA-1:Bertini.BCMSolve"
m (replaced <quotes> tag by real quotes) |
|||
(9 intermediate revisions by 3 users not shown) | |||
Line 1: | Line 1: | ||
+ | {{Version|1}} | ||
<command> | <command> | ||
<title>Bertini.BCMSolve</title> | <title>Bertini.BCMSolve</title> | ||
Line 7: | Line 8: | ||
<description> | <description> | ||
<em>Please note:</em> The function(s) explained on this page is/are using the <em>ApCoCoAServer</em>. You will have to start the ApCoCoAServer in order to use it/them. | <em>Please note:</em> The function(s) explained on this page is/are using the <em>ApCoCoAServer</em>. You will have to start the ApCoCoAServer in order to use it/them. | ||
+ | <par/> | ||
+ | This function solves a polynomial system of equations using multihomogeneous homotopy. The polynomial system of equations must be quadratic. If the system has <tt>N</tt> variables then multihomogeneous homotopy will introduce <tt>N</tt> homogeneous variables to solve the system. It uses total degree homotopy to find all isolated solutions and enables you to use all kinds of user configurations provided by Bertini. Please read about configuration settings provided in Bertini manual. The system of polynomials should be non-homogeneous. The output will be the list of all finite solutions. | ||
− | |||
− | |||
<itemize> | <itemize> | ||
<item>@param <em>P</em>: List of polynomials of the given system.</item> | <item>@param <em>P</em>: List of polynomials of the given system.</item> | ||
− | <item>@param <em>ConfigSet</em>: List of strings representing configurations to be used. Bertini uses multiple configuration settings. These configurations should be provided by the user. If you want to use default configurations then leave this list empty. The use of configuration settings is little tricky but very | + | <item>@param <em>ConfigSet</em>: List of strings representing configurations to be used. Bertini uses multiple configuration settings. These configurations should be provided by the user. If you want to use default configurations then leave this list empty. The use of configuration settings is little tricky but very useful. You can switch from one kind of solving to the other just by changing or adding some configurations. For details about configuration settings see Bertini mannual <tt>http://www.nd.edu/~sommese/bertini/BertiniUsersManual.pdf</tt>.</item> |
− | <item>@return A list of lists containing the finite solutions of the polynomial system P.</item> | + | <item>@return A list of lists containing the finite solutions of the polynomial system <tt>P</tt>.</item> |
</itemize> | </itemize> | ||
Line 22: | Line 23: | ||
Use S ::= QQ[x,y]; | Use S ::= QQ[x,y]; | ||
P := [x^2+y^2-5,xy-2]; | P := [x^2+y^2-5,xy-2]; | ||
− | ConfigSet := [ | + | ConfigSet := ["MPTYPE: 2"]; |
-- Then we compute the solution with | -- Then we compute the solution with | ||
Line 53: | Line 54: | ||
<example> | <example> | ||
− | --The same | + | --The same example as above but this time using regeneration. |
--Regeneration is an equation-by-equation method for finding the non-singular isolated solutions. | --Regeneration is an equation-by-equation method for finding the non-singular isolated solutions. | ||
-- We want to solve the system x^2+y^2-5=0,xy-2=0, using multi-homogenization, for adaptive precision. | -- We want to solve the system x^2+y^2-5=0,xy-2=0, using multi-homogenization, for adaptive precision. | ||
Line 89: | Line 90: | ||
<example> | <example> | ||
− | -- We want to solve the non-homogenous polynomial system (29/16)z[1]^3 - 2z[1]z[2], z[2] - z[1]^2, using multi-homogenization. | + | -- We want to solve the non-homogenous polynomial system (29/16)z[1]^3 - 2z[1]z[2]=0, z[2] - z[1]^2=0, using multi-homogenization. |
− | -- In addition suppose we want Bertini to follow the paths very closely (correct to 8 digits before the endgame and correct to 11 -- digits when running the endgame) and that you want the endgame to compute approximations of each solution until successive approximations agree to at least 14 digits. | + | -- In addition suppose we want Bertini to follow the paths very closely (correct to 8 digits before the endgame and correct to 11 |
+ | -- digits when running the endgame) and that you want the endgame to compute approximations of each solution until successive | ||
+ | -- approximations agree to at least 14 digits. | ||
Use S ::= QQ[z[1..2]]; | Use S ::= QQ[z[1..2]]; | ||
Line 97: | Line 100: | ||
-- Then we compute the solution with | -- Then we compute the solution with | ||
− | Bertini. | + | Bertini.BCMSolve(P,ConfigSet); |
-- And we achieve a list of lists containing finite solutions. | -- And we achieve a list of lists containing finite solutions. | ||
---------------------------------------- | ---------------------------------------- | ||
+ | [ | ||
+ | [ | ||
+ | Vector(-6272142330887219272991421734532181259967888828557/2500000000000000000000000000000000000000000000000000000000000000, | ||
+ | -16866585769960937068249054027382922938807741513177/100000000000000000000000000000000000000000000000000000000000000000), | ||
+ | Vector(8805522460698904717783577526969974058842568093851/20000000000000000000000000000000000000000000000000000000000000000000000000, | ||
+ | 371666764789855282623508352529304991388858172316/6250000000000000000000000000000000000000000000000000000000000000000000000000000) | ||
+ | ], | ||
+ | [ | ||
+ | Vector(2800994767392818954296004874687536350261142130669/2000000000000000000000000000000000000000000000000000000000000000, | ||
+ | -10442004545151043511524382737808469250701933034019/5000000000000000000000000000000000000000000000000000000000000000), | ||
+ | Vector(-16863840390941893071858292951853672304120979133893/100000000000000000000000000000000000000000000000000000000000000000000000, | ||
+ | -5137795605139651161672232162161044043786444504329/1250000000000000000000000000000000000000000000000000000000000000000000000000000) | ||
+ | ], | ||
+ | [ | ||
+ | Vector(1745556907758497619927983789627389730399846491963/1250000000000000000000000000000000000000000000000000000000000000, | ||
+ | 28437259304210599902269601024330245142947311467321/10000000000000000000000000000000000000000000000000000000000000000), | ||
+ | Vector(-43119800109410418470261587328237339500316468515037/100000000000000000000000000000000000000000000000000000000000000000000000, | ||
+ | 11161237954525426581510768500892187905929449123739/2000000000000000000000000000000000000000000000000000000000000000000000000000000) | ||
+ | ] | ||
+ | ] | ||
Line 116: | Line 139: | ||
<seealso> | <seealso> | ||
− | <see>Introduction to CoCoAServer</see> | + | <see>ApCoCoA-1:Introduction to CoCoAServer|Introduction to CoCoAServer</see> |
− | <see>Bertini.BSolve</see> | + | <see>ApCoCoA-1:Bertini.BSolve|Bertini.BSolve</see> |
− | <see>Bertini.BMSolve</see> | + | <see>ApCoCoA-1:Bertini.BMSolve|Bertini.BMSolve</see> |
− | <see>Bertini.BUHSolve</see> | + | <see>ApCoCoA-1:Bertini.BUHSolve|Bertini.BUHSolve</see> |
</seealso> | </seealso> | ||
Line 125: | Line 148: | ||
<key>bertini.bcmsolve</key> | <key>bertini.bcmsolve</key> | ||
<key>solve zero dimensional polynomial system</key> | <key>solve zero dimensional polynomial system</key> | ||
− | <wiki-category>Package_bertini</wiki-category> | + | <wiki-category>ApCoCoA-1:Package_bertini</wiki-category> |
</command> | </command> |
Latest revision as of 13:27, 29 October 2020
This article is about a function from ApCoCoA-1. |
Bertini.BCMSolve
Solves a zero dimensional non-homogeneous polynomial system of equations using multi-homogenization and user configurations.
Syntax
Bertini.BCMSolve(P:LIST, ConfigSet:LIST)
Description
Please note: The function(s) explained on this page is/are using the ApCoCoAServer. You will have to start the ApCoCoAServer in order to use it/them.
This function solves a polynomial system of equations using multihomogeneous homotopy. The polynomial system of equations must be quadratic. If the system has N variables then multihomogeneous homotopy will introduce N homogeneous variables to solve the system. It uses total degree homotopy to find all isolated solutions and enables you to use all kinds of user configurations provided by Bertini. Please read about configuration settings provided in Bertini manual. The system of polynomials should be non-homogeneous. The output will be the list of all finite solutions.
@param P: List of polynomials of the given system.
@param ConfigSet: List of strings representing configurations to be used. Bertini uses multiple configuration settings. These configurations should be provided by the user. If you want to use default configurations then leave this list empty. The use of configuration settings is little tricky but very useful. You can switch from one kind of solving to the other just by changing or adding some configurations. For details about configuration settings see Bertini mannual http://www.nd.edu/~sommese/bertini/BertiniUsersManual.pdf.
@return A list of lists containing the finite solutions of the polynomial system P.
Example
-- We want to solve the system x^2+y^2-5=0,xy-2=0, using multi-homogenization, for adaptive precision. Use S ::= QQ[x,y]; P := [x^2+y^2-5,xy-2]; ConfigSet := ["MPTYPE: 2"]; -- Then we compute the solution with Bertini.BCMSolve(P,ConfigSet); -- And we achieve a list of lists containing finite solutions. ---------------------------------------- [ [ Vector(1000000000000017/1000000000000000, 145132717332349/15625000000000000000000000000), Vector(49999999999999/25000000000000, -3537662353156057/100000000000000000000000000000) ], [ Vector(-62500000000003/62500000000000, 4415730565392687/100000000000000000000000000000), Vector(-499999999999983/250000000000000, -66866973306543/400000000000000000000000000) ], [ Vector(999999999999983/500000000000000, -1787591178181031/50000000000000000000000000000), Vector(1000000000000013/1000000000000000, 281412486737749/25000000000000000000000000000) ], [ Vector(-499999999999999/250000000000000, -3956938527452181/1000000000000000000000000000000), Vector(-9999999999999989/10000000000000000, -596634837824491/1250000000000000000000000000000) ] ] --For Bertini output files please refer to ApCoCoA directory/Bertini.
Example
--The same example as above but this time using regeneration. --Regeneration is an equation-by-equation method for finding the non-singular isolated solutions. -- We want to solve the system x^2+y^2-5=0,xy-2=0, using multi-homogenization, for adaptive precision. Use S ::= QQ[x,y]; P := [x^2+y^2-5,xy-2]; ConfigSet := ["USEREGENERATION: 1"]; -- Then we compute the solution with Bertini.BCMSolve(P,ConfigSet); -- And we achieve a list of lists containing finite solutions. ---------------------------------------- [ [ Vector(9999999999999999/10000000000000000, -643977180168769/1250000000000000000000000000000), Vector(2, 1660674691787513/5000000000000000000000000000000) ], [ Vector(-2000000000000001/1000000000000000, 584020313856301/500000000000000000000000000000), Vector(-9999999999999999/10000000000000000, 45486167963413/125000000000000000000000000000) ], [ Vector(2, 2989952880295369/1000000000000000000000000000000), Vector(9999999999999993/10000000000000000, 732258034227497/5000000000000000000000000000000) ], [ Vector(-1, -879366755419571/5000000000000000000000000000000), Vector(-2, 4460430333228999/10000000000000000000000000000000) ] ] --For Bertini output files please refer to ApCoCoA directory/Bertini.
Example
-- We want to solve the non-homogenous polynomial system (29/16)z[1]^3 - 2z[1]z[2]=0, z[2] - z[1]^2=0, using multi-homogenization. -- In addition suppose we want Bertini to follow the paths very closely (correct to 8 digits before the endgame and correct to 11 -- digits when running the endgame) and that you want the endgame to compute approximations of each solution until successive -- approximations agree to at least 14 digits. Use S ::= QQ[z[1..2]]; P := [(29/16)z[1]^3 - 2z[1]z[2], z[2] - z[1]^2]; ConfigSet := ["TRACKTOLBEFOREEG: 1e-8", "TRACKTOLDURINGEG: 1e-11", "FINALTOL: 1e-14"]; -- Then we compute the solution with Bertini.BCMSolve(P,ConfigSet); -- And we achieve a list of lists containing finite solutions. ---------------------------------------- [ [ Vector(-6272142330887219272991421734532181259967888828557/2500000000000000000000000000000000000000000000000000000000000000, -16866585769960937068249054027382922938807741513177/100000000000000000000000000000000000000000000000000000000000000000), Vector(8805522460698904717783577526969974058842568093851/20000000000000000000000000000000000000000000000000000000000000000000000000, 371666764789855282623508352529304991388858172316/6250000000000000000000000000000000000000000000000000000000000000000000000000000) ], [ Vector(2800994767392818954296004874687536350261142130669/2000000000000000000000000000000000000000000000000000000000000000, -10442004545151043511524382737808469250701933034019/5000000000000000000000000000000000000000000000000000000000000000), Vector(-16863840390941893071858292951853672304120979133893/100000000000000000000000000000000000000000000000000000000000000000000000, -5137795605139651161672232162161044043786444504329/1250000000000000000000000000000000000000000000000000000000000000000000000000000) ], [ Vector(1745556907758497619927983789627389730399846491963/1250000000000000000000000000000000000000000000000000000000000000, 28437259304210599902269601024330245142947311467321/10000000000000000000000000000000000000000000000000000000000000000), Vector(-43119800109410418470261587328237339500316468515037/100000000000000000000000000000000000000000000000000000000000000000000000, 11161237954525426581510768500892187905929449123739/2000000000000000000000000000000000000000000000000000000000000000000000000000000) ] ] --For Bertini output files refer to ApCoCoA directory/Bertini. ------------------------------------------
See also