Comptes Rendus
no. 1-2
Numerical Analysis
Computing the matrix sign and absolute value functions
[Calcul du signe et de la valeur absolue d'une matrice]

Présenté par : Olivier Pironneau

Michaël Ndjinga1, 2
1 Commissariat à l'énergie atomique, centre de Saclay, CEA/DEN/DM2S/SFME, 91191 Gif-sur-Yvette, France
2 Laboratoire de mathématiques appliquées aux systèmes, École centrale Paris, 92295 Châtenay-Malabry cedex, France
Comptes Rendus. Mathématique, Volume 346 (2008) no. 1-2, pp. 119-124.

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.

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.

Reçu le :
Accepté le :
Publié le :
DOI : 10.1016/j.crma.2007.11.028
Michaël Ndjinga 1, 2

1 Commissariat à l'énergie atomique, centre de Saclay, CEA/DEN/DM2S/SFME, 91191 Gif-sur-Yvette, France
2 Laboratoire de mathématiques appliquées aux systèmes, École centrale Paris, 92295 Châtenay-Malabry cedex, France 
@article{CRMATH_2008__346_1-2_119_0,
     author = {Micha\"el Ndjinga},
     title = {Computing the matrix sign and absolute value functions},
     journal = {Comptes Rendus. Math\'ematique},
     pages = {119--124},
     publisher = {Elsevier},
     volume = {346},
     number = {1-2},
     year = {2008},
     doi = {10.1016/j.crma.2007.11.028},
     language = {en},
}
TY  - JOUR
AU  - Michaël Ndjinga
TI  - Computing the matrix sign and absolute value functions
JO  - Comptes Rendus. Mathématique
PY  - 2008
SP  - 119
EP  - 124
VL  - 346
IS  - 1-2
PB  - Elsevier
DO  - 10.1016/j.crma.2007.11.028
LA  - en
ID  - CRMATH_2008__346_1-2_119_0
ER  - 
%0 Journal Article
%A Michaël Ndjinga
%T Computing the matrix sign and absolute value functions
%J Comptes Rendus. Mathématique
%D 2008
%P 119-124
%V 346
%N 1-2
%I Elsevier
%R 10.1016/j.crma.2007.11.028
%G en
%F CRMATH_2008__346_1-2_119_0
Michaël Ndjinga. Computing the matrix sign and absolute value functions. Comptes Rendus. Mathématique, Volume 346 (2008) no. 1-2, pp. 119-124. doi : 10.1016/j.crma.2007.11.028. https://comptes-rendus.academie-sciences.fr/mathematique/articles/10.1016/j.crma.2007.11.028/

[1] J.M. Delhaye; G. Hewitt; M. Zuber Numerical Benchmark Tests in Multiphase Science and Technology, vol. 3, Hemisphere, Washington, DC/New York, 1987

[2] D.A. Drew; S.L. Passman Theory of Multicomponents Fluids, Springer-Verlag, New York, 1999

[3] T. Gallouët; J.-M. Hérard; N. Seguin Numerical modeling of two-phase flows using the two-fluid two-pressure approach, Math. Models Methods Appl. Sci., Volume 14 (2004) no. 5

[4] J.M. Ghidaglia; A. Kumbaro; G. Le Coq On the numerical solution to two fluid models via a cell centered finite volume method, Eur. J. Mech. B/Fluids, Volume 20 (2001) no. 6

[5] G.H. Golub; C.F. Van Loan Matrix Computation, The John Hopkins University Press, 1996

[6] C.S. Kenney; A.J. Laub The matrix sign function, IEEE Trans. Automat. Control, Volume 40 (1995) no. 8, pp. 1330-1348

[7] R.J. LeVeque Finite Volume Methods for Hyperbolic Problems, Cambridge University Press, 2002

[8] J.M. Masella; I. Faille; T. Gallouët On an approximate Godunov scheme, Int. J. Comput. Fluid Dynamics, Volume 12 (1999), pp. 133-149

[9] M. Ndjinga Influence of interfacial pressure on the hyperbolicity of the two-fluid model, C. R. Acad. Sci. Paris, Ser. I, Volume 344 (2007), pp. 407-412

[10] M. Ndjinga, A. Kumbaro, P. Laurent-Gengoux, F. De Vuyst, 2006, Fast and robust computation of the matrix absolute value function. Application to Roe solver for the numerical simulation of two-phase flow models, in: 14th International Conference on Nuclear Energy, 17–20 July 2006, Miami, USA

[11] P.L. Roe Approximate Riemannn solvers, parameter vectors and difference schemes, J. Comput. Phys., Volume 43 (1981)

[12] J. Stoer; R. Bulirsch Introduction to Numerical Analysis, Springer-Verlag, 2002

[13] I. Toumi; A. Kumbaro An approximate linearized Riemann solver for a two-fluid model, J. Comput. Phys., Volume 2 (1996), p. 124

[14] J.H. Wilkinson; C. Reinsch Handbook for Automatic Computation, vol. II, Linear Algebra, Springer-Verlag, New York, 1971

Cité par Sources :

Commentaires - Politique

Ces articles pourraient vous intéresser

Hardware (DNA) circuits

Richard D'Ari; René Thomas

C. R. Biol (2003)

Propelling soft objects

Hakim Boukellal; Julie Plastino; Vincent Noireaux; ...

C. R. Phys (2003)

Revisiting poromechanics in the context of microporous materials

Gilles Pijaudier-Cabot; Romain Vermorel; Christelle Miqueu; ...

C. R. Méca (2011)