Discussion :

[Suntory]

Bonjour,

Je souhaite lire mon fichier all-out2.dat qui est constitué de deux colonnes de valeurs (x,y) avec deux boucles for et un tableau bidimensionnelle.
Je voulais faire : si x==15 écrit moi les valeurs x et y dans le fichier haut.
Mon code créer bien le fichier haut mais il est vide et je ne vois pas mon erreur.
J'aurais besoin d'un coup de pouce s'il vous plaît.
Merci

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
		if(t==EndTime){
			FILE * fy = fopen("all-out2.dat", "r");
			FILE * fh = fopen("Haut","w");
	  		for (int i=0;i<N;i++){
	   		 for (int j=0;j<2;j++){
	      			field[i][j]=0.;
	      			fscanf(fy,"%lf",&field[i][j]);
				if(field[i][0]>=(15.)){
					fprintf(fh, "%lf %lf %lf\n",field[i][0], field[i][1], 
							interpolate (u.x, field[i][0], Heau/2.));
					fflush(fh);
					//fclose(fh);
				}
	    		}
	 		}
		}

De plus je ne sais pas à l'avance si x==15.00000 mais plutôt 15.02 ....ou autre. Du coup j'ai mis un seuil >=15..
Est-ce qu'il y a une autre solution pour ce genre de problème ? car 0.<=x<=20.
Pour avoir la valeur la plus proche de 15. ?

[pjlan]

Bonjour,

Plusieurs remarques :
- Vous devriez tester les retour de fopen() avant de continuer tout traitement
- Il faudrait fermer les fichiers une fois qu'il ne sont plus utilisés
- Les doubles se comparent avec un epsilon, en principe la constante DBL_EPSILON dans l'entête float.h. dans votre cas vous pouvez faire une fonction qui gére un pas plus sensible.

Code :

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

1
2
3
4
 bool double_egal(double a,double b,double epsilon)
{
   return fabs(a - b) < epsilon;
}

- Je crois (à vérifier) que pour fprintf() c'est '%f' même pour un double (ceci dit ça ne devrait pas avoir d'impacts les deux devraient être acceptés)


Dans l'absolu je ne vois pas pourquoi il n'y aurai rien du tout dans le fichier mais le plus probable est que la condition n'est jamais vérifiée...

[Xelland]

Dans l'absolu je ne vois pas pourquoi il n'y aurai rien du tout dans le fichier mais le plus probable est que la condition n'est jamais vérifiée...

il suffit de commenter la condition
si on écrit dans le fichier, c'est que c'est la condition qui ne fonctionne pas

[Sve@r]

Bonjour

Mon code créer bien le fichier haut mais il est vide et je ne vois pas mon erreur.

Et tu n'as pas pensé à rajouter un printf() de trace quand tu écris le fichier ? C'est quand-même un peu l'évidence quoi. Le fichier ne s'écrit pas alors tu rajoutes un truc à cet endroit et tu vois si le truc s'exécute quoi. Tiens, par exemple printf("%f\n", field[i][0]) voir s'il est supérieur à 15...

[Suntory]

re-bonjour,
merci de vos reponses ;

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
			//double (*field)[2] = malloc(N * 2 * sizeof(double));
			double field[N][2];
			FILE * fy = fopen("all-out2.dat", "r");
			FILE * fh = fopen("Haut","w");
	  		for (int i=0;i<N;i++){
	   		 for (int j=0;j<2;j++){
	      			field[i][j]=0.;
	      			fscanf(fy,"%f",&field[i][j]);
				//if(field[i][0]>=(L0-L0/4.)){
					fprintf(fh, "%f %f %f\n",field[i][0], field[i][1], 
							interpolate (u.x, field[i][0], Heau/2.));
					fflush(fh);
					fclose(fh);
					fclose(fy);
				//}
				//free(field[i][j]);
				//field[i][j]=NULL;
	    		 }
	 		}
			free(field);
			//field=NULL;

mon probleme vient de l'allocation de memoire mais je ne sais pas comment resoudre cette erreur

Code :

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

double free or corruption (!prev)

[Bousk]

D'où sort N ?
Il manque la vérification que l'ouverture des fichiers a fonctionné.
Pourquoi les fichiers sont fermés à chaque itération ?

[Sve@r]

mon probleme vient de l'allocation de memoire mais je ne sais pas comment resoudre cette erreur

Code :

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

double free or corruption (!prev)

Pas de free() sur un truc qui na pas été alloué par malloc() ou une fonction dérivée (realloc(), calloc()).

Pourquoi les fichiers sont fermés à chaque itération ?

Ou du moins pourquoi ne sont-ils pas réouverts au moment de les lire ou de les écrire...

D'où sort N ?

C'est marrant comment au fur et à mesure que le topic avance on remonte petit à petit à l'origine des erreurs bribes par bribes. Parce qu'elle n'y était pas cette ligne dans son premier post. Et on est tous partis du principe que le tableau était correct.