Calcul rapide d'une gaussienne.

**Nebulix** · 08/10/2010, 11h29

Envoyé par ToTo13

Je cherche à avoir le code le plus rapide possible !
Un pré-calcul permet un accès direct à la valeur et évite des calculs forcément plus coûteux.

On peut aller plus vite en diminuant la précision. Pour des calculs longs et complexes celà peut être dangereux ! Il parait que les meilleurs compromis sont les approximants de Padé (et que le FPU les utilise)
Pour un gain de 20%, je ne le ferais pas.

**pseudocode** · 08/10/2010, 11h35

Envoyé par Nebulix

Voici le bout de code qui m'a permis d'évaluer la vitesse : 20 sec environ pour 200 millions d'appels
Je suis vraiment très curieux de savoir combien tu mesures sur ton système

Equivalent Java de ton bench:

Code java :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
int LOOP = (int)10E6;
long t0 = System.currentTimeMillis();
 
for(int j=0;j<LOOP;j++)
	for(int i=0;i<20;i++) exp(-i*1.33);
 
long t1 = System.currentTimeMillis();
 
System.out.println("duration : " + (t1-t0));

Sur mon PC du boulot (Java 6-u21 32bits), en utilisant l'appel natif Math.exp() j'arrive à 16.212 secondes. En utilisant l'interpolation de la table j'arrive à 2.89 secondes. Soit un gain de facteur 5.

Edit: et avec le hack sur les flottants, on tombe à 0.156 seconde

**ToTo13** · 08/10/2010, 11h38

Envoyé par Nebulix

Pour un gain de 20%, je ne le ferais pas.

C'est un gain sur l'ensemble global de mon projet et la perte de précision n'engendre par d'erreur dans le reste de mon programme.
Donc c'est au final une très bonne solution.

**Aleph69** · 08/10/2010, 20h40

Bonsoir,

au-delà de l'optimisation de ta fonction elle-même, il faut aussi que tu utilises ta fonction de manière optimale, en particulier dans les boucles.
Il ne faut pas non plus oublier de compiler ton code avec un tag d'optimisation.

**ToTo13** · 10/10/2010, 16h39

Envoyé par Aleph69

Il ne faut pas non plus oublier de compiler ton code avec un tag d'optimisation.

J'avoue ne pas les connaître en java :-(
En fait, je n'ai jamais réellement eut à les utiliser.

**ToTo13** · 10/10/2010, 16h40

Je marque cette discussion "Résolu" car l'optimisation proposée par PseudoCode est pertinente, mais je suis ouvert à toutes nouvelles propositions !!!