jordan
Modérateur : xcasadmin
-
- Messages : 1137
- Inscription : dim. mai 20, 2007 7:09 am
- Localisation : Paris
- Contact :
jordan
Salut,
j'ai des soucis de stabilite avec jordan selon le contexte.
s'il ne peut pas repondre, alors ca plante.
Ex:
A:=companion((x^2-2)^2,x)
il me semble jordan(A) plante si l'on n'active pas les racines carrees, mais fonctionne avec.
En revanche:
A:=companion((x^3-2)^2,x)
alors jordan(A) plante de toute facon, meme avec racines et complexes.
%%%%%%%%%%%%%%%
a part ca, tant que j'y pense:
les matrices avec des fractions ne passent pas bien vers pari:
Ex:
pari_charpoly(A) marche avec
A:=matrix([[2]])
mais pas avec A:=matrix([[2/3]])
a+
Frederic
j'ai des soucis de stabilite avec jordan selon le contexte.
s'il ne peut pas repondre, alors ca plante.
Ex:
A:=companion((x^2-2)^2,x)
il me semble jordan(A) plante si l'on n'active pas les racines carrees, mais fonctionne avec.
En revanche:
A:=companion((x^3-2)^2,x)
alors jordan(A) plante de toute facon, meme avec racines et complexes.
%%%%%%%%%%%%%%%
a part ca, tant que j'y pense:
les matrices avec des fractions ne passent pas bien vers pari:
Ex:
pari_charpoly(A) marche avec
A:=matrix([[2]])
mais pas avec A:=matrix([[2/3]])
a+
Frederic
-
- Messages : 1137
- Inscription : dim. mai 20, 2007 7:09 am
- Localisation : Paris
- Contact :
-
- Messages : 1137
- Inscription : dim. mai 20, 2007 7:09 am
- Localisation : Paris
- Contact :
c'est a cause du solver, il renvoie 2^(1/3) au lieu de -2^(1/3), je corrige ca lundi, si tu es impatient de tester tu peux faire le patch suivant:
diff solve.cc solve.cc~
376c376
< newv.push_back(delta_x + (negatif?-tmp:tmp));
---
<newv> newv.push_back(delta_x + tmp);
diff solve.cc solve.cc~
376c376
< newv.push_back(delta_x + (negatif?-tmp:tmp));
---
<newv> newv.push_back(delta_x + tmp);
Dernière modification par parisse le sam. déc. 15, 2007 4:23 pm, modifié 1 fois.
-
- Messages : 1137
- Inscription : dim. mai 20, 2007 7:09 am
- Localisation : Paris
- Contact :
-
- Messages : 1137
- Inscription : dim. mai 20, 2007 7:09 am
- Localisation : Paris
- Contact :
-
- Messages : 1137
- Inscription : dim. mai 20, 2007 7:09 am
- Localisation : Paris
- Contact :