Discussion :

[wintercold]

Bonjour,
J'ai un program qui n'a aucun prob lors de la compilation mé à l'exécution, un warning est affiché:



Mon compilateur (C++ Builder) met une flêche jaune a coté de la ligne où réside le problème (que je met en grave au dessous). Pourriez-vous m'éclairer les fautes? s'il vous plait c'est très urgent, je verifie souvent la page pour vos réponses.
Merci beaucoup.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
void qrdcmp(float **A, int N, float *c, float *d, int *sing)
{
	int i,j,k;
	float scale,sigma,sum,tau;
	
	*sing=0;
	for (k=0;k<N-1;k++) 
		{
		scale=0.0;
		for (i=k;i<N;i++) 
			{
			if (scale>=fabs(A[i][k])) 
				{
				scale=scale;
				} 
			else 
				{
				scale=fabs(A[i][k]);
				}
			}
		
		if (scale == 0.0) 
			{ 
			*sing=1;
			c[k]=d[k]=0.0;
			} 
		else 
			{ 
			for (i=k;i<N;i++) A[i][k] /= scale;
			for (sum=0.0,i=k;i<=N;i++) sum += (A[i][k])*(A[i][k]);			if (A[k][k]>=0) 
			{
				sigma=sqrt(sum);
			}
			else 
			{
				sigma=-sqrt(sum);
			}
			A[k][k] += sigma;
			c[k]=sigma*A[k][k];
			d[k] = -scale*sigma;
			for (j=k+1;j<N;j++) 
				{
				for (sum=0.0,i=k;i<=N;i++) sum += A[i][k]*A[i][j];				tau=sum/c[k];
				for (i=k;i<N;i++) A[i][j] -= tau*A[i][k];
				}
			}
		}
	d[N]=A[N-1][N-1];
	if (d[N-1] == 0.0) *sing=1;
}

Vous trouverez ci-joint le fichier du program complet si vous voulez l'éxécuter.

[Obsidian]

Segfault :

Code :

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

