Bonjour,
j'ai programmé la formule de BBP. Seulement j'aimerais programmer la formule de Fabrice Bellard car elle est presque 2 fois plus rapide (43%).
Le problème vient du (-1)^i. Comment on fait? J'ai essayé de le faire normalement mais ça ne marche pas.
Pole.
Merci d'essayer de résoudre mon problème
(-1)^i=exp(i*ln(exp(iPi))=exp(-Pi)!!
Salut,
ça c'est faux : on a pas forcémentEnvoyé par indian58
, on pourrait par exemple avoir
Ceic étant dans la formule de Fabrice Bellard i est simplement un indice de sommation (donc un entier) !
Cordialement.
J'adoore le quiproquo
Bonsoir,
On a plus simplement pour n entier
En programmation, ce n'est pas compliqué, il faut utiliser une variable initialisée à K=1 et à chaque étape poser
-K -> K
Le seul problème est que ça ne marche pas.Bonsoir,
On a plus simplement pour n entier
En programmation, ce n'est pas compliqué, il faut utiliser une variable initialisée à K=1 et à chaque étape poser
-K -> K
Fichier Maple 6 : ici
(il contient la méthode et 2 formules pour Pi : la première est celle de BBP et l'autre est une formule qui est plus lente mais c'est pour montrer que c'est le (-1)i qui pose problème)
J'ai essayé aussi de prendre une fois sur 2 le - modulo mais ça ne marche pas.
Fichier .c : ici
(J'ai fais ma conversion en base 16 car je ne savais pas que printf le faisais)
Pour une vitesse plus rapide, j'ai fais une procédure avec une exponentiation modulaire un peu plus rapide que l'autre. (Au lieu de diviser par 2, on divise par 4)
Fichier .exe : ici
(Sans SSE2, peut-être qu'il y a quelqu'un qui ne la pas...)
Autre fichier .exe : ici
(Avec SSE2, bien plus rapide)
Pole.
Personne ne peut répondre?
Pole
