Discussion :

[alain10025]

Bonjour,

Je travaille sur Matlab avec des matrices carrées ayant entre 32000 et 64000 lignes. J'utilise donc la fonction sparse.

Sauriez-vous comment on peut enregistrer une matrice sparse en tableau 3 colonnes (les deux premieres etant les coordonnées dans la matrice et la troisieme colonne la valeur à cette position), je cherche en fait un "inverse" de la fonction "spconvert" ...

Merci.

Alain

[Jerome Briot]

Essaie ceci :

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
A=sparse((rand(4)>.5).*rand(4))
 
[r,c]=find(A);
nz=nonzeros(A);
 
fid=fopen('sparse.txt','wt');
 
   fprintf(fid,'%d %d %.8f\n',[r c nz].');
 
fclose(fid);
 
% Juste pour la vérification
edit sparse.txt
 
D=load('sparse.txt','-ascii');
spconvert(D)

Malheureusement, cette solution risque d'être assez lente sur une vraie matrice creuse

[Jerome Briot]

J'ai quelque peu simplifié le code ci-dessus (mais il risque d'être gourmand en ressource)

[Jerome Briot]

Je viens de voir la réponse de Tim Davis (Mr Sparse Matrix) sur le CSSM :

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
A=sparse((rand(4)>.5).*rand(4))
 
[r,c,nz]=find(A);
 
fid=fopen('sparse.txt','wt');
 
   fprintf(fid,'%d %d %.8f\n',[r c nz].');
 
fclose(fid);
 
% Juste pour la vérification
edit sparse.txt
 
D=load('sparse.txt','-ascii');
spconvert(D)

C'est encore plus propre

[alain10025]

Merci beaucoup !

Je rajoute une solution qui m'a été donnée par un collègue :

Code :

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

1
2
[i j x] = find (A) ;
D = [i j x] ;

ou A est la matrice sparse et D le vecteur voulu.

Alain

[Jerome Briot]

Euhh, ça doit être la fatigue, mais je ne sais plus pourquoi je me suis mis en tête d'écrire le résultat dans un fichier

[Jerome Briot]

Juste par curiosité, à quoi cela va-t-il te servir ?

[alain10025]

Je travaille avec une matrice K 64000 par 64000 et j'utilise la commande suivante:
K(1:32000,1:32000)=K11;

K11 etant une matrice sparse 32000 par 32000 (un bloc donc de la matrice K)

Le problème est que cette commande demande un temps considérable à Matlab.

Il vaut mieux donc repasser par un tableau trois colonnes et modifier les indices du tableau correspondant à la matrice K11 et ajouter une derniere ligne:

64000 64000 1
K11(64000,64000)=0;

Alors la matrice K11 devient la matrice K en très peu de temps.

J'avoue ne pas être sûr de l'optimalité de la solution ...

[Jerome Briot]

Tu trouveras peut être quelques informations utiles sur le blog de Loren Shure : http://blogs.mathworks.com/loren/200...ces-in-matlab/