Discussion :

[booby]

Bonjour à tous,

je travaille actuellement sur un programme C qui prends en entrée un fichier binaire issu de fortran (créé avec un write(... form = 'unformatted')) qui contient un tableau de flottants de dimension 9000*9000.

Je me suis renseigné ici, ici et là mais j'ai du mal à voir ce que je vais obtenir une fois le fichier ouvert:
une fois que j'aurai ouvert le fichier avec un:

Code :

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

1
2
FILE * fichier;
fichier = fopen( "monfichier", "rb");

j'aurai bien un flux se rapportant au fichier concerné, sachant que je veux recuperer des flottants uniquement y a-t-il possibilité de les récupérer directement (j'entends par là sans avoir à passer par une chaine de caractère à caster en flottants ensuite)?

j'ai pensé utiliser fseek, celà vous semble-t-il correct?

Merci d'avance.

[zooro]

Bonjour,
Ton fichier ne contient pas des chaîne de caractères, ou des flottants, ou autres. Il contient des octets.
Si tu es certain que ces octets correspondent à des flottants, et que la représentation des flottants est identique en C et en Fortran (désolé, je ne connais pas le Fortran), alors tu peux lire les X octets correspondants à un flottant, et les stocker directement dans ta variable de type flottant.

[booby]

question bête comment connaît-on le nombre d'octets sur lequel est codé une variable?

Donc je peux utiliser fseek?

Code :

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

