Bonjour

Je cherche à calculer la racine carrée d'une matrice symétrique positive (les valeurs propres sont positives ou nulles) d'une certaine matrice A. Il y a deux définition équivalentes.

La racine carrée de A provient de la diagonalisation de A : on remplace ses valeurs propres par leur racine carrée. En utilisant cette définition, les temps de calculs sont énormes...

A=B*B, avec B matrice symétrique. Cette écriture de A est unique.

Matlab possède une routine qui s'appelle SQRTM, mais malheureusement les temps de calculs sont astronomiques dès que l'on manipule des matrices trop grosses (je bosse avec des matrices 4000 x 4000). J'ai lu le code mais je ne suis pas vraiment sûr de la méthode utilisée...

Existe-t-il d'autres algorithmes moins coûteux?

Merci beaucoup !