Discussion :

[doctorant_en_gallère]

bonjour a tous,

je suis doctorant et j'utilise matlab pour résoudre des équation de transfères de masse grâce a la loi de fick.

mon problème vient du fait que je réalise du calcul itératif qui stock mes vecteur (champs de température (100 noeud) ou champ de concentration (idem)) dans une matrice de sauvegarde à chaque pas de temps. celle-ci sert a afficher mes résultats en fin de calcul.

le problème ce situe dans le fait que mes matrices de stockage de données sont très volumineuse et ralentisse mon calcul (temps de calcul quasi exponentiel avec le temps de simulation).

je cherche donc a sauvegardé mes matrices de donnée vers un fichier externe types .txt par exemple. pour remplir ce fichier texte tout les "n" nombre de cycle et vidangé mes matrices de sauvegarde.

l'objectif étant a la fin du calcul de pouvoir rechargé ces fichier contenant la concaténation cyclique de mes matrice de données.

mon problème c'est que avec fpintf j'arrive a stocker les donnée mais je perds l'information matriciel (plus de différenciation ligne/colonne). mes donnée sont donc inutilisable.

avec save pas réussi a concaténé dans le fichier .mat les données de la matrice lors des sauvegarde successive.


voici un morceau de mon code de sauvegarde externe qui est au cœur de ma boucle temporelle qui réalise le calcul de simulation.

CODE :

Code :

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

1
2
3
4
5
6
7
8
9
10
	if mod(i,100)==0
        disp(sprintf('Temps de simulation : %8.2f s',t(i+1)))
        toc
 
        fid = fopen('Ts.txt','a+');
        fprintf(fid,'%12.8e %12.8e\n',Ts);
        fclose('all');
 
         Ts=[];
        end

[Invité]

Bonsoir,

Je te revoie à un sujet qui a été traité il n'y a pas si longtemps
http://www.developpez.net/forums/d10...fichier-texte/

[Jerome Briot]

Au sujet de l'écriture dans le fichier, je sortirais les fonction d'ouverture et de fermeture de la boucle :

Code :

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

1
2
3
4
5
6
7
8
9
10
 
fid = fopen(...,'wt');
 
for n = ...
 
   fprintf(...)
 
end
 
fclose(fid);

Ensuite, il faudrait en savoir un peu plus sur la quantité de données que tu traites. Quel est le nombre de valeurs stockées au final dans le fichier texte ?

As-tu envisagé la pré-allocation de mémoire ?
Voir la : Qu'est-ce que la préallocation de mémoire ?

[doctorant_en_gallère]

merci pour ses réponse si rapides,

je n'ai pas encore eu le temps de tester vue que j'étais en cours ce matin.
mais pour répondre a la question des quantité de données!!! c'est ça le problème!

a l'heure actuelles mon programme plante a cause d'erreur mémoire au bout d'une 10ène d'heures simuler car mes matrice de stockage qui pour l'instant n'était pas externalisé sature.

soit 100*3600/2*10= 1 800 000 données et ceci dans le cas de seulement 100 noeuds sauvegardé toute les deux seconde pendant 10h.

or dans l'absolus le nombre de nœuds peu augmentés le pas de temps de sauvegarde difficilement baissé et le temps de simulation être multiplier par 4 ou 5.

mon idée était de stoker mes donnée dans des variable tampon comme je le fait actuellement mais de les transférés a intervalle de temps régulier vers un fichier extérieur pour les vidés.
le but étant de ne pas avoir a écrire dans le fichier a chaque tour de boucle pour gagner du temps mais peu être est ce la mon erreur et devrais-je supprimé mes variable de stockage de donnée et écrire les résultat de chaque pas de temps directement dans un fichier externe.

je vais commencer par tester les solution que vous me proposé pour voir ce que ça change.

encore merci.

[doctorant_en_gallère]

en combinant la solution de Dut dans le liens de Winjerome j'ai fini par obtenir quelque chose comme ça :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
        fid = fopen('Ms.txt','at');
        fmt = [repmat('%12.8e ',1,length(M)) '\n'];
        fprintf(fid,fmt,Ms);
 
        fid = fopen('M1s.txt','at');
        fmt = [repmat('%12.8e ',1,length(M1)) '\n'];
        fprintf(fid,fmt,M1s);
 
        fid = fopen('Xms_sav.txt','at');
        fmt = [repmat('%12.8e ',1,length(Xms)) '\n'];
        fprintf(fid,fmt,Xms_sav);
fclose('all')

etc au coeur de ma boucle temporelle et seulement tout les X pas de temps.

cella créer bien des fichier text avec les bonnes données et une structure ligne colonne cohérente.

je pense donc que ça répond a ma question... sauf que quand j'appelle mon fichier text en fin de programme pour récupéré ma matrice complète avec la fonction fread je récupère des données incohérente :

Code :

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

1
2
3
4
5
6
7
8
9
10
 
        fid = fopen('Ms.txt','rt')
        Ms=fread(fid);
 
        fid = fopen('M1s.txt','rt')
        M1s=fread(fid);
 
        fid = fopen('Xms_sav.txt','rt')
        Xms_sav=fread(fid);
fclose('all')

j'imagine que c'est le même type de problématique c'est a dire un problème dans le format de lecture si quelqu'un a un tuyaux je serai donc preneur.

quoi qu'il en soit merci pour l'aide précieuse déjà apporté.

[Buide]

Essaie avec: fscanf.

Code :

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

doc fscanf

[doctorant_en_gallère]

merci a force d'acharnement j'ai finie par trouvé un code avec fscanf qui marche en effet très bien.