Bonjour,

j'ai un projet sur la résolution d'un système linéaire avec la méthode de Cholesky.

j'ai fait ceci mais ça ne marche pas. Une aide serait la bienvenue.

Merci par avance.

Voir code:

void Cholesky(float** A, int n,float* b){
	float **r;float s;
	//r = malloc(n * sizeof(float));
	r = malloc2(n,n);
	for(int i = 0; i < n; i++){
		for(int j = 0; j < i-1;j++){
			s = A[i][i] - pow(r[j][i],2);
			printf("val de s: %f",s);
		}
			if(s <= 0){
				printf("La matrice n'est pas definie positive\n");
				exit(EXIT_SUCCESS);
			}else{
				r[i][i] = sqrt(s);
			}
			for(int j = i+1;j < n; j++){
				for(int k = 0; k <= i-1; k++){
					r[i][j] = (A[i][j] - r[k][i]*r[k][j])/r[i][i];
				} 
			}
		
	}
}

"Ca marche pas" ne veut rien dire.

Un message d'erreur, un exemple ?

Salut,

Qu'est-ce qui ne marche pas ? On peut voir malloc2() ? Qu'est censé faire ton code ? N'hésites pas à donner plus précisions pour ceux qui ne connaissent pas ton exo et qui n'ont jamais entendu parler de cholesky.

Merci pour les remarques

la méthode de Cholesky est une méthode de résolution de système linéaires type Ax = b elle fait une décomposition de la matrice de base en deux sous matrice R et Rt où R est triangulaire inférieure et Rt sa transposée

en fait la méthode fait rien même pas un message voici une capture

voici la méthode malloc2

 7 float * * malloc2 ( int rows , int cols )
  8 {
  9     float * * ptr = ( float * * ) malloc ( sizeof ( float * ) * rows + sizeo    f ( float ) * cols * rows ) ;
 10     float *   dat = ( float *   ) ( ptr + rows ) ;
 11     int i ;
 12     if ( ptr == NULL ) exit( EXIT_FAILURE ) ;
 13     for ( i = 0 ; i < rows ; ++ i , dat += cols ) ptr [ i ] = dat ;
 14     return ptr ;
 15 }

-
Edité par Taserface 23 octobre 2017 à 0:57:01

-
Edité par robun 23 octobre 2017 à 12:38:17

@robun merci pour le conseil sur la fonction pow on me l'avait déjà dit mais la flemme de changer en plus l'optimisation viendra après l'important c'est que le programme marche pour l'instant

j'me suis encore penché sur l'algo et je suis arrivé à ceci elle m'affiche bien une matrice supérieure droite mais toutes les valeurs ne sont pas bonnes.

void Cholesky(float** A,float* b,int n){
    (void) b;//we specify to the compilator we don't need b
    float **r;float s;
    r = malloc2(n,n);
   
    for(int i = 0; i < n; i++){
        
        float sommeR = 0;
      
        for(int j = 0; j < i-1;j++){
          
            sommeR += pow(r[j][i],2);
          
        }
        s = A[i][i] - sommeR;
        //if s <= 0 the matrice is not positive definite
        if(s <= 0){
          
            printf("the matrice is not positive definite\n");
          
            exit(EXIT_SUCCESS);
          
        }else{ 
          
            r[i][i] = sqrt(s);
          
        } 
        
       for(int j = i+1; j < n; j++){
         
         float multiplicationRki = 1,sommeRki = 0;
         
          for(int k = 0; k < i-1; k++){
          
           multiplicationRki = r[k][i] * r[k][j];
          
           sommeRki += multiplicationRki; 
           
         } 
         
         r[i][j] = (A[i][j] - sommeRki ) / r[i][i];
       }
   }
  /* if we want wa can print matrice element for see the transformation
  //ITS START HERE
   for(int i = 0; i < n; i++){
     for(int j = 0; j < n; j++){
       printf("%f\t",r[i][j]); 
     }
     printf("\n");
   }
  */ // HERE IS THE END OF PRINT
}

pour la ligne d'appel de la méthode elle est simple et n n'est pas < 0 puisque c'est l'utilisateur(pour l'instant moi ) qui la renseigne

Cholesky(A,b,n);

A étant la matrice(3*3) je précise b le vecteur et n la taille (qui est égal à 3)

NB: je travaille déjà sur des matrices pré-programmées de dimension 3 (Lotkin, Hilbert carré,...)

Ah, je le savais qu'il suffisait de regarder l'appel de la fonction : tu as mis les arguments dans le mauvais ordre !

(Oui mais le compilateur l'aurait signalé. En fait tu as juste mal recopié l'appel de la fonction et ça n'a rien à voir ?)

-
Edité par robun 24 octobre 2017 à 22:07:12

void Cholesky(float** A,float* b,int n)

Cholesky(A,b,n);

Je crois qu'on ne se comprend pas.

Hier, tu as écrit que l'appel était :

Cholesky(A,b,n);

alors que la définition de cette fonction commence par :

void Cholesky(float** A, int n,float* b){

(je ne fais que copier ton programme). Donc je me suis dit : voilà ! c'est pour ça que ça marche pas : il n'a pas écrit les arguments dans bon ordre ! Et puis, comme l'erreur est trop simple pour être vraie, je me suis dit ensuite : ou alors, quand il a écrit l'appel, il l'a mal recopié, mais dans le code source il est bien dans le bon ordre.

Et tu réponds (je te cite) :

« les options dans le mauvais ordre

void Cholesky(float** A,float* b,int n)

»

Apparemment, tu as donc rectifié ? C'est bien ça ? Et donc maintenant ça marche ? Voilà quelle était ma question, car à vrai dire ta réponse n'est pas claire (relis-toi, on dirait que tu admets que les options sont dans le mauvais ordre, et là tu dis que c'est dans le bon ordre, donc tu as rectifié ou pas ? et est-ce que c'était bien ça l'erreur ?)

-
Edité par robun 25 octobre 2017 à 21:52:22

OK. Donc le problème n'est toujours pas résolu ? 

Dans ce cas, il va falloir regarder plus en profondeur. Puisque tu as modifié le code, on devrait peut-être tout reprendre à zéro ? Une bonne exposition du problème serait :

1) Tu donnes un exemple de système linéaire simple dont tu indiques la solution (tu l'as calculée à la main).

2) Tu indiques ce que fait le programme (pour qu'on sache ce qui ne marche pas) et à quel moment précis il déconne (il suffit de comparer ses résultats intermédiaires avec ce que tu as calculé à la main).

3) Et tu re-postes le programme, ou en tout cas les parties concernées par les erreurs (la fonction, les appels à la fonction, éventuellement la définition des variables).

D'ailleurs normalement, si tu fais tout ça, tu trouvera les erreurs.

-
Edité par robun 26 octobre 2017 à 13:57:54

Enfin l'heure a sonné voici l'ensemble du code que j'ai fait les méthodes marchent.

Merci pour vos réponses

Lien vers mon github(Gauss+Jacobi): https://gist.github.com/Cirediallo/368db7a39f09c64ed108d4c39637a298

-
Edité par Taserface 2 novembre 2017 à 22:40:09