Computing the matrix sign and absolute value functions

We present two algorithms for the computation of the matrix sign and absolute value functions. Both algorithms avoid a complete diagonalisation of the matrix, but they however require some informations regarding the eigenvalues location. The first algorithm consists in a sequence of polynomial iterations based on appropriate estimates of the eigenvalues, and converging to the matrix sign if all the eigenvalues are real. Convergence is obtained within a finite number of steps when the eigenvalues are exactly known. Nevertheless, we present a second approach for the computation of the matrix sign and absolute value functions, when the eigenvalues are exactly known. This approach is based on the resolution of an interpolation problem, can handle the case of complex eigenvalues and appears to be faster than the iterative approach. To cite this article: M. Ndjinga, C. R. Acad. Sci. Paris, Ser. I 346 (2008).

RésuméNous présentons deux algorithmes pour le calcul des fonctions signe et valeur absolue matricielles. Ces algorithmes évitent une diagonalisation complète de la matrice, mais nécessitent des estimations de ses valeurs propres. Le premier algorithme consiste en une suite d'itérations polynomiales construite à partir d'approximations des valeurs propres, et convergeant vers le signe de la matrice lorsque les valeurs propres sont réelles. La convergence s'obtient en un nombre fini d'itérations lorsqu'elles sont connues exactement. Nous présentons cependant une seconde approche s'appuyant sur la résolution d'un problème d'interpolation polynomiale, pour le calcul des fonctions signe et valeur absolue matricielles dans le cas où les valeurs propres sont connues. Ce second algorithme se révèle plus rapide que le premier, et permet la prise en compte de matrices ayant des valeurs propres complexes. Pour citer cet article : M. Ndjinga, C. R. Acad. Sci. Paris, Ser. I 346 (2008).