Discussion :

[f56bre]

Salut,

Je suis amené à appeler souvent les fonctions sinus et cosinus de la librairie math.h et j'aurais aimé savoir comment elles calculent? Existe t'il des algo qui permettent de gagner de la vitesse quitte à perde (un peu) en précision?

merci d'avance

[Jean-Marc.Bourguet]

Salut,

Je suis amené à appeler souvent les fonctions sinus et cosinus de la librairie math.h et j'aurais aimé savoir comment elles calculent?

Certains processeurs ont des instructions pour. Pour les autres, je connais deux familles d'algorithmes:
- CORDIC
- approximation polynomiale par morceau

Existe t'il des algo qui permettent de gagner de la vitesse quitte à perde (un peu) en précision?

Faire une table et se limiter au contenu de cette table ou interpoler linérairement peut gagner ou pas. Suivant la taille de la table, le processeur et une quantité d'autres facteurs. Le seul moyen d'être sûr est de mesurer.

[ToTo13]

Bonjour,

il y a les développements et certaines fonctions d'approximation (pour la fonction sqrt).

Pour gagner du temps : si tu appeles souvent les mêmes valeurs ou des valeurs entières en degrès, remplit un tableau avec ces valeurs en début de soft et pioche dedans (c'est le cas pour un de mes soft...).

[TNT89]

Bonjour,

Je m'excuse de remonter le post mais j'aurais besoin de savoir si les processeurs intègre des solutions en "dur" pour le calcul de exp, cos, sin...
Je viens d'essayer de produire un code à peu près optimisé en C++ à partir de la métaprogrammation et il ne peut tenir la comparaison avec la bibliothèque (je me fait battre d'un bon coefficient 2)...

Le code que j'applique :

Code :

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

1
2
3
 
exp(x) = Sum(x^n/n!,{n,0,Inf})
       ~  ( ( ( (x/n + 1) * x/(n-1) + 1) * x/(n-2) + 1) * x/(n-3) + 1) * ...

Je passe à côté de quelque chose?

[Mac LAK]

Je m'excuse de remonter le post mais j'aurais besoin de savoir si les processeurs intègre des solutions en "dur" pour le calcul de exp, cos, sin...

Oui, cela s'appelle le coprocesseur mathématique, coprocesseur arithmétique ou encore FPU (Floating-point unit). De plus, sur certains types de calculs, ton compilateur peut même générer du code SIMD qui vont accélérer les calculs par un facteur significatif (tu peux avoir un facteur 8 par rapport à un calcul "simple" utilisant le FPU sur certaines opérations).

Les processeurs de type Pentium actuels intègrent aussi bien un FPU que des registres SIMD (appelés "SSE" chez Intel)... Vouloir battre des fonctions câblées avec du logiciel "pur" est franchement illusoire, je préfère te prévenir.

La seule solution que je connais pouvant battre les circuits dédiés est le précalcul de tables : c'est monstrueusement efficace sur les fonctions trigonométriques, par exemple, car elles sont comme tu le sais périodiques.
C'est nettement moins simple avec des fonctions non-périodiques comme la racine carrée ou les logarithmes / exponentielles, mais en fonction de ton besoin et du domaine d'application, cela reste possible.

[TNT89]

Merci de ta réponse!

Pour les FPU et le multi-core : comment ça se passe? Chaque cœur a son unité? et pour les threads? (pour les Intel I7 par exemple...).