1
2
3
4
void qrdcmp(float **A, int N, float *c, float *d, int *sing)
{
  ...
*sing=0;

Tu essaies d'écrire zéro à l'adresse indiquée par le pointeur « sing », lequel t'est passé en paramètre. Fort bien, mais cela n'empêche pas la fonction appelante de te transmettre n'importe quoi.

Le problème se situe vraissemblablement au niveau d'au dessus.

[wintercold]

merci de votre réponse.
En fait, à part de cette fonction, le reste du program se compose:
- une fonction d'allocation mem pour matrice
- une autre pour sa libration
- une fonction saisir matrice
- une fonction afficher matrice
- fonction main

g deja testé av succès ces 5 fonctions avant d'y integer void qrdcmp(float **A, int N, float *c, float *d, int *sing)

(/* Les prototypes de fonctions utilisé */)

Code :

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

1
2
3
4
5
6
int Menu ( void );
float** Allocate_Mat_Memory ( int );
void Free_Mat_Memory ( float**, int );
void Affiche_Mat ( float**, int );
void Saisi_Mat ( float**, int );
void qrdcmp(float **A, int N, float *c, float *d, int *sing);

le prob doit se situer ds la fonction qrdcmp et en effet le compilateur a positionné sur la ligne

Code :

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

	*sing=0;

Vous pouvez le vérifier en compilant le fichier essaie.c (ci-joint) par C++ Builder.

J'attends tjrs vos réponses.

[Obsidian]

g deja testé av succès ces 5 fonctions avant d'y integer void qrdcmp(float **A, int N, float *c, float *d, int *sing)

le prob doit se situer ds la fonction qrdcmp et en effet le compilateur a positionné sur la ligne

Comme je te l'ai dit, le problème vient de ce qui appelle « qrdcmp » car la valeur du cinquième paramètre (« sing ») est invalide. Le code est tout-à-fait correct sur la forme, et c'est pour cela qu'il compile sans erreur, mais la valeur du paramètre ne l'est pas.

Je jette un œil sur ton source.

MàJ : dans main(), on trouve :

Code C :

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

1
2
3
4
5
int *sing = 0; /* Ligne 149 */
 
...
 
qrdcmp(matrice1, N, c, d, sing); /* Ligne 180 */

Donc, tu définis un pointeur sur int dans ta fonction main() et tu l'initialises avec l'adresse 0x00000000. Ensuite, tu passes ce pointeur à ta fonction, qui essaie donc d'écrire à cette adresse. Donc, ça plante. Et si tu lis le message d'erreur, tu t'aperçois qu'il te l'explique en clair : « Write of address 00000000 ».

Peut-être voulais-tu que ta fonction stocke son résultat dans une variable. Dans ce cas, il aurait fallu déclarer un entier et non un pointeur sur un entier, soit int sing = 0; sans l'astérisque « * », et passer l'adresse de cette variable à ta fonction, avec « &sing » plutôt que « sing » seul.

[wintercold]

peut-être g trouvé la sortie pour le souci sur la ligne *sing=0.
En fait, ds la declaration des variables ds la fonction main, javais initializé *sing=0.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
void main()//void)
{ 
    FILE *fpt;
    float *c=0, *d=0;
    int *sing=0;
    float **matrice1;
    int N=1;
    char caractere;

    fpt = fopen("notesma11.dat","r");
	if(fpt == NULL)
	{
		printf("\nError: File couldn't be opened (1).\n");
		exit(-1);
	}
	do
        {
            caractere = fgetc(fpt); // On lit le caractère
          
			if(caractere == '\n')
			{
				N=N+1;
			}
			

        } while (caractere != EOF);
    fclose(fpt);

          
                printf("Nombre de lignes: %d  \n",N);     
                matrice1 = Allocate_Mat_Memory( N );
                         
                Saisi_Mat( matrice1, N );
                printf("La Matrice saisie:\n\n");
                Affiche_Mat( matrice1, N);
		qrdcmp(matrice1, N, c, d, sing);
Affiche_Mat( matrice1, N);
....

Quand je met seulement int *sing; ce souci est enlevé mais g tjrs le même warning et maintenant la faute se trouve à la ligne en grave (ds la fonction qrdcmp):

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
...
if (scale == 0.0) 
			{ 
			*sing=1;
			c[k]=d[k]=0.0;
			} 
		else 
			{ 
			for (i=k;i<N;i++) A[i][k] /= scale;
			[B]for (sum=0.0,i=k;i<=N;i++) sum += (A[i][k])*(A[i][k]);[/B]
			if (A[k][k]>=0) 
			{
				sigma=sqrt(sum);
			}
			else 
			{
				sigma=-sqrt(sum);
			}
...

[gl]

Quand je met seulement int *sing; ce souci est enlevé

Non, tu n'as rien résolu ici. Le code invoque un comportement indéfini, il peut sembler fonctionner correctement, mais ce n'est qu'une impression.
Regarde plutôt l'explication de Obsidian qui fourni la solution à ce problème.

Au passage, tu as le même problème avec les variables c, d et matrice

[Obsidian]

peut-être g trouvé la sortie pour le souci sur la ligne *sing=0. En fait, ds la declaration des variables ds la fonction main, javais initializé *sing=0. Quand je met seulement int *sing; ce souci est enlevé.

NON !

D'abord, ce n'est pas un warning mais une authentique erreur à l'exécution. Ensuite, en initialisant le pointeur à zéro, tu écrivais à un endroit incorrect mais tu pouvais le localiser et tu étais à peu près sûr que ton programme s'arrêterait. À présent, ton pointeur contient une valeur indéfinie, et si ton programme ne plante plus, c'est qu'il écrit n'importe où au sein d'un segment de mémoire qui t'appartient.

C'est encore pire !

[wintercold]

Je suis d'accord avec vous mais désolé le program est tjrs "bugged" sur la ligne

Code :

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

for (sum=0.0,i=k;i<=N;i++) sum += (A[i][k])*(A[i][k]);

[diogene]

Apparemment, A a pour dimensions NxN. Pourquoi (si c'est le cas) la condition i<=N ?

[wintercold]

Apparemment, A a pour dimensions NxN. Pourquoi (si c'est le cas) la condition i<=N ?

tt à fait, il fo enlever "=". le bug passe maintenant sur

Code :

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

c[k]=sigma*A[k][k];

Je ne vois rien qui va pas sur cette ligne.

[diogene]

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
void main()//void)
{
    FILE *fpt;
    float *c=0, *d=0;
...
	  qrdcmp(matrice1, N, c, d, sing);
...
 
void qrdcmp(float **A, int N, float *c, float *d, int *sing)
{
...
	for (k=0;k<N-1;k++)
		{
...
		if (scale == 0.0)
			{
			*sing=1;
			c[k]=d[k]=0.0;
			}
		else
			{
...
			A[k][k] += sigma;
			c[k]=sigma*A[k][k];
			d[k] = -scale*sigma;
...

Normal, c et d ne sont pas des tableaux mais des pointeurs.
Le programme aurait pu planter également sur c[k]=d[k]=0.0;

[wintercold]

Je vous remercie beaucoup. tout est éclairé.