Discussion :

[ds6700]

bonjour;

j'ai Cherché sur wikipedia la décomposition QR selon la méthode de Householder et je l'ai implémenté avec mon type de matrice.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
typedef float E;
 
typedef struct matrix { 
	E * mat; 	
	int nb_rows,nb_columns;	
	}StrMatrix, *Matrix;

puis j'ai utilisé cette décomposition pour trouver les valeurs propres d’une matrice.

mon programme compile, et s'exécute mais mon problème et que la matrice que je saisi ne s'affiche pas et pourtant ma fonction "printMatrix(M)" me semble bien écrite.
et par la suite, toutes les fonctions qui viennent après ma fonction affiche dans mon "main.c" ne s'affichent pas aussi.

j'ai besoin de votre aide svp, car je doit avancer et il me reste à afficher les valeur propres de ma matrice puis le grand travail est de faire un affichage des ittérations avec gnuplot.

dans l'attente d'une réponse je vous remercie.


je ne comprends pas d'ou vient l'erreur

[Neckara]

Juste une question, si tu veux stocker des float dans une matrice à deux dimensions, pourquoi ne pas utiliser un E **Ainsi tu pourrais accéder à chacun de tes éléments comme ceci maMatrice[i][j];

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
 
    int n,i,j, p, q;
	E v;
    Matrix M=newMatrix(n,n);
    //Matrix V=newMatrix(n,1);
    Matrix Q=newMatrix(n,n);
    Matrix R=newMatrix(n,n);

n n'est pas initialisé il a donc une valeur indéterminée.

[ds6700]

la structure nous a été imposée.
je n'avais pas le choix

ce n'est pas le "n" qui pose pb, car n représente la taille de la matrice qui sera saisi par l'utilisateur (par moi).

[Neckara]

Si tu as des questions à poser, poses-les plutôt ici que par mp.

En effet, par mp tu poses la question à une seule personne alors qu'ici tu poses la question à tous tes lecteurs.

Pour retourner plusieurs matrices, tu as trois choix :
- Tu retournes une structure contenant plusieurs matrices
- Tu retournes un tableau de pointeur de matrices que tu aura alloué au préalable. La dernière valeur valant NULL pour marquer la fin du tableau.
- Tu passes en paramètre des pointeurs de pointeurs de matrices afin de pouvoir écrire sur l'adresse de la matrice sur le pointeur de matrice.

[diogene]

L'indexation est fausse

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
    for(i=1;i<=n;i++){
        for(j=1;j<=n;j++){
            // printf("            Ligne","i : 1, : " );
            printf("inserez mat(%d,%d):",i,j);
            scanf("%f", &v);
            setElt(M, i, j, v);
 
        }

setElt() fonctionne avec des indices basés sur 0, pas sur 1 -> on écrit en dehors de la zone allouée.

[ds6700]

Merci merci diogene, Oui c'est exactement ça le pb.
J'ai passé un temps fou à trouver la solution.
Je ne l'aurai jamais trouvé toute seule.
Encore merci.