Discussion :

[stchouan]

Bonjour,

j'aimerais stocker une matrice carré dense de taille très importante et pouvoir avoir accès aux valeurs de la matrice avec la meilleure performance possible. Actuellement j'utilise un stockage par représentation physique classique avec l'operator() ( M(i,j) ). Sauf que l'accès aux valeurs de la matrice est vraiment très lent .
Quelqu'un aurait une idée?
Merci.

Ste.

[koala01]

Salut, et bienvenue sur le forum.

Je ne doute absolument pas que beaucoup de gens ont de très nombreuses idées, dont certaines n'ont sans doute pas leur place sur ce forum

Seulement, pour ton problème particulier, je crains fort que personne ne puisse en avoir une réellement efficace...

Le fait est que, quand tu viens demander de l'aide, il est important d'aider le plus possible les gens à t'aider...

En effet, plus tu pourra nous expliquer clairement le problème auquel tu sera confronté ou, si tu arrive à nous donner un code minimal compilable qui permet de reproduire le problème auquel tu est confronté, plus tu auras une chance de recevoir une aide efficace et adaptée à ton problème...

Avec un peu de code pour nous permettre de nous rendre compte, je ne doute absolument pas que quelqu'un sera en mesure de t'aider

[stchouan]

Salut et merci pour ta réponse.

Je comprends ta remarque et tu as tout à fait raison.

En fait je travaille dans une librairie C++ de calcul scientifique ( LibMesh si jamais tu connais..) dans laquelle plusieurs classes pour la gestion des matrices sont implémentées. Je ne peux donc pas vraiment donner de bout de code compilable et en fait je n'ai ni de problème de compilation ni d'exécution mais de performance.
Voici ce qui se passe en gros.
Pour des raisons de calcul, j'ai besoin d'avoir accès à un certain nombre de valeurs discrètes. Pour cela je stocke ces valeurs dans une matrice pleine de taille (n,n) en faisant appel à une classe "DenseMatrix". Cette classe est basée sur une représentation physique matricielle avec accès aux éléments via l'operator(). Le problème est que plus n est grand et plus ma résolution est ralentie en raison de l'accès aux valeurs dans la matrice. Je pense que c'est tout à fait normal mais je me demande juste s'il existe un moyen technique pour optimiser cette performance.
Je sais que ma question est assez générale mais je ne sais pas trop comment donner plus de détails.

Je suis preneur d'idées diverses et variées à ce sujet.

Merci.

Ste

[Goten]

n est grand comment? Sinon t'a essayé d'autres bibliothèques spécialisé dans les matrices? blitz++ pour ne cité qu'elle...

[stchouan]

n vaut 100000.
Avec n=100, mon temps de résolution n'est pas altérée donc ca va. Mais quand je passe à n=1000, ca commence à vraiment ramer. Pour donner un ordre d'idée, avec n=100 j'ai un temps CPU de 15 secondes. Avec n=1000, je passe à 2500 secondes.
Je vais me renseigner sur blitz++..