Problème écriture dans un fichier

Version imprimable

Bonjour, j'ai un souci avec l'écriture dans un fichier de mes résultats.
C'est un projet qui utilise la fonction tridag pour la résolution d'une matrice tridiagonal.
Voici le code:

Code:

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
51
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include "tridag.h"
 
 
main()
{                                   
      double k,a,b,L,T0,Te,t,teta0;                 
      double A[100],B[100],C[100],D[100],u[100];     
      unsigned long N;
      int i,j;
      FILE *fp;                                                      
      fp=fopen("RESULTAT", "w");    
 
      printf("Entrez les valeurs de k,a,b,L,T0,Te et N\n");          
      scanf("%lf%lf%lf%lf%lf%lf%d",&k,&a,&b,&L,&T0,&Te,&N);
 
      t=(10*(2*a+2*b)*L*L)/(k*a*b);                            
      teta0=(T0-Te)/Te;
 
 
 
      for (i=1;i<N;i++){                                 
                B[i]=2+(1./N)*(1./N)*t;
                        }
                          B[N]=1+(1./N)*(1./N)*t;   
 
 
      for (i=2;i<=N;i++){                                 
                A[i]=-1;
                        }
      for (i=1;i<=(N-1);i++){                                 
                C[i]=-1;
                        }
      D[1]=teta0;
      for (i=2;i<=N;i++){                                 
                D[i]=0;
                        }
 
     tridag(A,B,C,D,u,N); 
 
 
     for (i=1;i<=N;i++){                                   
                        u[i]=Te*u[i]+Te;
                        //printf("T[%d]=%lf\n",i,u[i]);
                        fprintf(fp,"%d      %lf\n",i,u[i]);           
                            }
                        //getch();
                        fclose(fp);                                
          }

Voici également la fonction tridag:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
void tridag(double a[], double b[], double c[], double r[], double u[],
	unsigned long n)
{
	unsigned long j;
	double bet,*gam;
 
	gam=(double*)malloc(n*sizeof(double));
	if (b[1] == 0.0) printf("Error 1 in tridag");
	u[1]=r[1]/(bet=b[1]);
	for (j=2;j<=n;j++) {
		gam[j]=c[j-1]/bet;
		bet=b[j]-a[j]*gam[j];
		if (bet == 0.0)	printf("Error 2 in tridag");
		u[j]=(r[j]-a[j]*u[j-1])/bet;
	}
	for (j=(n-1);j>=1;j--)
		u[j] -= gam[j+1]*u[j+1];
 
}

Pour les valeurs a entrer, dans l'ordre:

Code:

Je ne comprend pas, peut-ètre que j'utilise mal les fopen, fprintf...
Enfin le fait est que l'affichage a l'écran s'effectue, j'ai les bons résultats, mais a l'écriture, la console s'arrète.
Merci bcp d'avance.

Pour info je tourne sous xp, j'utilise dev cpp, et c'est un projet de physique, je ne suis pas informatitien donc désolé d'avance si l'erreur est bête.

13/06/2007, 13h33
GyZmoO

Quelques remarques :

Dans le main :
- int main() et non pas main()
- pas de test sur l'ouverture du fichier
- scanf a éviter ----> préférer plutot fgets + la ou les fonctions de conversion qui vont bien
- tu demandes a l'utilisateur de rentrer N, qui est la borne de tes boucles, si l'utilisateur rentre 101 par exemple, t'es fichu car tes tableaux sont déclarés comme ça : B[100] . . .

Dans la fonciton :

- pas de test sur malloc
- pas besoin de caster le retour d'un malloc

Code:

1
2
 
                        fprintf(fp,"%d      %f\n",i,u[i]);

Code:

1
2
 
                        fprintf(fp,"%d      %g\n",i,u[i]);

13/06/2007, 19h57
bob_warrior89

re

