Fast Computation of Shifted Popov Forms of Polynomial Matrices via Systems of Modular Polynomial Equations release_gtdz4yfwjvflhdc3vw4rno76fm

by Vincent Neiger

Released as a article .

2016  

Abstract

We give a Las Vegas algorithm which computes the shifted Popov form of an m × m nonsingular polynomial matrix of degree d in expected O(m^ω d) field operations, where ω is the exponent of matrix multiplication and O(·) indicates that logarithmic factors are omitted. This is the first algorithm in O(m^ω d) for shifted row reduction with arbitrary shifts. Using partial linearization, we reduce the problem to the case d <σ/m where σ is the generic determinant bound, with σ / m bounded from above by both the average row degree and the average column degree of the matrix. The cost above becomes O(m^ωσ/m ), improving upon the cost of the fastest previously known algorithm for row reduction, which is deterministic. Our algorithm first builds a system of modular equations whose solution set is the row space of the input matrix, and then finds the basis in shifted Popov form of this set. We give a deterministic algorithm for this second step supporting arbitrary moduli in O(m^ω-1σ) field operations, where m is the number of unknowns and σ is the sum of the degrees of the moduli. This extends previous results with the same cost bound in the specific cases of order basis computation and M-Padé approximation, in which the moduli are products of known linear factors.
In text/plain format

Archived Files and Locations

application/pdf  279.7 kB
file_hbu6hn2fd5atvgfna2qvg3owbu
arxiv.org (repository)
web.archive.org (webarchive)
Read Archived PDF
Preserved and Accessible
Type  article
Stage   accepted
Date   2016-05-12
Version   v2
Language   en ?
arXiv  1602.00710v2
Work Entity
access all versions, variants, and formats of this works (eg, pre-prints)
Catalog Record
Revision: 73ab2130-63e9-4668-979d-e7359b90ed5e
API URL: JSON