Conversion d'un code FORTRAN en C

**stage75** · 16/09/2010, 10h10

Bonjour,
D'abord je suis pas informaticien mais j'utilise les outils informatique dans mon boulot, je travaille normalement sur fortran mais là je dois faire un petit programme sur C et je connais rien dans le monde de C, est ce que ya quelqun qui peut me convertir mon programme fortran en C :

Code fortran :

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

1
2
3
4
5
6
7
8
9
OPEN (UNIT=50,FILE='sourceO.txt') 
    		DO 1001 S=1,280
    		READ (50,*)LTEST, ATEST
	IF ((L.NE.LTEST).AND.(NUMSPC.EQ.3)) GOTO 1001
	ATERM = 0.0
	BTERM = 0.0
	GOTO 1002
1001	CONTINUE
1002	CLOSE(50)

merci

**Sve@r** · 17/09/2010, 21h52

Ben déjà si tu nous disait ce qu'il fait ce programme ? On sent bien qu'il y a un truc de lecture de fichier mais c'est pas parce qu'on est bons en C qu'on est forcément bons en Fortran !!!

**souviron34** · 18/09/2010, 15h14

Code FORTRAN :

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

1
2
3
4
5
6
7
8
9
OPEN (UNIT=50,FILE='sourceO.txt') 
    		DO 1001 S=1,280
    		READ (50,*)LTEST, ATEST
	IF ((L.NE.LTEST).AND.(NUMSPC.EQ.3)) GOTO 1001
	ATERM = 0.0
	BTERM = 0.0
	GOTO 1002
1001	CONTINUE
1002	CLOSE(50)

Deviendrait (si transformée en procédure) :

Code C :

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
#include <stdlib.h>
#include <stdio.h>
 
 
void MaFonction ( int L, int NUMSPC, double *ATERM, double *BTERM )
{
FILE *Fichier= NULL ;
char Ligne[1000] ;
int i ;
int LTEST ;
float ATEST ;
 
Fichier = fopen ( "sourceO.txt", "r" );
if ( Fichier != NULL )
  {
       while ( fgets ( Ligne, 999, Fichier) != NULL )
         {
            if ( sscanf ( Ligne,  "%d %g", &LTEST, &ATEST ) == 2 )
              {
                  if ( ! ((LTEST != L) && (NUMSPC == 3)) )
                    {
                        *ATERM = 0.0 ;
                        *BTERM = 0.0 ;
                        break ;
                    }
              }
         }
 
       fclose ( Fichier );  
  }

Fonction qui doit être appelée comme :

Code C :

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

MaFonction ( L, NUMSPC, &ATERM, &BTERM );

Par contre :

a) il est possible d'optimiser le if...
b) il serait souhaitable de passer le nom du fichier en paramètre...

du style :

Code C :

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

MaFonction ( "source0.txt", L, NUMSPC, &ATERM, &BTERM );

Avec une définition comme :

Code C :

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