[-mod- Merci de relire les regles du forums. La langue d'expression est le français et non le borborygme. Merci de te relire et de te corriger]

Tout d'abord merci pour ces premières réponces,
J'ai éssayé de changer %lf en pourcent g ou f, mais ça n'as rien changer, la console plante toujours.
Et pour ce que tu me suggère Gyzmo, comme par exemple fgets, je ne sais pas m'en servir.
Apriori je pense que c'est à cause du sous programme que sa ne fonctionne pas, en effet si je place le fprintf avant la fonction tridag sa marche( pour écrire B[] par exemple), dès que je place sa après, rien, même un texte ne sera pa écrit.
13/06/2007, 20h53
GyZmoO

Pour fgets : http://man.developpez.com/man3/fgets.3.php
Tu trouveras ton bonheur
13/06/2007, 21h49
bob_warrior89

C'est bon j'ai trouvé, c'est la mémoire alloué par le malloc qui n'était pas suffisante, et sans trop savoir pourquoi j'ai rajouté un facteur dix et ça a fontionner.
Code:

gam=(double*)malloc(n*sizeof(double)*10);

Merci d'avoir cherché!
++
13/06/2007, 21h56
droggo

Mau,

Citation:

Envoyé par bob_warrior89

C'est bon j'ai trouvé, c'est la mémoire alloué par le malloc qui n'était pas suffisante, et sans trop savoir pourquoi j'ai rajouté un facteur dix et ça a fontionner.
Code:

gam=(double*)malloc(n*sizeof(double)*10);

Merci d'avoir cherché!
++

à propos du malloc, où est le free ?
13/06/2007, 21h57
Emmanuel Delahaye

Citation:

Envoyé par bob_warrior89

C'est bon j'ai trouvé, c'est la mémoire alloué par le malloc qui n'était pas suffisante, et sans trop savoir pourquoi j'ai rajouté un facteur dix et ça a fontionner.

Code:

gam=(double*)malloc(n*sizeof(double)*10);

Programmation au hasard = comportement indéfini. Tu n'as rien résolu du tout.
13/06/2007, 23h00
souviron34
En plus de tout ce qu'on a dit au dessus, je signale que en C les tableaux commencent à 0...

Donc :
Code:

1 2 3 4 5 6 7 8 for (j=2;j<=n;j++) { gam[j]=c[j-1]/bet; bet=b[j]-a[j]*gam[j]; if (bet == 0.0) printf("Error 2 in tridag"); u[j]=(r[j]-a[j]*u[j-1])/bet; } for (j=(n-1);j>=1;j--)
Est FAUX et entraîne un dépassement de mémoire (il n'y a pas de gam[n], ni vraisemblement des autres tableaux... ). Et en plus ne prend pas en compte le premier point...

tableau[1] réfère au deuxième élément...
14/06/2007, 10h23
bob_warrior89

Désolé, mais je ne sais pas se qu'est un free.
Ensuite, pour tout se qui est du malloc ou du reste, en fait mon prof m'a juste donné la fonction et il m'a dit de l'utiliser sans expliquer son fonctionnement.
Donc je ne savais pas trop si je devait la modifier ou pas.
Mais du coup, est-ce que les résultat que j'obtient sont éronés?
14/06/2007, 12h23
Emmanuel Delahaye
Citation:

Envoyé par bob_warrior89

Désolé, mais je ne sais pas se qu'est un free.
Ensuite, pour tout se qui est du malloc ou du reste, en fait mon prof m'a juste donné la fonction et il m'a dit de l'utiliser sans expliquer son fonctionnement.
Donc je ne savais pas trop si je devait la modifier ou pas.
Mais du coup, est-ce que les résultat que j'obtient sont éronés?

J'ai quelques doutes sur les qualités pédagogiques de l'enseignement que tu reçois...

En principe, tu ne dois pas écrire une ligne de code que tu ne maitrises pas parfaitement. Pour la maitriser, il n'y a pas de secrets :
- Explications du prof
- Cours, livre de référence
- Exercices personnels sur le sujet (et uniquement celui-ci, en utilisant éventuellement d'autres connaissances maitrisées).
C'est mon point de vue de formateur en informatique...
14/06/2007, 15h37
souviron34

Citation:

Envoyé par bob_warrior89

...Mais du coup, est-ce que les résultat que j'obtient sont éronés?

OUI à cause de toutes les remarques ci-dessus...

Et re-OUI encore une fois à cause des indices des tableaux..