C++

[C++]Factorisation Cholesky qui marche pas

Recherche :

Mot : Pseudo : Filtrer
Bas de page
Auteur	Sujet : [C++]Factorisation Cholesky qui marche pas

danikool

Bonsoir,

J'ai essayé d'écrire une procédure qui me donne la factorisation Cholesky d'une matrice supposé symétrique définie positive ( factorise une matrice en 2 matrices triangulaire dont l'une est transposé de l'autre) mais il y a un problème: Le programme compile bien (0 erreur qu'il me dit) mais lorsque je lance le programme pour une matrice de dimension supérieur ou égal à 3, ça fait planter MS-DOS !

Je précise que j'utilise codeblocks sous windows.

Voici le bout de code correspondant à la procédure:

Code :

void Cholesky(int n, float A[N][N], float L[N][N], float Lt[N][N])
{
L[0][0]=sqrt(A[0][0]);
float somme;
for(int i=1;i<n;i++)
{
for(int j=0;j<i;j++)
{
somme=0;
for(int k=0;k<j;j++)
{
somme=somme+L[i][k]*L[j][k];
}
L[i][j]=(A[i][j]-somme)/L[j][j];
}
somme=0;
for(int k=0;k<i;k++)
{
somme=somme+L[i][k]*L[i][k];
}
L[i][i]=sqrt(A[i][i]-somme);
}
}

(j'ai défini la valeur de N à 10 au début du programme)

Voilà, j'espère que vous pourrez m'aider à trouver d'ou vient le problème.

Merci d'avance!

Publicité

Amonchakai

Salut !

tu est sûr, que c'est bien j++ qu'il faut mettre comme incrémentation dans ton for ligne 11 et pas plutôt k++ ? Car là le le vois pas sortir de cette boucle...

Message cité 1 fois

danikool

Amonchakai a écrit :

Salut !

tu est sûr, que c'est bien j++ qu'il faut mettre comme incrémentation dans ton for ligne 11 et pas plutôt k++ ? Car là le le vois pas sortir de cette boucle...

lol !
Merci, heureusement que tu me l'as fait remarqué. J'étais vraiment trop fatigué hier ..

Je viens de tester et ca plante plus en effet.

rui

Strike Out Looking..

le numerical recipes est ton ami sinon

FORUM HardWare.fr

Programmation

C++

[C++]Factorisation Cholesky qui marche pas

Sujets relatifs
Formulaire mailto ne marche pas sur IE	[PHP] Marche sous Easyphp mais pas sous wamp?
Script vbs qui marche pas	[SQL] noob, script arrêt/marche service SQL
[python] control+c ne marche pas ???	[C] fgets ne marche qu'une fois...
Avis aux flasheurs experts // comment ca marche PreviewCity	[Résolu]Client/serveur qui marche qu'en local
[C++] "typage dynamique" marche pas	boost::bind, pourquoi this ne marche pas?
Plus de sujets relatifs à : [C++]Factorisation Cholesky qui marche pas

Page générée en 0.060 secondes