β Scribed by Beckermann, Labahn.
No coin nor oath required. For personal study only.
In this paper we provide a taet, numerically stable algorithm to determine when two given polynomials a arid b are relatively prime and remain relatively prime even after small perturbations of their coefficients. Such a problem is important in ninny applications where input data are only available up to a certain precision.Our methodβan extension of the Cabay-Meleshko algorithm for Pade approximationβis typically Π»ΠΈ order of magnitude faster than previously known stable methods. As such it may lie used an an inexpensive t.-sl which may lie applied before attempt ing to compute a ''numerical GCIJ'', in general a much more difficult task. We prove that the algorithm is numerically stable and give experiments verifying the numerical behaviour. Finally, we discuss possible extensions of our approach that can be applied to tin- problem of actually computing a numerical GCD.
π SIMILAR VOLUMES
Three new algorithms for multivariate polynomial GCD (greatest common divisor) are given. The first is to calculate a GrSbner basis with a certain term ordering. The second is to calculate the subresultant by treating the coefficients w.r.t, the main variable as truncated power series. The third is