FONCTIONS USUELLES (TRIGO) EN METAPROGRAMMATION

Commentaire de BruNews le 07/06/2007 15:50:03 administrateur CS

récursion = risque explosion de la pile + lenteur.

Commentaire de coucou747 le 07/06/2007 18:44:53 administrateur CS

gcc ne fait pas tout ce qu'il pourrait faire a la compilation, je testerais sous vs ce soir... (mais ma factorielle reste inferieure en temps a une boucle)

Commentaire de coucou747 le 07/06/2007 21:07:07 administrateur CS

... j'etais sous windows a l'instant pour voir ce que ca rendait :
Suite a quelques conseils de NitRic, j'ai pris cette fonction main
int main(){
int a;clock_t begin, end;begin = clock();
for (int i=0;i<9000000;i++){
a=1;
for (int j=1;j<=15;j++)
a*=j;
}
end = clock();
std::cout << (unsigned long long)(begin - end) << " - " << a << std::endl;
begin = clock();
for (int i=0;i<9000000;i++)
a=Math<int, 15>::factorielle();
end = clock();
std::cout << (unsigned long long)(begin - end) << " - " << a << std::endl;
return 0;
}

sous visual en release:
18446744073709551148 - 2004310016
0 - 2004310016

sous linux (g++), ou sous visual en debug, j'ai quelquechose de different :
18446744073708611616 - 2004310016
18446744073708691616 - 2004310016

il te faudrait l'asm fourni par quel compilateur ? (pardonne moi, je ne sais pas utiliser vs, je ne sais pas lui faire cracher un assembleur, ni meme a g++ d'ailleur, et le man page de g++ fait 8000 ligne :()

Commentaire de coucou747 le 08/06/2007 05:32:08 administrateur CS

t'as les asm dans les fichiers html, ils correspondent aux mains commentes dans tcml.cpp

Commentaire de coucou747 le 09/06/2007 01:55:35 administrateur CS

sous g++ faut compiler avec l'option -O2 et on a les memes optimisations

Commentaire de coucou747 le 09/06/2007 12:08:06 administrateur CS

en debug, ca ne donne rien de bon a ce sujet sous vc apparement...

la syntaxe, tu reproches le class T ? tu prefererais un truc genre Math<50>::... directement ?

Commentaire de coucou747 le 09/06/2007 16:55:09 administrateur CS

personellement, je ne suis pas assez maso pour ecrire moi meme un DL50 de cos(x) en 0... on a des machines pour ca... autant faire une lib qui nous epargne ce travail (chose faite ici)

Commentaire de BruNews le 09/06/2007 23:15:11 administrateur CS

Bien d'accord coucou, faisons bosser les machines.
Mon idée est juste de mettre ce code dans un outil qui donne la constante à insérer dans le code, ça évite de se trimballer ces templates partout.

Commentaire de BruNews le 10/06/2007 01:34:28 administrateur CS

"version template bat l'asm" n'a pas de sens, le compilo produit un asm... donc au moins imitable.
Le 'metaprogrammé' donne les mêmes résultats précalculés avec des valeurs non connues à la compil ?
Bien des lustres que je n'ai pas abordé ce sujet, faudrait voir.

Commentaire de coucou747 le 10/06/2007 08:01:49 administrateur CS

juste dans le but de virer les templates ? tu lacerais mapple ou autre ? pour avoir un dl50...

apres t'aurais un code incomprehensible (enfin ca, on a l'habitude...), et impossible a modifier

Commentaire de coucou747 le 10/06/2007 11:16:34 administrateur CS

math.h te renvoi souvent trop de decimales pour cos et sin... tu t'en fous generalement que ton float ai une valeur legerement erronee car finalement, quand tu convertis en int a la fin pour placer ton point a l'ecran, bah tu perds cette precision, et meme si t'avais un decalage, decaler d'un pixel ou deux, c'est pas forcement genant...

Calculer un DL8, ca donne Math<float, 4>::cos(...) de cosX, et la, X peut etre variable, et on a une meilleur vitesse que cos dans <math.h>

18446744073709441616 - 1
18446744073708791616 - 1


pour float a;
clock_t begin, end;
float b=3;
begin = clock();
for (int i=1;i<900000;i++){
a=cos(b/i);
}
end = clock();
std::cout << (unsigned long long)(begin - end) << " - " << a << std::endl;
begin = clock();
for (int i=1;i<900000;i++){
a=Math<float, 6>::cos(b/i);
}
end = clock();
std::cout << (unsigned long long)(begin - end) << " - " << a << std::endl;

(en compilant avec -O2 -O3 sous g++, on est plus rapide que math.h) et un DL8 sur des fonctions genre cos(x), c'est deja d'une bonne precision...
bon, et sans parler de ca, pour les constantes, c'est du bonheur :)
De plus, calculer sur toute classe qui overload les operations + / * et -, je trouve ca plutot utile personellement...