réflexion sur le calcul matriciel

Citation:

---

Envoyé par regis.portalez

Et quand je dis que les expressions templates sont à avoir à la fin, c'est parce qu'avant de les faire, il vaut mieux savoir faire le reste... Surtout si c'est pour rédiger une bibliothèque utilisée en local.

---

Le problème, c'est que les expression template ne sont pas factorisable dans un code déjà écrit s'ils n'a pas été pensé pour. Si on veut s'en servir, on est obligé de s'en servir dès le début.

Citation:

---

Envoyé par regis.portalez

Les expressions templates sont là pour corriger les défauts d'utilisation de la lib par ses utilisateurs.

---

Non. Ils sont là pour corriger des problèmes d'instanciation de classes et de combinaison des opérateurs associés.

Citation:

---

Envoyé par regis.portalez

Par exemple, il faut être naïf pour écrire :

A = (B*C) + D - 3*D.norm()*E...

---

Naïf ? Je dirais qu'il faut être normal pour écrire un tel code. Le principe d'une lib, c'est justement de permettre d'écrire le code de manière claire et lisible.

Citation:

---

Envoyé par regis.portalez

Quand on écrit une lib pour soi, on sait qu'on ne va pas écrire ça.

---

Ben si, justement. C'est le but même :) Sinon, autant ne pas écrire de lib :)

Citation:

---

Envoyé par regis.portalez

On va faire soi même la boucle en dur, comme un grand

---

Pourquoi pas ? Ca n'est qu'une duplication de code comme une autre :)

Citation:

---

Envoyé par regis.portalez

Quand on la livre à qqn, il faut prévoir le cas où le mec va pas s'embêter à réfléchir...

---

Le client a droit à une lib simple d'utilisation et efficace, quand moi je devrais m'ennuyer avec un truc odieux ? C'est pas très logique tout ça :)

Youhou les gens, le type du début a pas dit qu'il voulait distribuer sa lib au monde entier, il veut juste avoir un truc qui marche et qui lui fasse LU en à peu près autant de temps que MATLAB.

C'est pas la peine de le gaver avec des templates de fous.

Il voudrait faire un vrai lib distribuable (ou même commercialisable) je dis pas, mais là c'est pas vraiment le cas...

à voir le deuxième message du mec : "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."

Si vous voulez vous faire un débat sur la meilleur lib de matrice de niveau pro allez-y, mais ça ne répond pas au sujet.

Citation:

---

Le problème, c'est que les expression template ne sont pas factorisable dans un code déjà écrit s'ils n'a pas été pensé pour. Si on veut s'en servir, on est obligé de s'en servir dès le début.

---

Oui, mais il faut d'abord savoir comment on rédige des classes et des templates avant de le faire.
Sinon, c'est comme si tu disais à un type qui sait pas comme marche un moteur de fabriquer une voiture en commençant par la boîte automatique...

Bonsoir,

je pense au contraire que chacun cherche à répondre à la problématique initiale qui spécifie bien un problème de temps de calcul pour une matrice de Fourier comparativement à ce que fait Matlab. Pour approcher les performances de ce dernier, il faut s'inspirer des bibliothèques externes qu'il utilise (blas & lapack optimisés pour la réutilisation des mémoires caches + umfpack si ma mémoire est bonne). La moindre des choses c'est bien d'inviter les "candidats à faire aussi bien en C++" à utiliser les expression templates pour faire du plus haut niveau que fortran tout en maintenant des performances de calcul honnêtes. En ce qui concerne la factorisation LU proprement dite, imbriquer trois boucles ne suffira pas.

Citation:

---

Envoyé par regis.portalez

Youhou les gens, le type du début a pas dit qu'il voulait distribuer sa lib au monde entier, il veut juste avoir un truc qui marche et qui lui fasse LU en à peu près autant de temps que MATLAB.

<snip>

à voir le deuxième message du mec : "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."

<snip>

Oui, mais il faut d'abord savoir comment on rédige des classes et des templates avant de le faire.

---

Ouai sauf que c'est totalement contradictoire, si tu maitrises pas ce genre de chose alors l'objectif fixé est pas réalisable.

Petite parenthèse : les NRC ont une licence fort pénible. Je ne sais pas comment s'interprète le modèle d'allocation vis à vis de la licence, mais c'est un détail auquel il faut faire attention.
J'ai eu un exemple de lib opensource (OTB) dernièrement qui a dégagé ses dépendances aux NRC.

pour avoir eu le probleme icic, les algos numeriques sont sous licenses, l'allocateur memoire ets "Public Domain"

Merci pour l'info. C'est bon à savoir.