Allocation memoire pour grande matrice : tableau de tableau ?

Bonjour à tous, j'ai des questions concernant la façon dont c++ gère certaines allocation de mémoire.

J'ai besoin de créer une matrice de petits entiers relativement grande (16000*16000) et d'en calculer les puissances selon un algèbre Max+(je remplace l'addition par min et la multiplication par l'addition).

J'en suis pour l'instant à calculer la première matrice (puissance 1). Évidemment j'ai testé un short Tab[16000][16000] qui me donnait un stack overflow, du coup j'ai remplacé ça en inscrivant chaque case dans un fichier txt, mais les entrées sorties en mémoire sont réputés comme étant super lent.

J'aimerais donc garder la matrice totale en mémoire vive, j'ai remplacé les inscriptions dans le fichier txt par un truc dans ce genre :

Code:

---

1 2 3 4

short **Matr= new short*[16824]; for(i=0;i<16824;i++){ Matr[i]=new short[16824]; }

---

mais j'ai l'impression que le programme ne va pas plus vite, loin de là.

Les pointeurs pointent-ils bien vers la ram?

Y a-t-il un autre moyen pour calculer la totalité de ma matrice plus rapidement?

Y a-t-il un autre moyen que les threads pour utiliser la puissances de mes 4 procs en meme temps (pour l'instant un seul mouline pendant l'exécution) ?

Merci

boost.arrays me parait être une bonne solution mais je vois pas plus rapide que des pointeurs...

Bonjour,

Citation:

---

Envoyé par lyxthe

j'ai testé un short Tab[16000][16000] qui me donnait un stack overflow...

Code:

---

1 2 3 4

short **Matr= new short*[16824]; for(i=0;i<16824;i++){ Matr[i]=new short[16824]; }

---

---

ne placez tout simplement pas votre tableau dans la pile, mettez le dans le tas via un new short Tab[16000][16000] ou tout simplement dans une variable globale (éventuellement une variable de classe). Il y a à cela trois avantages, 1) vous n'aurez pas a 'fabriquer' celui-ci via la boucle précédente, 2) vous pourez l'utiliser en tant que tableau sans devoir faire vous même le calcul de l'accès à un élément, 3) les accès seront légèrement plus rapide avec une indirection en moins

Citation:

---

Envoyé par bruno_pages

Bonjour,



ne placez tout simplement pas votre tableau dans la pile, mettez le dans le tas via un new short Tab[16000][16000] ou tout simplement dans une variable globale (éventuellement une variable de classe). Il y a à cela trois avantages, 1) vous n'aurez pas a 'fabriquer' celui-ci via la boucle précédente, 2) vous pourez l'utiliser en tant que tableau sans devoir faire vous même le calcul de l'accès à un élément, 3) les accès seront légèrement plus rapide avec une indirection en moins

---

Dans l'état actuel, le tableau est déjà dans le tas, seul le pointeur qui en donne l'accès est sur la pile. Sinon T** t=new T[X][Y] ne compilera pas, il faudrait préciser une dimension.

Enfn, la bonne solution c'est blitz++, boost::ublas, ....

Citation:

---

Envoyé par Davidbrcz

Dans l'état actuel, le tableau est déjà dans le tas

---

:nono: actuellement il y a des vecteurs alloués dans le tas, pas un tableau ;)

Citation:

---

Envoyé par bruno_pages

:nono: actuellement il y a des vecteurs alloués dans le tas, pas un tableau ;)

---

Peut que je suis mal reveillé mais je ne vois pas ce que tu veux dire. Dans le code du PO on a un tableau de tableau alloué dans le tas, point.

PO >>J'ai oublié de te demander si ta ta matrice est creuse ?

un tableau à 1 dim est un vecteur, mais je n'ai pas été clair,il aurait fallut dire actuellement il y a des vecteurs alloués dans le tas, pas un tableau à deux dimentions

La première est peu danse, mais au fur est à mesure mon but est d'en calculer les puissance. A priori plus je vais avancé dans la série des puissances de la matrice plus elle sera danse. Du coup je suis pas certain que ce soit intéressant de la coder de façon particulière.

Je vais me concentrer plus particulièrement sur des soucis de performance dans cette discution : http://www.developpez.net/forums/d96...ds-appropriee/

En général, pour les matrices on conseille plutôt d'allouer un seul tableau de n*m éléments plus qu'un tableau de n tableau de m éléments.

Ensuite, on fait une interface d'accès avec l'opérateur ()(int abs, int ord), ce qui permet de structurer les éléments au sein de sa table sans se préoccupper de l'interface (chose impossible avec [][]).

Ici, ton tableau fait 512Mo. Je pense que tu peux l'allouer en ram contigüe avec un bête new short[16000*16000], sur une machine moderne ça devrait passer si elle n'est pas trop chargée.

Après, ce qui va être important, c'est l'agencement en mémoire de tes éléments, par rapport aux parcours que tu vas faire (en règle générale, il faut toujours minimiser la "distance" parcourue en mémoire lors de ton parcours).