SNAP
EuclideanReduction
compute the smallest degree pair of univariate polynomials by Euclidean-like unimodular reduction
Calling Sequence
Parameters
Description
Examples
References
EuclideanReduction(a, b, z, tau = eps, out)
a, b
-
univariate numeric polynomials
z
name; indeterminate for a and b
tau = eps
(optional) equation where eps is of type numeric and non-negative; stability parameter
out
(optional) equation of the form output = obj where obj is 'UR' or a list containing one or more instances of this name; select result objects to compute
The EuclideanReduction(a, b, z) command returns the last numerically well-conditioned basis accepted by the Coprime algorithm [2]. This corresponds to the smallest degree pair of polynomials in the sequence of numerically well-behaved polynomial remainders that can be obtained from (a,b) by unimodular reduction.
It thus provides the user with a pair of polynomials that generates the same ideal generated by (a,b) but with degrees that are, in general, much smaller. Furthermore, the highest degree component of such a reduced pair is a good candidate for an epsilon-GCD of (a,b).
The optional stability parameter tau can be set to any non-negative value eps to control the quality of the output. Decreasing eps yields a more reliable solution. Increasing eps reduces the degrees of the returned basis.
As specified by the out option, Maple returns an expression sequence containing the following:
* UR contains a 2 by 2 unimodular matrix polynomial U in z such that a,b.U=a'⁢,b'⁢ where (a', b') is the last basis accepted by the algorithm of [2].
with⁡SNAP:
a≔z6−12.4⁢z5+50.18112+62.53⁢z4−163.542⁢z3+232.9776⁢z2−170.69184⁢z
b≔z5−17.6⁢z4+118.26⁢z3−372.992⁢z2−274.09272+538.3333⁢z
EuclideanReduction⁡a,b,z
4.⁢z4−45.3201452919812⁢z3+182.643498183852⁢z2−301.305647387541⁢z+164.902292707462,3.48999306545107⁢z3−25.4412393835560⁢z2+55.5055044834605⁢z−34.9173360478195
EuclideanReduction⁡a,b,z,τ=1.×10−8
0.250000000000000⁢z2−0.875000000000161⁢z+0.659999999999995,−1.07691633388640×10−14⁢z−5.54556400800266×10−14
Beckermann, B., and Labahn, G. "A fast and numerically stable Euclidean-like algorithm for detecting relatively prime numerical polynomials." Journal of Symbolic Computation. Vol. 26, (1998): 691-714.
Beckermann, B., and Labahn, G. "When are two numerical polynomials relatively prime?" Journal of Symbolic Computation. Vol. 26, (1998): 677-689.
See Also
SNAP[DistanceToCommonDivisors]
SNAP[EpsilonGCD]
Download Help Document