English

Variable projection methods for approximate (greatest) common divisor computations

Optimization and Control 2015-11-05 v2 Numerical Analysis Numerical Analysis

Abstract

We consider the problem of finding for a given NN-tuple of polynomials (real or complex) the closest NN-tuple that has a common divisor of degree at least dd. Extended weighted Euclidean seminorm of the coefficients is used as a measure of closeness. Two equivalent representations of the problem are considered: (i) direct parameterization over the common divisors and quotients (image representation), and (ii) Sylvester low-rank approximation (kernel representation). We use the duality between least-squares and least-norm problems to show that (i) and (ii) are closely related to mosaic Hankel low-rank approximation. This allows us to apply to the approximate common divisor problem recent results on complexity and accuracy of computations for mosaic Hankel low-rank approximation. We develop optimization methods based on the variable projection principle both for image and kernel representation. These methods have linear complexity in the degrees of the polynomials for small and large dd. We provide a software implementation of the developed methods, which is based on a software package for structured low-rank approximation.

Keywords

Cite

@article{arxiv.1304.6962,
  title  = {Variable projection methods for approximate (greatest) common divisor computations},
  author = {Konstantin Usevich and Ivan Markovsky},
  journal= {arXiv preprint arXiv:1304.6962},
  year   = {2015}
}

Comments

32 pages, 4 figures

R2 v1 2026-06-22T00:06:28.263Z