COS ET SIN PLUS RAPIDE QUE MATH.H
coucou747
Messages postés
12303
Date d'inscription
mardi 10 février 2004
Statut
Membre
Dernière intervention
30 juillet 2012
-
24 juil. 2004 à 13:44
BoboLaricot Messages postés 4 Date d'inscription vendredi 6 juin 2008 Statut Membre Dernière intervention 25 novembre 2008 - 16 oct. 2008 à 12:20
BoboLaricot Messages postés 4 Date d'inscription vendredi 6 juin 2008 Statut Membre Dernière intervention 25 novembre 2008 - 16 oct. 2008 à 12:20
Cette discussion concerne un article du site. Pour la consulter dans son contexte d'origine, cliquez sur le lien ci-dessous.
https://codes-sources.commentcamarche.net/source/24821-cos-et-sin-plus-rapide-que-math-h
https://codes-sources.commentcamarche.net/source/24821-cos-et-sin-plus-rapide-que-math-h
16 oct. 2008 à 12:20
16 oct. 2008 à 12:20
Pourquoi stocker une table des sinus quand on sait que sin(x)=cos(x-PI/2) ?
A ta place je mettrai ça dans la fonction sinus :
float sinus(float)
{
return cosinus[x+PI/2]
}
Ca fait quand même une grosse économie de mémoire, non ?
14 avril 2008 à 20:29
1 oct. 2005 à 21:03
Pour DevCpp 4.9.9.2, avec mon Duron 1 GHz :
Editeur de liens : retirer les symboles de l'exe : [YES]
Génération du code : générer instructions : [i586]
Génération du code : utiliser les fonctions spécifiques : [MMX]
Optimisations du code : Toutes
En moyenne, les fonctions classiques sont 3 fois plus efficaces ! (170ms au lieu de 510ms)
Sans changer les options, les fonctions classiques (1500ms) sont 2 fois moins rapide que les fonctions modifiées (750ms).
14 déc. 2004 à 01:58
pour la fonction cosinus son devellopement limité en 0 est :
1-(x²/2)+((x²)²/24).......
au plus tu ajoutera de therme au plus tu sera preci!voila jespere tavoir aidé!
1 oct. 2004 à 22:31
pas assez de mémoire pour arriver a tout calculer...
expérimentation sur mon site
1 oct. 2004 à 21:48
Avec 3 000 000 de valeurs, j'ai un facteur de x3 et avec 30 000 000, j'arrive à un facteur 9.
7 août 2004 à 18:54
29 juil. 2004 à 12:56
int long unsigned n[1000]
for (x=0 x<1000;x++){
for (y=0 y<1000;y++){
if (x*x+y*y==1000000){
n[x]=y;
}
int diametre=3141592;
on a obtenu un quart de cercle de rayon 1000 (1000 pour être précis)
ensuite, on cherche a combien doit être la distance entre deux points du cercle séparés d'un degré...
diametre=(diametre-diametre%360)/360;
ensuite, caclul basique :
int cos[90];
int sin[90];
x=0;
y=0;
j=0
for (i=0;i<1000;i++){
if ( (n[i]-y)*(n[i]-y)+(i-x)*(i-x)==diametre*diametre)
x=i;
y=n[i];
sin[j]=y;
cos[i]=x;
}
}
29 juil. 2004 à 10:00
J'essayerais d'en faire une ebauche qd j'aurais un moment ...
djl -> oui, puis ca me paraissait bizarre que ce soit le role du préprocesseur d'évaluer des expressions en paramètres templates ...
et oui, je suis d'accord avec toi, du point de vue conceptuel on pourrait voir ca comme un préprocesseur c++
28 juil. 2004 à 22:11
28 juil. 2004 à 22:07
28 juil. 2004 à 21:55
if a=sqrt(b)
c'est pareil que
if (a*a=b)
28 juil. 2004 à 21:46
28 juil. 2004 à 21:44
ce serais plus simple
Moi, je vais commencer par expérimenter ce soir sur tibasic, puis demain sur pc, si j'y arrives...
28 juil. 2004 à 21:39
28 juil. 2004 à 21:34
ça te permet de faire vraiment ta fonction...
28 juil. 2004 à 21:32
dans la logique ca pourait etre fais par un preprocesseur mais seul un compilateur permet de faire des controle statique (verification de code)
je suis pas sur que parler dans ce cas de preprocesseur c++ soit totalement faux mais je retir quand meme ce que j'ai dit ;)
coucou747 >merci, je vais voir avec ta solution, mais au niveau perf je sais pas si ca sera plus rapide qu'un tableau
28 juil. 2004 à 15:25
c t a djl que je parlais :p
28 juil. 2004 à 14:27
moi, je ne sais rien du préprocesseur, mai en math, je me débrouille
28 juil. 2004 à 11:06
as-tu un lien pour montrer ta 'certitude' ?
27 juil. 2004 à 23:10
i=0;
for (x=-rayon,x<rayon, x++){
for (y=-rayon,y<rayon, x++){
/*formule permetant de calculer la longueur séparant l'origine de (x;y)*/
if (sqrt( y*y + x*x)=rayon){
pointx[i]=x;
point[i]=y;
i++;
}}}
ici, t'as tout les points du cercle, selon le rayon demandé, tu remplace le rayon par 1, et t'as cosinus en fonction du sinus, tu cherches un peu tu peux trouver une fonction permetant de faire ceci...
27 juil. 2004 à 18:50
pour le
"oui, ca va faire bosser le preprocesseur c++ ;)
=> oulala :o"
désolé mais c'est le preprocesseur qui travaille, c'est une certitude, le compilo ne connais pas les templates et le code doit etre genere avant la compilation
les template sont un mechanisme de generation de code, quel rapport avec la compilation ?
a ton avis pourquoi on appel ca meta programmation ?
27 juil. 2004 à 15:44
=> oulala :o
non ca fait bosser le compilo, pas le préprocesseur ...
l'exemple de djl se contente de déplier la boucle, pas le calcul du cosinus ...
il faudrait faire une opération comme :
template
int f()
{
return i;
}
int my_cos()
{
return f</* calcul du cosinus */>();
}
bon c pas terrible, c la 1ere solution qui me passe par la tete ...
ensuite, coupler cette solution a celle de djl pour obtenir un tableau de valeur entierement résolu a la compilation
27 juil. 2004 à 00:07
26 juil. 2004 à 22:55
double v_cos[100];
template void init()
{
v_cos[N] = cos((double)N);
init<N-1>();
}
template<> void init<0>()
{
v_cos[0] = cos(0.0);
}
et l'appel de init<100>() genera le code completement deroulé
je te conseil le addison wesley sur les template, ya que ca
26 juil. 2004 à 22:37
26 juil. 2004 à 19:42
et comme je l'ai dit, attention à la taille de l'exe
26 juil. 2004 à 14:41
24 juil. 2004 à 20:15
J'avais totallement oublié fmod.
Mais bon, corrigez moi si j'ai tord, je pense pas que fmod soit plus rapide que ce que je présente : je ne crois pas qu'on puisse faire plus cours... (sauf peut etre avec de l'asm, hein brunews)
Donc en prennant ma méthode tu economises quelques cycles (genre ceux d'appel a fmod). :D
24 juil. 2004 à 17:50
sinon j'ai aussi bien preciser que c'etait 10 fois plus rapide dans l'intervalle et precalculer 20000 fois cos et sin c'est rien du moment qu'apres c'est rapide , ya que ca dans les jeux (dans certains, les 3/4 du temps de chargements c'est du precalcul)
si tu ve jeux te sort la version c++ qui fait l'initialisation à la compilation (par le preprocesseur) mais gar à la taille du prog
Xs > je vais essayer de mettre en place de cette maniere, en fait ce qu'il faut c'est une fonction qui calcul le modulo d'un reel plus rapide que fmod (ou alors j'utilise fmod ?)
24 juil. 2004 à 16:38
24 juil. 2004 à 15:50
garslouche > c'est la méthode al plus rapide et la plus courante dans un jeu : on stock et on prend après. Les 20000 fois au départ tu t'en fiches : elles se font dans l'initialisation de ton jeu (et si c'est en thread, c'est encore plus invisible).
djl > j'ai lu ton code. Si je fais tr_cos(95851) je crois que je vais devoir patienter un bout de temps avant d'avoir l'abscisse curviligne de 95851 dans [-Pi;Pi].
Donc si tu veux calculer le sinus d'un réel "b" tel que b > 0 (sinon y'a une petite modif a faire) tu dois résoudre
b - k*Pi > -Pi/2 & b - k*PI < Pi/2 (la 1er equation est utile : ca n'induit pas b < 0 )
Je trouve que b/Pi -1/2 < k < b/Pi+1/2.
Donc pour le sinus de b (b > 0), tu calculs b/Pi+1/2 que tu castes en "int" (pour obtenir aisement la partie entiere) puis tu déduis k*Pi de b . Et la zoup, en 1 division, une addition, une soustraction, une multiplication, tu obtiens l'indice du tableau v_sin.
Pour le cos, tu résouds b-k*Pi < PI & b-k*pi > PI (peut etre qu'il y a une différence pour b < 0).
Bon, je peux me tromper, mais je suis sur qu'il y a NETTEMENT plus rapide que des incrémentations ( tu te rends compte ? au lieu d'un grand nombre d'inc tu fais juste v_sin[b-((int)(b/Pi+1/2)*Pi]
ca semble barbare mais c'est rapide)
cordialement
24 juil. 2004 à 15:15
Je ne vois pas en quoi c'est plus rapide que math.h puisque tu commences par appeler 20000 fois la fonction cos de math.h !
Si tu veux savoir comment un ordi calcule un cos tu peux regarder ma source "math.h reprogrammé"
24 juil. 2004 à 13:44
tu pourrais détailler un peu ? expliquer pourquois ça marche ? En général, je suis plutot pas mauvais en math, mais dès que je vais sur ce site, je déprime... lol