An extended GCRD algorithm for parametric univariate polynomial matrices and application to parametric Smith form

The first extended greatest common right divisor (GCRD) algorithm for parametric univariate polynomial matrices is presented. The starting point of this GCRD algorithm is the free property of submodules over univariate polynomial rings. We convert the computation of GCRDs to that of free basis for modules and prove that a free basis of the submodule generated by row vectors of input matrices forms just a GCRD of these matrices. The GCRD algorithm is obtained by computing a minimal Gröbner basis for the corresponding submodule since a minimal Gröbner basis of submodules is a free basis for univariate cases. While the key idea of extended algorithm is to construct a special module by adding the unit vectors which can record the representation coefficients. This method based on modules can be naturally generalized to the parametric case because of the comprehensive Gröbner systems for modules. As a consequence, we obtain an extended GCRD algorithm for parametric univariate polynomial matrices. More importantly, we apply the proposed extended GCD algorithm for univariate polynomials (as a special case of matrices) to the computation of Smith normal form, and give the first algorithm for reducing a univariate polynomial matrix with parameters to its Smith normal form.