Kaltofen's division-free determinant algorithm differentiated for matrix adjoint computation

Gilles Villard 1, 2
2 ARENAIRE - Computer arithmetic
Inria Grenoble - Rhône-Alpes, LIP - Laboratoire de l'Informatique du Parallélisme
Abstract : Kaltofen has proposed a new approach in 1992 for computing matrix determinants without divisions. The algorithm is based on a baby steps/giant steps construction of Krylov subspaces, and computes the determinant as the constant term of a characteristic polynomial. For matrices over an abstract ring, by the results of Baur and Strassen, the determinant algorithm, actually a straight-line program, leads to an algorithm with the same complexity for computing the adjoint of a matrix. However, the latter adjoint algorithm is obtained by the reverse mode of automatic differentiation, hence somehow is not ''explicit''. We present an alternative (still closely related) algorithm for the adjoint thatcan be implemented directly, we mean without resorting to an automatic transformation. The algorithm is deduced by applying program differentiation techniques ''by hand'' to Kaltofen's method, and is completely decribed. As subproblem, we study the differentiation of programs that compute minimum polynomials of lineraly generated sequences, and we use a lazy polynomial evaluation mechanism for reducing the cost of Strassen's avoidance of divisions in our case.
Type de document :
Article dans une revue
Journal of Symbolic Computation, Elsevier, 2011, 46 (7), pp.773-790. 〈10.1016/j.jsc.2010.08.012〉
Liste complète des métadonnées

https://hal-ens-lyon.archives-ouvertes.fr/ensl-00335918
Contributeur : Gilles Villard <>
Soumis le : vendredi 31 octobre 2008 - 10:39:48
Dernière modification le : vendredi 20 avril 2018 - 15:44:24
Document(s) archivé(s) le : lundi 7 juin 2010 - 20:52:16

Fichiers

detdiff.pdf
Fichiers produits par l'(les) auteur(s)

Identifiants

Collections

Citation

Gilles Villard. Kaltofen's division-free determinant algorithm differentiated for matrix adjoint computation. Journal of Symbolic Computation, Elsevier, 2011, 46 (7), pp.773-790. 〈10.1016/j.jsc.2010.08.012〉. 〈ensl-00335918〉

Partager

Métriques

Consultations de la notice

269

Téléchargements de fichiers

230