compute the triangular decomposition of a polynomial system
Triangularize(F, H, R)
Triangularize(F, rc, R)
Triangularize(F, H, rc, R)
Triangularize(F, R, 'normalized'='yes')
Triangularize(F, R, 'radical'='yes')
Triangularize(F, R, 'output'='lazard')
Triangularize(F, R, 'probability'='prob')
list or set of equations
(optional) list or set of inequations
(optional) regular chain
(optional) boolean flag
(optional) numerical value
The command Triangularize(F,H,rc,R) returns a triangular decomposition of the set of the zeros of F that are zeros of rc but not zeros of H.
If H is not specified, it is set to 1.
If rc is not specified, it is set to the empty regular chain.
If H= and rc is empty, this command returns a triangular decomposition of the set of the zeros of F.
Each element of this triangular decomposition is a regular chain.
The concepts of a regular chain and a triangular decomposition are defined in the page RegularChains.
The algorithm for the command Triangularize is described in the paper "On Triangular Decompositions of Algebraic Varieties" by Marc Moreno Maza, available on the author's web page.
This command is part of the RegularChains package, so it can be used in the form Triangularize(..) only after executing the command with(RegularChains). However, it can always be accessed through the long form of the command by using RegularChains[Triangularize](..).
If 'normalized'='yes', each of the returned regular chains is normalized. If 'normalized'='strongly', each regular chain is strongly normalized. By default, each regular chain may not be normalized.
If 'radical'='yes', the saturated ideal of each regular chain is radical. By default, the saturated ideal of each regular chain may not be radical.
This option is currently only supported if the characteristic of the polynomial ring R is zero.
If 'output'='lazard', the decomposition is the sense of Lazard; otherwise, it is made in the sense of Kalkbrener. The latter is weaker but less expensive to compute. When F has only a finite number of solutions, both senses coincide.
As mentioned above, if this option is specified and if inequations (the input polynomial list H) are present, then the output is a constructible set. See the commands of the subpackage ConstructibleSetTools for operations on constructible sets.
If a non-empty regular chain rc is provided, then the decomposition is the sense of Lazard, whether 'output'='lazard' is specified or not.
If 'probability'='prob', the probabilistic algorithm of Dahan, Moreno Maza, Schost, Wu, and Xie (ISSAC 2005) is used. This algorithm applies only to square systems in characteristic zero. In this implementation, it will fail if two primes such that the system modulo those primes is radical cannot be found after ten tries. This option is not compatible with any other options.
The options 'output'='lazard' and 'normalized'='yes' cannot be mixed in the current implementation of this command. (This limitation will be solved in the next version of the RegularChains library).
The commands BivariateModularTriangularize and ChangeOfOrder implement algorithms for computing triangular decompositions under particular circumstances. When they apply, they are likely to outperform Triangularize.
Define a ring of polynomials.
R ≔ PolynomialRing⁡z,y,x
Define a set of polynomials of R. Each of them will be viewed as an equality to 0.
sys ≔ x2+y+z−1,x+y2+z−1,x+y+z2−1
Ideally, you would like to decompose the set of the common solutions of sys into a list of points. The Triangularize command does this by using symbolic expressions. Sometimes several points are grouped together in a generic one, as in this example. These groups of points are called regular chains, and they are grouped together because they share some mathematical properties.
dec ≔ Triangularize⁡sys,R
Since regular chains may contain large expressions, their output form is just a word. To view their members, use the Equations command.
The last three regular chains are very simple: each of them clearly corresponds to a single point. The first regular chain corresponds to two points, because its univariate polynomial in the "smallest" variable x has two roots.
Consider now another polynomial ring and another polynomial system.
R ≔ PolynomialRing⁡x,y,a,b,c,d,g,h
sys ≔ a⁢x+b⁢y−g,c⁢x+d⁢y−h
In the polynomial ring, the ordering on the variables is such that x>y>a>b>c>d>g>h. Solving sys with this ordering implies that you want to express x and y as functions of the other variables. Hence you can view the system sys as a parametric linear system with two equations and two unknowns, x and y. Applying RegularChains[Triangularize] displays the generic solution, which is similar to the solution given by Groebner[Solve].
dec ≔ Triangularize⁡sys,R;map⁡Equations,dec,R
This generic solution assumes that the determinant of the system is not zero. With the option output=lazard, Triangularize gives all of the solutions, including those that cancel the determinant of the input system.
decl ≔ Triangularize⁡sys,R,output=lazard
You already know that each regular chain is associated with a set of equations. It is also associated with a set of inequations.
rc ≔ decl1;Equations⁡rc,R;Inequations⁡rc,R
The inequations of a regular chain rc are the set of the initials of the polynomials of rc. In the first regular chain above, the inequations are d⁢a−b⁢c and c. Hence, for this regular chain, none of these two polynomials should vanish. The solutions of sys that cancel either c or the determinant d⁢a−b⁢c are given by the other regular chains of decl. Below, for each regular chain of decl, we print its list of equations together with its set of inequations.
Assume now that you want to see g and h as transcendental quantities; that is, quantities that cannot satisfy any polynomial equations. Then you need to redefine the polynomial ring as follows.
R2 ≔ PolynomialRing⁡x,y,a,b,c,d,g,h
dec2 ≔ Triangularize⁡sys,R2,output=lazard
Now, you can obtain five regular chains, none of them imposing a condition on g or h. The following example uses the option probability, by which a triangular decomposition is done by using a modular algorithm.
R ≔ PolynomialRing⁡x,y,z
sys ≔ −x5+y5−3⁢y−1,5⁢y4−3,−20⁢x+y−z
Moreno Maza, M. "On Triangular Decompositions of Algebraic Varieties." MEGA-2000 conference. Bath, UK, England.
Download Help Document