Page 1 sur 1

tsimplify

Publié : jeu. janv. 07, 2016 6:29 pm
par frederic han
J'ai un pb avec tsimplify:

Code : Tout sélectionner

0>> tmp:=sqrt(1-cos(u)^2)
sqrt(1-cos(u)^2)
// Time 0
1>> tsimplify(tmp)
i*abs(sin(2*u/2))*exp((-i)*u)

Re: tsimplify

Publié : ven. janv. 08, 2016 11:44 am
par parisse
en fait c'est normal qui posait probleme avec un input complexe en mode reel. Ca devrait etre mieux maintenant (j'ai mis a jour le source et les paquets debian, mais je n'ai pas le temps de faire win et mac pour le moment).

Re: tsimplify

Publié : ven. janv. 08, 2016 12:42 pm
par frederic han
Cette fois c'est correct mais les erreurs d'arrondis dans les complexes sont toujours un peu surprenantes:

Code : Tout sélectionner

0>> tmp:=sqrt(1-cos(u)^2)
sqrt(1-cos(u)^2)
// Time 0
1>> tsimplify(tmp)
((-i)*exp(i*u)^2+i)/2*exp((-i)*u)
// Time 0.01
2>> tcollect(tmp)
abs(sin(2*u/2))
// Time 0.01
3>> regroup(tcollect(tmp))
abs(sin(u))
// Time 0.01
4>> u:=19*pi/180
19*pi/180
// Time 0
5>> approx(tsimplify(tmp))
0.325568154457+7.72715225139e-14*i
// Time 0
6>> approx(tcollect(tmp))
0.325568154457

Re: tsimplify

Publié : ven. janv. 08, 2016 1:24 pm
par parisse
c'est inevitable, puisque tsimplify passe en complexe : tsimplify n'est pas une commande pour debutants...