Page 1 sur 1

solve et flottants

Publié : lun. oct. 15, 2007 9:51 pm
par frederic han
Bonjour a tous,

J'ai des erreurs de segmentation avec solve sur l'exemple suivant:
(En Digits:=15, ca n'a pas l'air de le faire en Digits inferieurs)
(je suis en mode maple)
0>> P:=t^2*u^2/2+t*u-1;
t^2*u^2/2+t*u-1

// Time 0.01
1>> Q:=normal(2*P);
t^2*u^2+2*t*u-2

// Time 0.01
2>> nSolve(subs(t=0.1,P),u);
7.3205080756887729352744632

// Time 0
3>> solve(subs(t=0.1,Q),u);
7.3205080756888,-27.320508075689

// Time 0.01
4>> solve(subs(t=0.1,P),u);
Erreur de segmentation


A part ca factor est plus joli, il ne me reste jusqu'a present que
factor(x^4+1,exp(I*pi/4) qui ne se scinde pas autant qu'il le devrait:

0>> factor(x^4+1,exp(I*pi/4))
(x^2+(-(sqrt(2)))*x+1)*(x^2+sqrt(2)*x+1)

// Time 0.03
1>> factor(x^4+1,exp(I*pi/4)+1)
(x+(-(sqrt(2)))/(1-I))*(x+(I)*sqrt(2)/(1-I))*(x+(-I)*sqrt(2)/(1-I))*(x+sqrt(2)/(1-I))


Merci pour tout

Frederic

Publié : jeu. oct. 25, 2007 7:09 am
par frederic han
Dans le meme genre, en 16 chiffres j'ai des Pb avec:
A:=[[1,2],[3,4]]
B:=approx(A)
inv(B) plante

Vous aussi?

sinon en approx_mode, les tableaux sont difficilement utilisables, C[i,j] ne marche plus, je dois mettre des exact partout

A+
Frederic

Publié : jeu. oct. 25, 2007 8:52 am
par parisse
j'ai corrige le factor (esperons que c'est la bonne ce coup-ci!), le inv en approx avec + de 15 chiffres, par contre je n'ai pas reussi a reproduire le plantage du solve. J'ai aussi ajoute qu'en mode approx les indices marchent sur les matrices.
tout ca devrait etre dispo en version xcas_root/user dans environ 1/2 heure.

Publié : jeu. oct. 25, 2007 11:26 am
par frederic han
Oui, ca marche, super, le plantage du solve c'etait aussi avec plus de 15 chiffres, mais dans cette nouvelle version ca ne le fait plus, dans celle d'hier si.

Merci

Publié : jeu. oct. 25, 2007 12:36 pm
par parisse
ca reste mysterieux quand meme puisque je n'ai rien corrige pour le solve, ca ressortira peut-etre un jour, en attendant, tant que ca marche, je n'y touche pas!