Discussion :

[wurst]

bonjour,

j'ai l'ambition de créer une bibliothèque en c++ pour faire du calcul matriciel. Dans l'idéal cette bibliothèque devra s'utiliser avec une syntaxe proche de celle de matlab.

Cette bibliothèque utilise les template.

Dans une première version, je me suis inspiré du conteneur vector pour la possibilité d'écrire les surcharges d'opérateur indépendamment de la dimension de la matrice créé.

Ainsi les vecteur<float> utilisent les mêmes opérateurs (en théorie) que les vecteurs<vecteurs<float>> ... peu importe la taille.

Un gros défaut de ma méthode actuelle : la mémoire.

En effet j'ai fait l'expérience de créer la matrice de transformé de fourrier. Ma matrice crée est de la taille 8192 au delà, le programme plante. cette matrice occupe 1Go de mémoire et met beaucoup plus de temps à être calculée que sous matlab.

ma question: quel est selon vous la meilleur façon de bien partir pour créer la bibliothèque?
quelle classe? s'appuyer sur des classes standard existante?
bref je vous propose d'y réfléchir avec moi.

[Joel F]

Spa trivial le calcul matriciel en C++.

Deja vector<vector<T>> est tres tres mauvais comme conteneur de base car il alloue des bloc discontinue et donc empeche le cache et les pages d'etre accédé proprmeent. COmme d'hab, c'est dans les vieux pots etc... donc regarde du coté de l'allcoation façon NRC.

Ensuite, expression tempaltes ou alors tu va crouler sous le point des temporaires.

[wurst]

je cherche pas quelque chose de trivial .

Le coté mathématique de l'affaire ne me fait absolument pas peur. C'est plus la partie c++ j'ai découvert pas mal de chose coté c++ depuis que j'ai commencé. Et je réalise que je suis pas suffisamment au point pour faire quelque chose de propre seul.


NRC => ok petite recherche rapide dessus et je tombe sur pleins de doc. Avant de tout lire, ce qui me prendra pas mal de temps, est-ce que c'est allouer un grand tableau et simuler une matrice2D en calculant les offset des différentes lignes?

"Ensuite, expression tempaltes ou alors tu va crouler sous le point des temporaires." => j'ai découvert le principe des Template que récement, je suis pas sur d'en connaitre assez pour comprendre ce que tu veux me dire.


lien en rapport NRC (à lire demain )

http://www.nrbook.com/a/bookcpdf.php

[Joel F]

quand je dis pas trivial c'ets au niveau implantation pas math. Pour bosser depuis 5 ans sur un projet simialire ... je t'assure que c'ets pas le truc qu'on design sur un coin de nappe le dimanche au pique-nique.

Les expressions tempaltes sont une technique qui permet de travailler à un grain plus gros que l'operateur et de proposer des optimizations global.

NRC => ok petite recherche rapide dessus et je tombe sur pleins de doc. Avant de tout lire, ce qui me prendra pas mal de temps, est-ce que c'est allouer un grand tableau et simuler une matrice2D en calculant les offset des différentes lignes?

Oui, et ca a l'avantage de s'écrire récursivement

Dasn ton lien, le chapitre débutant page 15 doit contenir ce que tu cherches

[Flob90]

Pour les template, Il y a un cours d'introduction à la méta-prog sur ce site où il est question à un moment de matrice (écrit par Laurent il me semble).

L'idée est d'utilisé la méta-prog pour que le compilateur génère le plus de calcul possible, ca allègera le traville de l'ordinateur au moment d la compilation (en gros il n'aura plus que des addition et multiplication à faire pour calcul chaque élément de la matrice et pas à faire les calcul un à un.

[omc24]

Bonjour,
La majorité des lib d'algèbre linéaire se basent sur blas (level 1,2 et 3) et lapack. Ce sont des librairies plutôt bas niveau écrites en fortran. Ces librairies sont des dinosaures du calcul matriciel, elles sont terriblement rapides et débuggées.
Ne pas les utiliser serait - à mon avis - une faute de goût . Il y a des passerelles avec le C++ via les numeric bindings et boost uBlas.
Bon courage,
Damien.