apply the Weierstrass preparation theorem to a univariate polynomial over power series
univariate polynomial over power series generated by this package
The command WeierstrassPreperation(u) applies the Weierstrass preparation theorem to the univariate polynomial over power series u. It returns two univariate polynomials over power series p and a, described below.
Let d be the smallest integer such that the d-th coefficient of u is a unit. If there is no such coefficient, in other words, if each coefficient of u has a constant coefficient equal to 0, then Maple signals an error.
The polynomial p is a Weierstrass polynomial of degree d, thus p is monic and every coefficient of p, other than the leading one, evaluates to zero at the origin, that is, it belongs to the maximal ideal of the ring of power series. The polynomial a regarded as a power series (even with respect to its main variable) is invertible. Furthermore, u = p * a.
The command WeierstrassPreperation(u) is based on lazy evaluation techniques. To be precise, this command returns as soon as all coefficients (which are power series) of the univariate polynomials p and a are initialized, that is, as soon as the generators of all those coefficients have been determined.
This command is supported for univariate polynomials over power series, but not for univariate polynomials over Puiseux series.
When using the MultivariatePowerSeries package, do not assign anything to the variables occurring in the power series, Puiseux series, and univariate polynomials over these series. If you do, you may see invalid results.
We create a cubic univariate polynomial over power series and apply the Weierstrass preparation theorem to it. The quadratic coefficient of f is the first unit, so p1 will be quadratic and a1 will be linear.
f ≔ UnivariatePolynomialOverPowerSeries⁡PowerSeries⁡x,PowerSeries⁡y,PowerSeries⁡1,1PowerSeries⁡1+x+y,z
p1,a1 ≔ WeierstrassPreparation⁡f
p1,a1≔UnⅈvarⅈatⅇPolynomⅈalOvⅇrPowⅇrSⅇrⅈⅇs: 0+…+0+…⁢z+1⁢z2,UnⅈvarⅈatⅇPolynomⅈalOvⅇrPowⅇrSⅇrⅈⅇs: 1+…+1+…⁢z
The terms of p1 and a1 can only be computed in tandem, so if we update the precision of p1, then the precision of a1 will also be updated.
We compute the product of p1 and a1, and verify that its coefficients are equal at precision 10.
h1 ≔ p1⁢a1
This cubic univariate polynomial over power series has non-unit constant, linear, and quadratic coefficients of its main variable. Only the cubic coefficient is a unit. Hence if we apply the Weierstrass preparation theorem, p2 will be cubic and a2 will be independent of z.
g ≔ UnivariatePolynomialOverPowerSeries⁡PowerSeries⁡x,PowerSeries⁡y,PowerSeries⁡x,GeometricSeries⁡x,y,z
p2,a2 ≔ WeierstrassPreparation⁡g
p2,a2≔UnⅈvarⅈatⅇPolynomⅈalOvⅇrPowⅇrSⅇrⅈⅇs: 0+…+0+…⁢z+0+…⁢z2+1⁢z3,UnⅈvarⅈatⅇPolynomⅈalOvⅇrPowⅇrSⅇrⅈⅇs: 1+x+y+…
We multiply the factors together and verify that the coefficients are equal to gs at precision 10.
h2 ≔ Multiply⁡p2,a2
Alexander Brandt, Mahsa Kazemi, Marc Moreno Maza. "Power Series Arithmetic with the BPAS Library." Computer Algebra in Scientific Computing (CASC), Lecture Notes in Computer Science. Vol. 12291 (2020): 108-128.
Mohammadali Asadi, Alexander Brandt, Mahsa Kazemi, Marc Moreno Maza, and Erik Postma. "Multivariate Power Series in Maple." Corless R.M., Gerhard J., Kotsireas I.S. (eds) Maple in Mathematics Education and Research. MC 2020. Communications in Computer and Information Science (CCIS). Vol. 1414 Springer (2021): 48-66.
The MultivariatePowerSeries[WeierstrassPreparation] command was introduced in Maple 2021.
For more information on Maple 2021 changes, see Updates in Maple 2021.
The MultivariatePowerSeries[WeierstrassPreparation] command was updated in Maple 2023.
Download Help Document