1
2
3
4
5
6
7
8
9
10
11
void MaFonction ( char *Filename, int L, int NUMSPC, double *ATERM, double *BTERM )
{
FILE *Fichier= NULL ;
char Ligne[1000] ;
int i ;
int LTEST ;
float ATEST ;
 
Fichier = fopen ( Filename, "r" );
if ( Fichier != NULL )
..

**stage75** · 22/09/2010, 16h37

Pour répondre à ta question:
ce que je voudrais faire dans ce programme c'est d'ouvrir un tableau qui s'appelle "sourceO.txt", ce tableau il contient deux colonnes que j'ai nommé (LTEST: première colonne et ATEST: deuxièmement colonne) et chaque valeur de la deuxieme colonne correspond au maille de la premiere colonne et le but c'est d'injecter chaque valeur de la deuxieme colonne dans le numero de maille de la première colonne voilà j'espere que j'ai bien expliquer le probleme

**Sve@r** · 22/09/2010, 22h37

Envoyé par stage75

et le but c'est d'injecter chaque valeur de la deuxieme colonne dans le numero de maille de la première colonne voilà j'espere que j'ai bien expliquer le probleme

Presque. C'est quoi le séparateur qui sépare la colonne 1 de la colonne 2 ?
Les nombres sont-ils des entiers ou ont-ils une virgule ?
C'est quoi un "maille" ??? T'aurais-pu rajouter le début de ton fichier qu'on voit mieux.

Bon, déjà le code pour ouvrir le fichier et le lire ligne à ligne

Code c :

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
#include <stdio.h>
#include <errno.h>
 
int main()
{
    FILE *fp;
    char ligne[81];
 
    fp=fopen("sourceO.txt", "r");
    if (fp == NULL)
    {
         perror("Erreur ouverture sourceO.txt");
         exit(1);
    }
 
    while (fgets(ligne, 81, fp) != NULL)
    {
         printf("Ligne lue: %s", ligne);
    }
 
    fclose(fp);
}

Bon, j'ai tapé à capella donc peut-être que j'ai oublié un include mais j'espère que tu vois l'idée...

**stage75** · 23/09/2010, 10h20

Merci Sve@r,

EN FAIT MON JE SUIS PAS INFORMATICIEN C'EST POUR çA J'EXPLIQUE MAL, MOI J4UTILISE JUSTE LES OUTILS INFORMATIQUE POUR FAIRE DE LA SIMULATION

1 maille ça vient du maillage

par exemple je devise mon domaine de calcul à 200 petit carreau et moi je veux injecter des termes sources que dans des cellules définies, c'est pour ça mon tableau il contient deux colonnes la 1ere colonnes contiens des entiers (1, 2, 3,....100) et en face dans la 2éme colonne j'ai les valeurs de termes sources corresponds à chaque cellule (1E-12, 2E-11,...2.4E-10) et moi ce que je veux c'est de faire un petit programme qui permis de lire ce tableau et de dire a ce programme que la valeur 1E-12 il faut la mettre dans la cellule 1 et la valeur 2E-11 faut la mettre dans la cellule 2....et la valeur 2.4E-10 faut la mettre dans la cellule 100.. voilà, voiçi une partie de mon programme fortran

Code fortran :

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
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
                                                                   C
C----------------------------------------------------------------------C
C
#include "IMPLICIT.INC"
C
#include "SIZE.INC"
#include "SPECSZ.INC"
C
C----------------------------------------------------------------------C
C     COMMON BLOCKS...                                                 C
C----------------------------------------------------------------------C
C
#include "BOUNDS.INC"
#include "GRID.INC"
#include "FLOVEL.INC"
#include "FLOPRO.INC"
#include "FLPEXV.INC"
#include "FLPGEO.INC"
#include "SPCONS.INC"
#include "FTURB.INC"
#include "MPHASE.INC"
#include "MPPROP.INC"
#include "MPVOF.INC"
#include "CSHEAR.INC"
#include "DNITER.INC"
#include "TDFLOW.INC"
#include "USFUNS.INC"
#include "POLAR.INC"
#include "GAPS.INC"
#include "AREAS.INC"
#include "LOCATE.INC"
#include "LOCATV.INC"
#include "LCTYPE.INC"
#include "LBOUND.INC"
#include "WSFLOW.INC"
#include "TDTEMP.INC"
#include "CHANEL.INC"
#include "NCHANL.INC"
C
C----------------------------------------------------------------------C
C     ARGUEMENT TYPE DECLARATIONS...                                   C
C----------------------------------------------------------------------C
C
      INTEGER IVBLE
      INTEGER IPHASE
      REAL    ATERM
      REAL    BTERM
C
C----------------------------------------------------------------------C
C     LOCAL VARIABLE TYPE DECLARATIONS...                              C
C----------------------------------------------------------------------C
C
      INTEGER NUMSPC
      INTEGER NUMVOF
      INTEGER LTEST,S
      REAL  ATEST,A_N
      REAL  A_T,A_T2
      REAL  A_V
 
SINGLE PHASE FLOW :
C
 
     	OPEN (UNIT=40,FILE='source.txt') (#source.txt c'est mon fichier tableau#)
 
     		DO 1005 S=1,276  (# S sont le nombre des ligne (de 1 jusqu'a 276#)
     		READ (40,*)B, A_T (# B c'est le nom de la 1ere colonne et A_T c'est le nom de la 2eme colonne#)
 
	IF (L.NE.B) GOTO 1005  (# L c'est le nom du maille (chaque petite cellule s'appel L #)
 
	ATERM = A_T*YACOB(L)   (# ATERM represente le parametre terme source#)
 
	GOTO 1006
 
1005	CONTINUE

Conversion d'un code FORTRAN en C

C

Vue hybride

Discussions similaires

Partager

Partager