RegularChains
RealTriangularize
compute all the solutions of a semi-algebraic system
LazyRealTriangularize
compute the solutions of a semi-algebraic system interactively
SamplePoints
compute sample points for a semi-algebraic system
Calling Sequence
Parameters
Options
Description
Examples
References
Compatibility
RealTriangularize(sys, R, options)
RealTriangularize(F, N, P, H, R, options)
LazyRealTriangularize(sys, R, options)
LazyRealTriangularize(F, N, P, H, R, options)
SamplePoints(sys, R, options)
SamplePoints(F, N, P, H, R, options)
SamplePoints(C, R)
sys
-
list of polynomial equations, inequations, and inequalities over R
F
list of polynomials over R representing equations
N
list of polynomials over R representing non-negativity conditions
P
list of polynomials over R representing positivity conditions
H
list of polynomials over R representing inequations
C
object of type cadcell
R
polynomial ring
options
(optional) equation of the form 'output'=format, where format is either 'piecewise', 'record', 'list' , 'zerodimensional', or equation of the form 'isolation'='algorithm', where 'algorithm'='Discoverer' or 'algorithm'='VincentCollinsAkritas' , or equation of the form 'method'='algorithm', where 'algorithm'='recursive' or 'algorithm'='incremental' , or equation of the form 'relaxation'=algorithm, where 'algorithm'='explicit' or 'algorithm'='implicit' .
The output option controls the form of the output.
If 'output'='list' is specified, the output is a list of objects of type regular_semi_algebraic_system, for RealTriangularize and LazyRealTriangularize, and a list of objects of type box, for SamplePoints. This is the default for RealTriangularize and for SamplePoints.
If 'output'='piecewise' is specified, then the output format is a case discussion tree encoded by a piecewise object. This is not a valid option for SamplePoints, and it is the default for LazyRealTriangularize.
If 'output'='record' is specified, then the output is a sequence of Maple records.
If 'output'='zerodimensional' is specified, then the output is a list of regular semi-algebraic sets. (See SemiAlgebraicSetTools for the definition of a regular semi-algebraic system.) The usage of this format is restricted to the case where the constructible set associated to S has finitely many solutions. This condition holds in particular when the semi-algebraic system S itself has finitely many solutions.
The isolation option controls the method used for isolating the real roots of a zero-dimensional regular chain. (See RealRootIsolate for the subject of real root isolation of a regular chain.)
If 'isolation'='VincentCollinsAkritas' is specified, then the method is a generalization of the Vincent-Collins-Akritas algorithm which isolates real roots for univariate polynomials. The techniques are very close to the ones used by Renaud Rioboo in his ISSAC 1992 paper (see reference below).
If 'isolation'='Discoverer' is specified, then the method implements the algorithm of the paper "An Algorithm for Isolating the Real Solutions of Semi-algebraic Systems." by Bican Xia and Lu Yang. The default method is Discoverer.
The method option controls the decomposition algorithm.
If 'method'='recursive' is specified, then the algorithm presented in the ISSAC 2010 paper by Changbo Chen et al. (see reference below) is used.
If 'method'='incremental' is specified, then the algorithm presented in the ISSAC 2011 paper by Changbo Chen et al. (see reference below) is used. This option is the default.
The relaxation option controls the use of a relaxation technique aiming at reducing the number of components in the output. This strategy is described in the ISSAC 2011 paper by Changbo Chen et al. (see reference below).
If 'relaxation'='explicit' is specified, then criteria are used to check whether relaxation techniques apply. These criteria have some computational cost. However, relaxation techniques might speed up the decomposition process dramatically in some cases.
If 'relaxation'='implicit' is specified, then relaxation techniques are applied in a light way, at essentially no cost. However, some opportunities for relaxation can only discovered in the 'explicit' mode. The option 'relaxation'='implicit' is used by default.
Consider the semi-algebraic system S whose unknowns are the variables of R and whose equations, non-negative inequalities, strictly positive inequalities, and inequations are given either by sys or by F, N, P, and H, respectively. The solutions of S are the real numbers which satisfy simultaneously the equations, inequalities, and inequations of S.
The commands RealTriangularize(F, N, P, H, R) or RealTriangularize(sys, R) return a triangular decomposition of the solutions of S in the following sense. They compute a list of regular semi-algebraic systems such that the union of their solution sets is exactly the solution set of S. (See SemiAlgebraicSetTools for the definition of a regular semi-algebraic system.)
The commands LazyRealTriangularize(F, N, P, H, R) or LazyRealTriangularize(sys, R) return regular semi-algebraic systems and unevaluated recursive calls in a piecewise format, unless a different output format is requested. The zero sets of those regular semi-algebraic systems form a "lazy" real triangular decomposition of S in the following sense.
The zero sets of those regular semi-algebraic systems form a subset U of the solution set of S;
The dimension of the solution set of S∖U is strictly less than the dimension of the constructible set associated with S.
By evaluating the recursive calls in finitely many steps, one can eventually obtain a triangular decomposition of S.
When the option 'output'='list' is specified, the recursive calls will be discarded and only the regular semi-algebraic systems are returned as a list.
The commands SamplePoints(F, N, P, H, R) and SamplePoints(sys, R) return at least one sample point per real connected component of S.
The output format is the same as for the command RealRootIsolate. That is, a list of so-called boxes, each box containing one and only one sample point. (See SemiAlgebraicSetTools or RealRootIsolate for the definition of a box.)
The BoxValues command retrieves a detailed description of a box.
In addition, each box can be refined with an arbitrary small precision. See the commands RefineBox and RefineListBox.
The command SamplePoints(C,R) returns a sample point of the CAD cell given by C. The type cadcell encodes any CAD cell and is used for one of the output formats of the command CylindricalAlgebraicDecompose. See SemiAlgebraicSetTools for the definition of a CAD cell.
with⁡RegularChains:
Define a ring of polynomials.
R≔PolynomialRing⁡x,y,z:
Define a set of equations.
F≔x3+y+z−1,x+y3+z−1,x+y+z3−1
F≔x3+y+z−1,y3+x+z−1,z3+x+y−1
Define a set of non-negative polynomials.
N≔
Define a set of positive polynomials.
P≔
Define a set of inequations.
H≔x
Compute the solutions of F for which x is nonzero.
dec≔RealTriangularize⁡F,N,P,H,R
dec≔regular_semi_algebraic_system,regular_semi_algebraic_system,regular_semi_algebraic_system,regular_semi_algebraic_system,regular_semi_algebraic_system
Display⁡dec,R
x−z=0y−z=0z3+2⁢z−1=0,x−1=0y+1=0z−1=0,x+1=0y−1=0z−1=0,x−1=0y=0z=0,x−1=0y−1=0z+1=0
Compute the solutions of F for which x≠1.
dec≔RealTriangularize⁡F,N,P,x−1,R
dec≔regular_semi_algebraic_system,regular_semi_algebraic_system,regular_semi_algebraic_system,regular_semi_algebraic_system
x−z=0y−z=0z3+2⁢z−1=0,x=0y=0z−1=0,x+1=0y−1=0z−1=0,x=0y−1=0z=0
Compute the solutions of F for which x is negative.
dec≔RealTriangularize⁡F,N,−x,,R
dec≔regular_semi_algebraic_system
x+1=0y−1=0z−1=0
For input systems which are known to have finitely many solutions like F, the option 'output'='zerodimensional' provides an alternative output format based on isolation boxes, as for RealRootIsolate.
dec≔RealTriangularize⁡F,R,output=zerodimensional
dec≔regular_semi_algebraic_set,regular_semi_algebraic_set,regular_semi_algebraic_set,regular_semi_algebraic_set,regular_semi_algebraic_set,regular_semi_algebraic_set,regular_semi_algebraic_set
map⁡Display,dec,R
x=38,12y=38,12z=2964,117256,x=1y=1z=−1,x=0y=1z=0,x=0y=0z=1,x=1y=−1z=1,x=−1y=1z=1,x=1y=0z=0
Consider the unit circle in the Euclidean plane.
R≔PolynomialRing⁡y,x;sys≔x2+y2=1
R≔polynomial_ring
sys≔x2+y2=1
Compute its real points.
RealTriangularize⁡sys,R,output=record
y2+x2−1=0x<1x+1>0,y=0x−1=0,y=0x+1=0
Consider the generic equation of degree two.
R≔PolynomialRing⁡x,c,b,a;sys≔a⁢x2+b⁢x+c=0
sys≔a⁢x2+b⁢x+c=0
Compute a triangular decomposition of the 4-variable hypersurface it defines.
dec≔RealTriangularize⁡sys,R
a⁢x2+b⁢x+c=0−4⁢c⁢a+b2>0anda≠0,2⁢a⁢x+b=04⁢a⁢c−b2=0a≠0,b⁢x+c=0a=0b≠0,c=0b=0a=0
Consider the piecewise output format.
RealTriangularize⁡sys,R,output=piecewise
x2⁢a+x⁢b+c=00<−4⁢a⁢c+b2∧a≠0x⁢b+c=0,a=0b≠0c=0,b=0,a=0b=0a=02⁢a⁢x+b=0,4⁢a⁢c−b2=0a≠0c=0,b=0,a=0a=0−4⁢a⁢c+b2=0otherwise
Consider the record output format.
a⁢x2+b⁢x+c=0−4⁢c⁢a+b2>0a≠0,2⁢a⁢x+b=04⁢a⁢c−b2=0a≠0,b⁢x+c=0b≠0a=0,c=0b=0a=0
Use LazyRealTriangularize to start the decomposition.
dec≔LazyRealTriangularize⁡sys,R
dec≔a⁢x2+b⁢x+c=00<−4⁢c⁢a+b2∧a≠0LazyRealTriangularize⁡a=0,a⁢x2+b⁢x+c=0,polynomial_ringa=0LazyRealTriangularize⁡−4⁢c⁢a+b2=0,a⁢x2+b⁢x+c=0,polynomial_ring−4⁢c⁢a+b2=0otherwise
Go one step further, computing components in lower dimension.
dec2≔value⁡dec
dec2≔x2⁢a+x⁢b+c=00<−4⁢a⁢c+b2∧a≠0x⁢b+c=0,a=0b≠0LazyRealTriangularize⁡a=0,b=0,x⁢b+c=0,x2⁢a+x⁢b+c=0,polynomial_ringb=0a=02⁢a⁢x+b=0,4⁢a⁢c−b2=0a≠0LazyRealTriangularize⁡a=0,−4⁢a⁢c+b2=0,4⁢a⁢c−b2=0,2⁢a⁢x+b=0,x2⁢a+x⁢b+c=0,polynomial_ringa=0−4⁢a⁢c+b2=0otherwise
Go the last step, computing components in dimension zero.
value⁡dec2
Use the list output format.
dec≔LazyRealTriangularize⁡sys,R,output=list
a⁢x2+b⁢x+c=0−4⁢c⁢a+b2>0anda≠0
Use the record output format.
dec≔LazyRealTriangularize⁡sys,R,output=record
dec≔a⁢x2+b⁢x+c=0−4⁢c⁢a+b2>0a≠0,LazyRealTriangularize⁡a=0,a⁢x2+b⁢x+c=0,polynomial_ring,output=record,LazyRealTriangularize⁡−4⁢c⁢a+b2=0,a⁢x2+b⁢x+c=0,polynomial_ring,output=record
Go one step further.
dec2≔a⁢x2+b⁢x+c=0−4⁢c⁢a+b2>0a≠0,b⁢x+c=0b≠0a=0,LazyRealTriangularize⁡a=0,b=0,b⁢x+c=0,a⁢x2+b⁢x+c=0,polynomial_ring,output=record,2⁢a⁢x+b=04⁢a⁢c−b2=0a≠0,LazyRealTriangularize⁡a=0,−4⁢a⁢c+b2=0,4⁢a⁢c−b2=0,2⁢a⁢x+b=0,a⁢x2+b⁢x+c=0,polynomial_ring,output=record
Go the last step.
a⁢x2+b⁢x+c=0−4⁢a⁢c+b2>0a≠0,b⁢x+c=0b≠0a=0,c=0b=0a=0,2⁢a⁢x+b=04⁢a⁢c−b2=0a≠0,c=0b=0a=0
Consider again the unit circle in the Euclidean plane.
Compute sample points of it.
P≔SamplePoints⁡sys,R
P≔box,box,box,box
Display⁡P,R
y=−1x=0,y=1x=0,y=0x=1,y=0x=−1
Request points with positive abscissa.
sys≔x2+y2=1,0<x
P≔box,box,box
y=−111128,−5564x=12,y=5564,111128x=12,y=0x=1
Return to the generic equation of degree two.
Compute sample points of the 4-variable hypersurface it defines.
P≔box,box,box,box,box,box,box,box,box
x=−1c=12b=0a=−12,x=1c=12b=0a=−12,x=−1c=−12b=0a=12,x=1c=−12b=0a=12,x=0c=0b=−12a=0,x=0c=0b=12a=0,x=0c=0b=0a=0,x=0c=0b=0a=−12,x=0c=0b=0a=12
The relaxation option.
with⁡RegularChains:with⁡SemiAlgebraicSetTools:
R≔PolynomialRing⁡x,a,b
sys≔x2+a+b⁢x+a=0,0<x−a
st≔time⁡:rtd0≔RealTriangularize⁡sys,R:time⁡−st
0.311
st≔time⁡:rtd1≔RealTriangularize⁡sys,R,relaxation=explicit:time⁡−st
0.652
evalb⁡nops⁡rtd1<nops⁡rtd0
true
Accessing a sample point from a CAD cell.
cad≔CylindricalAlgebraicDecompose⁡x2−a=0,R,output=cadcell
cad≔cad_cell,cad_cell,cad_cell
cc≔cad1;Display⁡cc,R
cc≔cad_cell
x=0a=0b=b
SamplePoints⁡cc,R,output=record
x=0a=0b=0
F. Boulier, C. Chen, F. Lemaire, M. Moreno Maza, "Real root isolation of regular chains." ASCM'2009, Math-for-Industry, Lecture Note Series Vol. 22.
C. Chen, J. H. Davenport, J. P. May, M. Moreno Maza, B. Xia, R. Xiao, "Triangular Decomposition of Semi-algebraic Systems." Proceedings of ISSAC 2010, ACM Press, 2010.
C. Chen, J. H. Davenport, M. Moreno Maza, B. Xia, R. Xiao, "Computing with semi-algebraic sets represented by triangular decomposition". Proceedings of the 2011 International Symposium on Symbolic and Algebraic Computation (ISSAC 2011), ACM Press, pp. 75--82, 2011.
R. Rioboo, Computation of the real closure of an ordered field. Proceedings of ISSAC'92, Academic Press, San Francisco.
L. Yang, X. Hou, B. Xia, A complete algorithm for automated discovering of a class of inequality-type theorems. Sci. China Ser. F, Vol. 44 (2001): 33-49.
B. Xia, L. Yang, An Algorithm for Isolating the Real Solutions of Semi-algebraic Systems. J. Symb. Comp., Vol. 34 No. 5 (2001): 461-477.
B. Xia, T. Zhang, Real Solution Isolation Using Interval Arithmetic. Computers and Mathematics with Applications, Vol. 52 (2006): 853-860.
The RegularChains[LazyRealTriangularize] and RegularChains[SamplePoints] commands were introduced in Maple 15.
For more information on Maple 15 changes, see Updates in Maple 15.
The C parameter was introduced in Maple 16.
The and options were introduced in Maple 16.
For more information on Maple 16 changes, see Updates in Maple 16.
See Also
BorderPolynomial
BoxValues
CylindricalAlgebraicDecompose
Display
RealRootClassification
RealRootIsolate
RefineBox
RefineListBox
SemiAlgebraicSetTools
Triangularize
Download Help Document