compute the PuiseuxSeries solution of a regular chain
RegularChainBranches(rc, R, L)
RegularChainBranches(rc, R, L, V, coeff_opt, acc_opt)
regular chain of R
list of polynomials of R
(optional) variable to be used in the result
(optional) equation of the form coefficient = complex or coefficient = real (default: complex)
(optional) equation of the form accuracies = list_of_integers
The command RegularChainBranches(rc, R, L) returns the Puiseux series solutions of the regular chain rc around each of the roots of the polynomials in L.
The returned Puiseux series are actually truncated at a precision which is high enough to separate the different branches associated with the same root of a polynomial in L.
It is assumed that the coefficient field of R is the field of rational numbers.
It is assumed that rc is a one-dimensional strongly normalized regular chain.
This implies that every initial of a polynomial f in rc is either constant or univariate in a variable, say v, of R which is not algebraic w.r.t. rc.
It is assumed that the polynomials in L are univariate in that variable v. Moreover, it is assumed that every root of a polynomial in L is a root of the initial of a polynomial in rc.
Finally, it is assumed that V is an unassigned variable name.
The output of the command RegularChainBranches(rc, R, L) is a list of tuples. Each tuple is a Puiseux series solution of rc (associated with a root of the initial of a polynomial in rc) using an auxiliary variable to represent each Puiseux series as a pair of two polynomials; the first polynomial in such a pair is of a power of a variable (which is V if it is specified, otherwise T) and the second polynomial is univariate in that same variable.
In each Puiseux series solution of rc, each variable of R is given as a Puiseux series and all those Puiseux series share the same first polynomial in their representation; therefore, this polynomial appears only once as the first item in the tuple encoding this Puiseux series solution of rc.
If the option coefficient=real is present, then only the Puiseux series solutions with real coefficients are returned. Note those real Puiseux solutions are independent of the precision which is used in the returned truncated Puiseux series. In other words, the underlying algorithm guarantees that the returned Puiseux solutions represent all the real branches.
The accuracies option is for advanced users. If the option accuracies = L is present, where L is a list of integers, then L needs to have the same number of entries as the number of polynomials in rc. The effect is that the number of terms used for the ith polynomial in rc is L[i]. Any value used must comply with a certain set of requirements set out in the paper by Alvandi, Chen, and Moreno Maza listed in the References section.
If the accuracies option is not present, then a default value is computed. This value complies with the requirements listed in the paper.
This command is part of the RegularChains[AlgebraicGeometryTools] package, so it can be used in the form RegularChainBranches(..) only after executing the command with(RegularChains[AlgebraicGeometryTools]). However, it can always be accessed through the long form of the command by using RegularChains[AlgebraicGeometryTools][RegularChainBranches](..).
R ≔ PolynomialRing⁡x,y,z
rc ≔ Chain⁡−z2+y,x⁢z−y2,Empty⁡R,R
br ≔ RegularChainBranches⁡rc,R,z
rc ≔ Chain⁡y2⁢z+y+1,z+2⁢z⁢x2+y+1⁢x+1,Empty⁡R,R
rc ≔ Chain⁡y3−2⁢y3+y2+z5,z4⁢x+y3−y2,Empty⁡R,R
br ≔ RegularChainBranches⁡rc,R,z,coefficient=real
Parisa Alvandi, Changbo Chen, Marc Moreno Maza "Computing the Limit Points of the Quasi-component of a Regular Chain in Dimension One." Computer Algebra in Scientific Computing (CASC), Lecture Notes in Computer Science - 8136, (2013): 30-45.
Parisa Alvandi, Masoud Ataei, Mahsa Kazemi, Marc Moreno Maza "On the Extended Hensel Construction and its application to the computation of real limit points." J. Symb. Comput. 98: 120-162 (2020)
The RegularChains[AlgebraicGeometryTools][RegularChainBranches] command was introduced in Maple 2020.
For more information on Maple 2020 changes, see Updates in Maple 2020.
Download Help Document