fseek( FILE *monfichier, long le_nombre_d'octets, int position);

ça me chiffone un peu fseek retourne un int non?

le code fortran sort des valeurs de ce type: -0.100009274483E+01
Un flottant n'est pas aussi précis, le type double suffira-t-il?

[homeostasie]

question bête comment connaît-on le nombre d'octets sur lequel est codé une variable?

Utilise la fonction sizeof qui retourne le nombre d'octets.
ex: octets = sizeof(double) -> 8 octets dans mon implémentation

[Jean-Marc.Bourguet]

Il faut que tu saches si tes flottants sont des simple (REAL*4 =C float) ou des double (REAL*8=C double) précisions (je ne plus sur de la syntaxe fortran). Si l'écriture a eu lieu sur la même machine, tu vas pouvoir lire ton tableau en faisant quelque chose du genre

Code :

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

1
2
3
4
5
6
7
int rowcount = 9000;
int colcount = 9000;
typedef double real; /* ou typedef float real; */
real *table = malloc(rowcount*colcount*sizeof(real));
if (fread(table, sizeof(real), rowcount*colcount, fp) != rowcount*colcount) {
   /* traitement de l'erreur */
}

Et alors tu peux accéder à tes éléments avec: table[c*rowcount+l] (note Fortran stocke les matrices transposées par rapport à l'usage C).

[zooro]

Donc je peux utiliser fseek?

Attention, fseek ne sert pas à lire des données dans le fichier, mais à s'y déplacer !

[homeostasie]

ça me chiffone un peu fseek retourne un int non?

Prototype de la fonction:

Code :

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

int fseek( FILE *stream, long depl, int posit);

Valeurs renvoyées:
-> 0 si déroulement correct
-> sinon une valeur non nulle est renvoyée. Exemple d'erreur: le flux n'a pas été ouvert.

Ca te permet de tester la réussite de ton positionnement dans le fichier et non de te retourner ton float ou double etc...
Il te faudra faire une lecture avec fread().

[booby]

ok, merci pour vos réponses, je teste ça et je vous tiens au courant.

[Emmanuel Delahaye]

question bête comment connaît-on le nombre d'octets sur lequel est codé une variable?

En C, on ne parle pas d'octets mais de char. Un char fait au moins 8 bits. Il fait exactement CHAR_BIT (défini dans <limits.h>).

sizeof retourne la taille d'un objet ou d'un (type) en nombre de char.

En combinant astucieusement ces informations, tu peux savoir combien de bits il y a dans un objet ou un type.

Mais attention, à part dans les [[un]signed] char, tous les bits ne sont peut être pas utilisables. Une certaine combinaison de bits pourrait très bien être une 'trap representation' et donc être invalilde. C'est pourquoi la plage des valeurs possibles d'un type est définie par des macros telles que INT_MIN, ULONG_MAX, toutes définis dans <limits.h>, indépendamment du nombre de bits d'un type.

[booby]

Code :

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

1
2
3
4
5
6
7
int rowcount = 9000;
int colcount = 9000;
typedef double real; /* ou typedef float real; */
real *table = malloc(rowcount*colcount*sizeof(real));
if (fread(table, sizeof(real), rowcount*colcount, fp) != rowcount*colcount) {
   /* traitement de l'erreur */
}

Et alors tu peux accéder à tes éléments avec: table[c*rowcount+l] (note Fortran stocke les matrices transposées par rapport à l'usage C).

ok, en regardant d'un peu plus pres pour adapter ça à mon code, je me suis aperçu que quelques subtilités m'échappent ;-):

sauf erreur de ma part, table est censé receuillir les nombres lus dans le fichier. Ce que je ne comprends pas c'est la ligne de test qui suit l'allocation mémoire:

Code :

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

if (fread(table, sizeof(real), rowcount*colcount, fp) != rowcount*colcount)

que teste-t-on exactement?
il me semble qu'à ce moment table est "vide" non? (je ne pense pas avoir compris comment le remplir :-( )

[zooro]

Code :

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

if (fread(table, sizeof(real), rowcount*colcount, fp) != rowcount*colcount)

La fonction fread lit "rowcount*colcount" éléments (de taille "sizeof(real)" chacun) dans le fichier décrit par "fp" et les stocke dans le tableau "table".
Le test sert à vérifier que tous les éléments attendus ont bien été lus.

[homeostasie]

que teste-t-on exactement?
il me semble qu'à ce moment table est "vide" non? (je ne pense pas avoir compris comment le remplir :-( )

On teste le retour de la fonction fread() au nombre d'éléments lus. Faut jeter un coup d'oeil au prototype de la fonction!

C'est l'équivalent de cela:

Code :

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

1
2
3
4
5
6
7
8
 
int octetsLus = 0;
 
octetsLus = fread(table, sizeof(real), rowcount*colcount, fp);
if (octetsLus != rowcount*colcount)
{
   /*traitement de l'erreur*/
}

[zooro]

C'est l'équivalent de cela:

Code :

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

1
2
3
4
5
6
7
8
 
int octetsLus = 0;
 
octetsLus = fread(table, sizeof(real), rowcount*colcount, fp);
if (octetsLus != rowcount*colcount)
{
   /*traitement de l'erreur*/
}

Euh... ben en fait non, ou alors le nom de ta variable est mal choisi

fread et fwrite renvoient le nombre d'éléments correctement lus ou écrits (et non pas le nombre d'octets).

[Emmanuel Delahaye]

Euh... ben en fait non, ou alors le nom de ta variable est mal choisi

champsLus ...

[homeostasie]

Euh... ben en fait non, ou alors le nom de ta variable est mal choisi

Oui...Le nom de la variable est mal choisi! Sorry.

[booby]

Code :

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

if (fread(table, sizeof(real), rowcount*colcount, fp) != rowcount*colcount)

La fonction fread lit "rowcount*colcount" éléments (de taille "sizeof(real)" chacun) dans le fichier décrit par "fp" et les stocke dans le tableau "table".
Le test sert à vérifier que tous les éléments attendus ont bien été lus.

ok, donc cette ligne fait deux opérations d'un coup!... vachement bien fait tout de même.

Merci encore de vous etre penchés sur mon problème, j'essaie et vous tiens au courant.

[Emmanuel Delahaye]

ok, donc cette ligne fait deux opérations d'un coup!... vachement bien fait tout de même.

C'est fread() qui fait ça (lecture par blocs), et je ne vois pas comment elle pourrait faire autrement...

[booby]

ok, tres bien, ça marche! merci à tous pour votre aide! chapeau!