Voila, je cherche un algo de calcul simple et non récursif pour obtenir une approximation à une précision arbitraire des fonctions : sin,cos,tan,acos,asin,atan et leur version hyperbolique sachant que je dispose seulement des fonctions : a+b,a-b,a*b,1/a,1/sqrt(a) , c-a*b et a*b+c.
 
j'ai cherché le forum et google mais je trouve que les developpements limités et si je les utilisa à un ordre tres grand ben l'approx est merdique.
 
Merci d'avance
 
<<OOOOPS j'avais pas vu le post un + plus bas ...>>>
apparement c à coup de série entiére donc ? qqs précisions ?

Qu'est ce que t'entant par ordre très grand?
 
T'as essayé l'approx par polynome de lagrange?

 
Ben pour la précision dont j'ai besoin je dois allez à l'ordre 6/7 ou + ... et ca me prend trop de temps.
 
Lagrange c polynomial donc en x puissance y ... je peut pas me le permettre trop couteux.

Demande à Harko de te le coder en assembleur magique qui va plus vite que tout!
Plus sérieusement je vois pas trop...

 
 
Je crois que t'as de bonnes approx avec des poly d'orde inferrieur à 5 mais je suis pas sur.
 

C déjà ce qu'il fait

ah euh oui ca explique les restrictions dans l'énoncé...

Bon le DL 3 pour les fonctions hyperboliques semblent suffire ....
me retse les normales et les inverses

T'as essayé les série entières?
 
Ca sera sasn doute pas très différents des dl mais bon  

faut que je retrouve les formule mais si je m'abuse, y a des x^n qui traienent non ? Si oui, I'm doomed car les exponentiation me coutent la peau des baloches

Pas grave avec ton nouveau compilo qui optimise 2 fois plus...

 
 x^3 ca coute plus ou moin cher que de faire x*x*x, ou c ce que ça fait derriere?

D.L. ?

Développement limité

marde, j'avaispas ouvert mon oeil gauche, je proposais en fait, j'avais po vu

 
ha ok

 
x^3 mon compilateur ne connait pas..
 
C'est en quel langage ?
 
LeGreg
edit: je ne parle évidemment pas du exclusive or..

 
 
Ct un racourcis pour désigner la fonction qui mettrait à la puissance n  

 
oui mais si on ne connait pas cette fonction c'est un peu dur de repondre a ta question..
 
LeGreg

 
Ben si joel il la connait puisqu'il travaille avec  
 
EDIT: moi je peux pas inventer un nom de fonctino que je connais pas  

 
D'apres son post il ne connait que l'addition et la multiplication..
Ca m'etonnerait qu'il puisse faire un cube autrement et de maniere plus efficace qu'en faisant x*x*x.
 
LeGreg

 
Ca fait presque une semeine j'avais oublié moi  

c'est quoi le problème avec les D.L., je fais le boolay là, mais je comprends pas trop pourquoi t'es pas satisfait? et puis si je me souviens, les séries entières admettent des DL, donc, ça change rien.
 
sinon tu te tapes une tapes une table, à l'ancienne

Joel, tu as droit au branchements ?
 
LeGreg

Bon je recentre :
je bosse en altivec donc pas de branchement, tout inliné.
j'ai à ma disposition un cretain nombre de fonctions ultra-rapides précablés :
 
a+b
a*b+c
a-b
1/a
1/racine(a)
 
maintenant je dois évalué sin(x) (x quelconque) avec un mix de ces fonctions
sans sauts ni boucles.
 
J'ai implementé le sin avec les séries entiéres mais le pb ce que les x*x*x c trop lent
 
 
EDIT :
 
j'ai trouvé ca :
 
sin(x) = (-0.352173719 + 1.845955302*x - 0.005761952*x*x)/100
pour x <= 45
 
pour x > 45 degrees:
sin(x)  = sqrt(10000-100*sin(90-x))
 
ca a l'air correct ...

taylor expansion ?
 
a mince si x est quelconque ben pouf (dommage que tu ne puisse pas le restreindre

 
ben quand je disais branchement je pensais a un truc comme ca: tirer parti des symetries du probleme..
 
LeGreg

Put1 merde g pas fait gaffe
ca me force a faire un if a chaque fois
 
Bon back to the drawing board ...

c'est si lent que ça les fonctions de bases?

les fonctions de bases sont pas lentes. Seules la multiplication répété de la même donnée pose probleme. Mon but est d'effectuer cette approche avec un minimum d'opérations (en evitant les exponentiation sauvage) et d emanière linéaire (cad sans if ou while)

Juste une petite remarque sur le calcul des polygones :
 
((((-0.005761952*x) + 1.845955302) * x) - 0.352173719) /100
 
ca donne la même chose que ce qui suit  
 
(-0.352173719 + 1.845955302*x - 0.005761952*x*x)/100  
 
et ça enlève une multiplication.

OUI !!
j'y avais pas pensez !
Merci
 
/me va tester

Me stupid mais c'est quoi un altivec
 
on dirait le nom d'un de ces dinosaures de l'info
(quand tous ces macro ordinateurs s'appelaient eniac, univac )
 
Quoique le jeu d'instructions me fait penser
aux combiner de ma Geforce..
 
Alors c'est une machine differentielle on chip ?
 
LeGreg

SIMD pour G4

AltiVec est une extension SIMD pour les processuers G4,G5 de Apple/Motoroal/IBM. du MMx en beaucoup mieux quoi

Bon si je déroule les DL de maniére à ls transformées en a+b*c ca a l'aire de marcher ... Je vais tester le tout je verrais bien ...

le MMX c'est de la blague à côté... ça n'en a que le nom de SIMD