Matrices sparse avec Boost

Bonjour!

J'utilise la library boost pour creer une matrice sparse 15 000 x 3000. En faisant tourner le programme, la memoire utilisee dans le programme augmente a toute vitesse... Aussi, j'aurais quelques questions concernant la construction d'une matrice sparse. Mon constructeur est le suivant:

boost::numeric::ublas::mapped_matrix<double> m(15000, 3000);

Pendant la simulation, je la remplis en utilisant les commandes:

Code:

---

m(i, j) = m(i, j) + 1;

---

J'ai aussi vu qu'il existe une fonction

Code:

---

m.insert_element (i, j, t);

---

ou t est l'element a inserer.

1) Est-ce que l'une ou l'autre des operations est plus couteuse?

2) Faut-il specifier le nombre estime d'elements non-nuls dans le constructeur de mapped_matrix?

3) Quelle est la difference entre mapped_matrix et compressed_matrix de la library boost?

Merci d'avance pour les suggestions!!

La classe 'mapped_matrix' utilise un map entre le la valeur de 'i*size+j' et le élément que tu spécifie.

L'opérateur () et insert_element sont très similaires, tout les deux appellent 'insert' dans le map.

La mémoire utilise peut être estimé avec la formule donné ici:

http://stackoverflow.com/questions/7...sage-of-stlmap

Pour compressed_matrix la méthode utilisé pour mémoriser est différent:

http://netlib.org/linalg/html_templates/node91.html