Formule

Réponses à toutes vos questions après le Bac (Fac, Prépa, etc.)
MacManus
Membre Irrationnel
Messages: 1365
Enregistré le: 28 Avr 2008, 14:41

formule

par MacManus » 06 Jan 2011, 14:23

Bonjour

Eh bien voila, je cherche à calculer le nombre de multiplications effectuées dans l'algorithme de factorisation de Cholesky (pour le calcul de la matrice triangulaire inférieure L).

En regardant de plus près dans le livre de P. Ciarlet (Introduction à l'analyse numérique matricielle et à l'optimisation, Masson) à la page 89 :

la formule affirme qu'il y a au total (n-1)+2(n-2)+...+(n-2)2+(n-1) = [n^(3) - n]/6 multiplications. Soit !

en y mettant la main à la patte, je trouve en calculant de mon côté :

( on sait que 1+2+...+(n-1) = n(n-1)/2 (ok) )

voici ce que je trouve :

n(n-1)/2 + (n-2)1 + (n-3)2 + (n-4)3 + ...+ ? +... + 3(n-4) + 2(n-3) + 1(n-2)

On s'aperçoit qu'il y a une symétrie entre le nb d'itérations et le nb de multiplications, ex : (n-2)1 et 1(n-2)

Ces formules coïncident-elles ?? (je ne parviens pas à le montrer...)

Merci bcp !!



XENSECP
Habitué(e)
Messages: 6387
Enregistré le: 27 Fév 2008, 19:13

par XENSECP » 06 Jan 2011, 14:39

Ta formule c'est :


Si tu développes :


Avec

Je te laisse terminer :)

MacManus
Membre Irrationnel
Messages: 1365
Enregistré le: 28 Avr 2008, 14:41

par MacManus » 06 Jan 2011, 14:57

Ok je suis tout à fait d'accord. En fait tu pars de de la première formule.
Avec ma formule, on aurait un truc du même genre ou presque :

je devrais alors retrouver la même chose , je vais vérifier.

merci bcp
j'ai manquer d'intuition sur ce coup là !!

MacManus
Membre Irrationnel
Messages: 1365
Enregistré le: 28 Avr 2008, 14:41

par MacManus » 06 Jan 2011, 15:02

Ok j'ai vérifié on retombe bien sur ta formule : un grand merci :)

 

Retourner vers ✯✎ Supérieur

Qui est en ligne

Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 113 invités

Tu pars déja ?



Fais toi aider gratuitement sur Maths-forum !

Créé un compte en 1 minute et pose ta question dans le forum ;-)
Inscription gratuite

Identification

Pas encore inscrit ?

Ou identifiez-vous :

Inscription gratuite