Discussion :

[RMATH]

Bonjour,
Après avoir discrétiser un problème sur FreeFem (j'ai pris laplacien avec les conditions au bord de Dirichlet) j'ai enregistré la matrice A et le vecteur B sous format txt afin de les utiliser sur scilab (pièce jointe), mais j'ai pas réussie à le lire

Code :

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

1
2
3
4
 
A=fscan(A,txt)
!--error 4 
Undefined variable: A

J'ai contacté un professeur, il m'a envoyé le script suivant mais je rencontre toujours le meme problème

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
 
function A=readFreeFemExport(f)
  u=mopen('f','r')
  t=stripblanks(mgetl(u,1));
  if t=="# Sparse Matrix (Morse)" then
    mgetl(u,2);
    sz=mfscanf(1,u,"%d %d %d %d\n");
    is_sym=sz(3)<>0;
    nel=sz(4);
    d=mfscanf(nel,u,"%d %d %lg");
    mclose(u)
    A=sparse(d(:,1:2),d(:,3),sz(1:2));
  else
    if execstr("n="+t,"errcatch")<>0 then
      error("Unknown file format")
    end
    A=mfscanf(n,u,"%lg")
    mclose(u)
  end
endfunction

voici ce qui s'affiche quand je l’exécute:

Code :

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

1
2
3
4
5
6
 
-->A=readFreeFemExport('A.txt')
 !--error 999 
mopen: Cannot open file f.
at line       2 of function readFreeFemExport called by :  
A=readFreeFemExport('A.txt')

quand j'utilise mopen

Code :

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

1
2
3
4
5
6
7
8
 
-->[fd,err]=mopen('A.txt','r')
 err  =
 
  - 2.  
 fd  =
 
    1.

j'ai essayé avec le format csv mais je trouve:

Code :

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

1
2
3
4
 
-->A=csvRead('A.csv');
                   !--error 999 
csvRead: A.csv does not exist.

aidez moi s'il vous plait !!

[Jerome Briot]

Dans readFreeFemExport, change :

Code :

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

u = mopen('f', 'r');

par :

Code :

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

u = mopen(f, 'r');

[RMATH]

Bonjour,
Merci Dut pour votre réponse.
J'ai changé mopen et voici ce qui s'affiche

Code :

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

1
2
3
4
 
-->U=readFreeFemExport(A)
                      !--error 4 
Undefined variable: A

je ne sais plus quoi faire??

[Jerome Briot]

Comme ceci :

Code :

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

U = readFreeFemExport("A.txt")

Si le fichier A.txt est bien dans le dossier courant de Scilab.

Sinon, il faut alors passer le chemin complet vers le fichier

Par exemple :

Code :

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

U = readFreeFemExport("D:\data\A.txt")

[RMATH]

Bonjour,
Merci infiniment Dut, ça marche !!
en fait je voulais savoir si j'enregistre la matrice sous format txt ou csv?? et pourquoi??
merci encore.

[Jerome Briot]

en fait je voulais savoir si j'enregistre la matrice sous format txt ou csv?? et pourquoi??

Je ne comprend pas bien ta question...

[RMATH]

Après la discrétisation sous FreeFem on peut sauvegarder la matrice sous format txt et csv...
Je voudrais savoir le quel vous me conseillez ? sachant que les matrices sont de grande taille.
Merci

[Jerome Briot]

Si les matrices sont de grandes tailles, il est peut être préférable de s'orienter vers un format binaire pour diminuer la taille des fichiers.

Mais tout dépend de ce que tu comptes faire avec ces fichiers

[RMATH]

Je veux résoudre le système Ax=b par la méthode multigrille.

[Jerome Briot]

Je veux résoudre le système Ax=b par la méthode multigrille.

Et quel est le rapport avec un quelconque format de fichier ?

[RMATH]

je n'ai pas bien compris, je voulais savoir si j'enregistre les fichiers sous format csv est ce que je gagne un peut de temps lorsque j'appelle le fichier sur scilab ou le format txt est mieux??? généralement les matrices sont de grandes tailles.

[Jerome Briot]

Le choix entre txt et csv ne devrait pas beaucoup influer sur le temps de calcul puisque généralement, dans le cas des matrices de grandes dimensions, c'est la résolution du système qui prend du temps (et non pas le pré- et post-traitement)

[RMATH]

Merci beaucoup Dut.
Mon problème est résolu !!
Bonne journée !!