pari matdet
Modérateur : xcasadmin
-
- Messages : 1137
- Inscription : dim. mai 20, 2007 7:09 am
- Localisation : Paris
- Contact :
pari matdet
Salut, j'ai un probleme avec la 0.9.0 sous linux
(Je suis en 32 bits)
(ca le fait avec une version que j'ai compile, et avec le paquet debian aussi)
k:=50;A:=matrix(k,k,(i0,j0)->rand(21)-10):;
d1:=matdet(A,0)
d2:=det(A,bareiss);
d3:=det(A)
alors d1 est tres different des 2 autres.
en 0.8.6 les 3 coincident bien et c'est d1 qui semble faux.
Fred
Ex:
[[-3,10,-7,-2,0,-3,5,-8,-6,5,8,1,-9,9,-6,7,-4,3,-3,2,-9,-2,1,1,-7,-4,-3,-8,-5,-3,-8,1,-7,-1,-4,7,10,5,1,2,5,-6,8,8,-1,-4,-7,8,-5,-2],[-4,-7,2,-8,-2,4,-3,10,-4,0,-4,0,4,-6,3,-10,-3,-5,-3,2,7,1,-9,0,-1,-4,2,1,2,-6,7,-6,-1,3,7,9,-8,-7,0,-5,8,-3,-8,10,8,2,-9,-3,-4,-2],[2,10,6,3,9,-9,-10,-2,6,-2,-9,-2,-9,-4,1,-1,-1,1,4,-1,-9,-2,-6,9,4,-2,9,8,-6,0,3,-2,-7,-2,-9,-9,-2,10,0,-6,9,-1,4,-1,-1,3,8,1,3,-3],[-1,-3,10,-8,8,-7,4,-10,5,7,-2,-8,-3,-10,10,-10,7,-2,-5,3,-6,-3,-4,1,10,9,-8,1,-8,-1,-6,-8,-1,-2,1,5,-8,-5,-7,-1,1,6,-2,-5,10,-2,2,-5,8,7],[2,-9,-9,9,7,-8,-10,6,1,9,2,-9,-3,-8,-9,-3,6,8,-1,-8,6,9,-6,-1,-1,3,0,-4,9,-5,-4,9,-8,6,-1,-1,6,-8,-2,2,5,-9,-9,7,5,-6,-5,4,-9,8],[-9,-8,2,-9,-9,-8,-1,-4,-1,10,-1,0,8,10,8,5,-6,-9,-9,10,-7,-4,-5,9,5,6,4,4,0,-5,-4,-2,2,10,-7,-9,9,2,1,9,6,-10,-7,-8,5,-1,3,-1,5,-4],[-10,1,8,1,-6,6,-2,-5,7,-4,3,0,8,5,-8,-8,4,2,-9,-9,3,0,-9,5,1,2,2,3,1,7,9,-4,-1,-3,3,-6,-10,1,-5,-3,-2,-3,-6,-5,6,3,-4,-3,3,8],[-4,-3,10,5,-5,6,1,10,2,-8,-5,0,-8,2,5,-5,10,9,6,1,-8,-2,-7,5,-8,10,-4,-5,-5,-10,-6,0,10,5,-4,5,8,-3,3,8,0,-10,8,-4,-5,-4,5,4,10,4],[-10,4,-6,-1,-3,3,6,0,-1,-2,3,8,10,5,9,10,9,-8,8,-7,7,-4,-1,-10,-8,7,4,10,6,2,-6,10,5,-4,-4,8,3,10,-6,-5,2,-4,1,-3,-10,-7,-1,-4,4,10],[8,0,1,6,-7,-5,-7,-10,-5,4,5,-2,4,7,1,-5,4,4,-5,9,6,-8,3,4,3,-5,10,-9,-4,1,6,8,5,10,5,-1,7,-3,2,-2,-3,-9,9,-10,-9,5,-2,10,-10,-3],[10,-3,-1,0,-10,6,-10,8,8,-10,9,9,-10,-10,-8,6,10,5,-8,0,3,9,9,-7,-6,8,-4,6,8,-6,2,7,10,2,4,-4,3,-3,-4,-8,-3,-5,-10,-5,5,10,8,5,-7,6],[-5,-4,9,-8,-5,2,9,8,-9,-1,-10,-2,-6,1,9,-3,-2,9,9,3,4,-6,0,-1,10,6,4,8,2,9,6,-7,9,8,3,9,-4,-3,9,1,-2,6,-3,-2,-9,-1,-7,-10,-9,-6],[5,-3,5,-9,9,5,1,-5,-7,-4,-4,-7,6,7,-4,10,-2,7,-2,-7,1,4,7,-5,8,1,9,-10,-3,-3,1,2,0,5,2,-10,5,7,-8,10,9,-2,-2,2,4,-2,-9,3,-1,7],[-1,10,4,6,0,8,-7,5,-9,-6,5,3,-2,9,-7,0,8,8,1,9,-1,4,-5,6,-10,9,-1,8,-1,1,-4,0,6,-2,3,5,8,-5,-5,-4,-5,-1,10,5,-1,4,-8,10,6,10],[-3,3,8,4,10,1,2,-3,1,8,-9,7,-5,4,-3,2,-5,10,-9,-3,4,-8,4,6,-3,8,-5,5,-6,7,7,-8,8,-7,-6,5,-3,-5,-8,3,9,-2,-6,-6,9,10,-1,-10,9,0],[3,3,-1,2,-7,1,2,9,3,-2,10,-4,-5,3,-8,-6,7,-8,-7,-8,-9,9,-9,-10,1,4,0,-4,-10,8,-1,3,4,-5,-7,0,9,6,9,-7,-7,3,-6,-8,2,-1,5,9,-1,1],[4,5,1,10,9,6,-6,-10,1,-6,-1,4,-7,-3,5,3,5,-5,5,9,-10,-3,4,8,-2,2,-10,6,7,-10,-1,-3,6,-7,-10,10,-9,-5,9,-3,-4,3,2,1,10,-9,8,6,7,6],[7,9,-8,7,3,9,-2,-3,-6,1,10,2,3,-5,2,-10,-1,2,-3,-1,8,-5,-1,-8,7,-10,-6,4,3,-5,-2,-4,10,-5,10,-2,-10,-5,4,5,-6,-10,-5,-2,-4,2,-9,-1,-3,-9],[-4,9,6,9,-9,2,-10,-6,5,4,-7,-8,-10,6,-2,8,6,-4,-6,4,-9,-1,4,-9,-4,2,-4,-9,-3,-9,5,-7,8,6,6,-10,8,2,1,4,7,10,10,6,4,2,8,-6,-5,-1],[-5,4,-9,4,6,10,-6,7,8,10,-4,10,-7,10,2,0,3,3,10,-9,-3,1,-2,10,-8,-3,10,0,4,-6,3,-10,-1,-8,0,-5,-9,3,0,-5,10,-9,-6,6,1,-4,-3,-7,0,-1],[-8,0,5,-3,6,2,-7,-4,9,10,1,-2,4,9,8,-3,8,-3,5,-3,-2,3,6,3,7,3,-1,-4,4,-3,-9,2,-2,9,2,-6,6,-5,-6,5,-10,7,5,4,1,-2,8,-4,10,5],[-2,10,-2,-8,-6,-1,-1,-9,7,5,1,-9,0,5,0,6,-4,3,7,8,7,-10,-3,-1,10,0,10,-5,-1,-2,-3,2,-6,7,3,-4,0,6,7,-4,-5,5,-3,8,1,10,1,-8,6,7],[4,10,-10,-4,-1,3,1,8,-6,6,10,0,6,5,5,4,5,-1,0,-2,-3,8,4,-9,-6,8,-3,-10,-5,-3,-2,-5,6,-3,-2,-9,-4,7,-6,-7,10,7,-1,10,6,-3,-8,10,6,-6],[1,1,4,8,2,8,-1,1,2,9,6,1,-10,-9,8,0,-9,-4,0,8,-7,-1,0,3,-3,5,-6,-6,-10,-1,-4,5,-1,-10,-6,-10,-9,9,-8,9,1,-1,5,-8,10,7,-2,-3,0,2],[-10,-8,-5,-8,-1,-9,0,-5,0,1,1,0,0,5,0,-8,3,-3,1,-2,4,3,7,3,-4,-8,-5,2,-7,-5,8,8,-1,-2,3,-9,-8,-3,-1,-4,0,-10,-3,9,9,-3,0,0,-4,-1],[-9,-10,-8,5,3,-7,-8,10,-8,-10,0,-5,3,5,6,-8,-7,8,-7,-3,-9,-3,4,-10,1,-9,3,-4,7,1,10,0,-8,-10,8,5,4,-6,-2,5,1,-9,-5,5,10,-7,2,-7,-10,-1],[8,-4,7,-5,10,-7,6,5,-4,5,-5,-10,9,-8,4,-10,9,3,-7,-3,6,7,7,-1,7,1,-3,4,8,7,10,-6,9,-8,4,-1,5,4,-8,-4,10,-3,-8,-8,10,-3,-3,-4,5,-9],[6,7,-2,-5,9,-7,-10,-6,8,-3,-6,-10,-5,-7,1,10,-9,7,8,-1,-5,2,-5,-5,-6,3,-3,7,-10,1,-3,-9,4,-7,9,1,-7,7,-7,2,-8,2,2,-8,4,5,-1,4,9,5],[8,1,-4,8,-5,-6,8,9,8,5,-6,9,-7,-7,2,3,1,4,-2,0,-7,-7,-7,1,2,10,-1,9,8,3,4,-8,1,6,9,-9,2,-3,9,3,7,2,10,-5,-7,-2,4,6,9,-8],[-8,6,-4,6,0,3,-6,8,-3,-7,2,2,-4,-3,-3,-2,-7,0,-4,-3,3,-8,-10,1,7,3,7,-3,3,6,-7,8,-9,6,0,10,-4,9,-10,6,2,-3,5,7,-7,-4,-6,-7,4,9],[-9,-2,8,-9,4,-5,4,6,-7,-8,4,-4,10,-2,4,10,-9,2,-9,9,-10,-9,10,-10,-9,-7,-4,-2,0,-6,3,5,8,-6,9,7,3,6,-7,10,2,-8,0,0,0,5,8,2,7,-9],[8,3,6,0,5,-8,7,9,-9,1,-7,-2,4,5,-1,8,3,-7,-6,5,-7,7,-8,0,-5,9,-2,-4,-10,-9,-7,-4,-4,2,-4,0,3,-8,5,7,0,-8,2,8,8,-5,6,-9,-4,-4],[-9,-5,-5,3,-6,7,3,7,7,-8,9,3,9,7,-2,8,-8,9,9,8,-10,-5,10,-10,4,-1,8,-7,0,-9,5,-9,3,10,-1,-9,6,-9,9,6,-1,5,-9,5,-1,-1,6,8,-4,-7],[-3,10,-3,-10,-2,2,8,-4,7,6,5,0,-5,3,2,9,-7,-4,-2,7,4,-10,-1,-10,-5,1,1,3,-2,5,-3,-9,-7,-4,10,-3,-4,7,-8,-6,-4,9,-9,0,4,8,-5,6,0,-9],[-10,-7,-7,-9,-1,2,8,2,3,-3,10,9,-3,9,-7,-3,3,-9,2,-10,4,-2,0,9,2,5,-3,10,-2,4,9,2,-9,3,-9,-7,3,-4,6,-10,0,-10,-1,-1,6,5,1,-10,1,0],[3,-8,3,2,-9,3,9,-1,-9,6,-3,-6,-1,2,0,-5,2,8,3,6,-1,-4,1,0,-1,7,4,5,-9,7,-1,-1,0,1,-5,-10,-10,-3,-8,-4,-8,0,8,-1,1,-6,8,2,-2,-1],[3,9,4,3,1,-6,-10,-9,3,8,10,5,0,-7,2,-7,2,-1,-8,9,6,-1,-3,10,-5,-10,-5,-6,-8,4,4,-1,-3,-6,8,2,9,-6,5,6,2,-1,-9,-8,0,0,6,-4,9,-7],[-3,4,-9,2,-9,-2,-5,-9,-9,5,4,-9,-7,-10,7,4,8,-7,0,6,3,8,-6,-3,2,-10,3,6,1,10,-6,6,0,0,0,8,-9,1,-7,-6,10,2,1,8,-10,7,6,7,-8,-10],[2,9,-3,9,-2,10,-4,8,3,5,8,0,8,3,2,3,-6,10,5,5,6,1,0,9,-6,2,9,-9,-2,2,-1,-7,-9,-10,8,9,6,-7,-2,2,5,4,-7,6,7,-1,-5,10,7,-9],[7,5,-2,3,3,-1,3,-5,-9,-4,-4,0,-4,-5,-4,-10,7,-1,-3,-5,4,4,3,-1,-3,-5,0,-1,4,7,-2,8,3,5,7,4,8,-3,-4,-1,-7,4,4,0,-1,-3,-8,-2,-5,-1],[3,-7,-2,-8,1,-4,9,6,5,-3,-1,-9,1,0,-9,2,-4,-5,-5,-2,-7,-6,-2,9,-7,-8,-3,-10,7,10,5,0,5,-7,-3,-4,-6,-8,5,5,-10,-2,4,9,-3,4,-3,1,2,0],[-4,5,5,-9,-1,-9,9,-4,8,-4,-3,-5,-4,1,-8,-6,4,-6,8,-8,-4,-10,-7,-3,-9,-3,-2,1,-3,-3,-9,10,9,-9,9,-8,1,8,3,2,-3,8,-10,-7,2,-7,7,3,-5,-4],[-6,-8,2,-9,4,0,4,2,7,7,-2,-4,-2,6,4,3,-6,9,-4,-2,-2,10,5,6,3,-8,2,-2,-2,-4,-6,1,0,-2,9,6,-1,-1,-7,1,-5,3,-5,-7,10,2,2,-7,3,-1],[10,-4,10,8,-3,2,2,10,8,5,9,2,-1,-2,-2,7,-3,0,7,1,3,0,4,4,-5,-6,1,1,-7,7,2,-3,-5,6,-1,-3,-8,-2,4,-5,1,-7,7,-4,-9,6,-7,4,-4,-8],[7,0,5,-5,6,-7,-5,-6,-3,-7,-3,-3,-8,2,-3,-1,0,7,-3,5,-4,7,-7,-6,-5,7,10,-6,1,5,3,-9,5,1,3,-6,4,4,0,-3,-5,2,5,8,3,6,8,-2,-4,-3],[-10,6,-4,-7,6,10,4,-8,-4,2,4,-4,9,-6,1,-2,8,-6,8,0,8,1,0,0,-6,-6,-7,9,-9,-3,0,-2,-2,-10,-10,-8,10,8,-7,-2,7,-2,-2,-6,10,10,4,-3,3,0],[1,-5,4,10,8,9,-3,0,-8,3,3,-2,-5,9,-1,2,-4,-7,7,-10,-5,2,8,2,8,4,-4,-5,10,-8,2,10,-5,-8,4,-2,3,0,-6,0,0,8,5,-5,-5,-1,1,-9,0,-5],[8,7,9,-8,-8,3,7,3,8,-2,9,8,6,10,-8,6,-5,-7,2,-2,-7,-3,7,-4,-9,10,2,9,0,0,-2,6,-6,-2,6,5,7,0,0,-5,-9,-3,-3,-1,8,-5,5,-2,-6,9],[-9,-7,-3,6,10,8,8,-5,10,8,6,1,-1,3,-3,-5,6,9,6,2,1,-2,7,0,5,-10,10,-5,-7,3,8,-9,6,9,-7,-10,0,-7,6,-6,2,-4,-7,-6,-1,-3,-2,-3,1,-1],[-5,8,-5,-3,4,5,10,-2,6,5,-10,5,10,-4,-9,-4,4,6,-3,-5,10,1,-10,3,-8,-9,5,5,9,-5,-9,4,9,8,-3,10,-1,-2,-1,5,-3,-9,-1,-5,9,8,5,-2,-3,-6]]
(Je suis en 32 bits)
(ca le fait avec une version que j'ai compile, et avec le paquet debian aussi)
k:=50;A:=matrix(k,k,(i0,j0)->rand(21)-10):;
d1:=matdet(A,0)
d2:=det(A,bareiss);
d3:=det(A)
alors d1 est tres different des 2 autres.
en 0.8.6 les 3 coincident bien et c'est d1 qui semble faux.
Fred
Ex:
[[-3,10,-7,-2,0,-3,5,-8,-6,5,8,1,-9,9,-6,7,-4,3,-3,2,-9,-2,1,1,-7,-4,-3,-8,-5,-3,-8,1,-7,-1,-4,7,10,5,1,2,5,-6,8,8,-1,-4,-7,8,-5,-2],[-4,-7,2,-8,-2,4,-3,10,-4,0,-4,0,4,-6,3,-10,-3,-5,-3,2,7,1,-9,0,-1,-4,2,1,2,-6,7,-6,-1,3,7,9,-8,-7,0,-5,8,-3,-8,10,8,2,-9,-3,-4,-2],[2,10,6,3,9,-9,-10,-2,6,-2,-9,-2,-9,-4,1,-1,-1,1,4,-1,-9,-2,-6,9,4,-2,9,8,-6,0,3,-2,-7,-2,-9,-9,-2,10,0,-6,9,-1,4,-1,-1,3,8,1,3,-3],[-1,-3,10,-8,8,-7,4,-10,5,7,-2,-8,-3,-10,10,-10,7,-2,-5,3,-6,-3,-4,1,10,9,-8,1,-8,-1,-6,-8,-1,-2,1,5,-8,-5,-7,-1,1,6,-2,-5,10,-2,2,-5,8,7],[2,-9,-9,9,7,-8,-10,6,1,9,2,-9,-3,-8,-9,-3,6,8,-1,-8,6,9,-6,-1,-1,3,0,-4,9,-5,-4,9,-8,6,-1,-1,6,-8,-2,2,5,-9,-9,7,5,-6,-5,4,-9,8],[-9,-8,2,-9,-9,-8,-1,-4,-1,10,-1,0,8,10,8,5,-6,-9,-9,10,-7,-4,-5,9,5,6,4,4,0,-5,-4,-2,2,10,-7,-9,9,2,1,9,6,-10,-7,-8,5,-1,3,-1,5,-4],[-10,1,8,1,-6,6,-2,-5,7,-4,3,0,8,5,-8,-8,4,2,-9,-9,3,0,-9,5,1,2,2,3,1,7,9,-4,-1,-3,3,-6,-10,1,-5,-3,-2,-3,-6,-5,6,3,-4,-3,3,8],[-4,-3,10,5,-5,6,1,10,2,-8,-5,0,-8,2,5,-5,10,9,6,1,-8,-2,-7,5,-8,10,-4,-5,-5,-10,-6,0,10,5,-4,5,8,-3,3,8,0,-10,8,-4,-5,-4,5,4,10,4],[-10,4,-6,-1,-3,3,6,0,-1,-2,3,8,10,5,9,10,9,-8,8,-7,7,-4,-1,-10,-8,7,4,10,6,2,-6,10,5,-4,-4,8,3,10,-6,-5,2,-4,1,-3,-10,-7,-1,-4,4,10],[8,0,1,6,-7,-5,-7,-10,-5,4,5,-2,4,7,1,-5,4,4,-5,9,6,-8,3,4,3,-5,10,-9,-4,1,6,8,5,10,5,-1,7,-3,2,-2,-3,-9,9,-10,-9,5,-2,10,-10,-3],[10,-3,-1,0,-10,6,-10,8,8,-10,9,9,-10,-10,-8,6,10,5,-8,0,3,9,9,-7,-6,8,-4,6,8,-6,2,7,10,2,4,-4,3,-3,-4,-8,-3,-5,-10,-5,5,10,8,5,-7,6],[-5,-4,9,-8,-5,2,9,8,-9,-1,-10,-2,-6,1,9,-3,-2,9,9,3,4,-6,0,-1,10,6,4,8,2,9,6,-7,9,8,3,9,-4,-3,9,1,-2,6,-3,-2,-9,-1,-7,-10,-9,-6],[5,-3,5,-9,9,5,1,-5,-7,-4,-4,-7,6,7,-4,10,-2,7,-2,-7,1,4,7,-5,8,1,9,-10,-3,-3,1,2,0,5,2,-10,5,7,-8,10,9,-2,-2,2,4,-2,-9,3,-1,7],[-1,10,4,6,0,8,-7,5,-9,-6,5,3,-2,9,-7,0,8,8,1,9,-1,4,-5,6,-10,9,-1,8,-1,1,-4,0,6,-2,3,5,8,-5,-5,-4,-5,-1,10,5,-1,4,-8,10,6,10],[-3,3,8,4,10,1,2,-3,1,8,-9,7,-5,4,-3,2,-5,10,-9,-3,4,-8,4,6,-3,8,-5,5,-6,7,7,-8,8,-7,-6,5,-3,-5,-8,3,9,-2,-6,-6,9,10,-1,-10,9,0],[3,3,-1,2,-7,1,2,9,3,-2,10,-4,-5,3,-8,-6,7,-8,-7,-8,-9,9,-9,-10,1,4,0,-4,-10,8,-1,3,4,-5,-7,0,9,6,9,-7,-7,3,-6,-8,2,-1,5,9,-1,1],[4,5,1,10,9,6,-6,-10,1,-6,-1,4,-7,-3,5,3,5,-5,5,9,-10,-3,4,8,-2,2,-10,6,7,-10,-1,-3,6,-7,-10,10,-9,-5,9,-3,-4,3,2,1,10,-9,8,6,7,6],[7,9,-8,7,3,9,-2,-3,-6,1,10,2,3,-5,2,-10,-1,2,-3,-1,8,-5,-1,-8,7,-10,-6,4,3,-5,-2,-4,10,-5,10,-2,-10,-5,4,5,-6,-10,-5,-2,-4,2,-9,-1,-3,-9],[-4,9,6,9,-9,2,-10,-6,5,4,-7,-8,-10,6,-2,8,6,-4,-6,4,-9,-1,4,-9,-4,2,-4,-9,-3,-9,5,-7,8,6,6,-10,8,2,1,4,7,10,10,6,4,2,8,-6,-5,-1],[-5,4,-9,4,6,10,-6,7,8,10,-4,10,-7,10,2,0,3,3,10,-9,-3,1,-2,10,-8,-3,10,0,4,-6,3,-10,-1,-8,0,-5,-9,3,0,-5,10,-9,-6,6,1,-4,-3,-7,0,-1],[-8,0,5,-3,6,2,-7,-4,9,10,1,-2,4,9,8,-3,8,-3,5,-3,-2,3,6,3,7,3,-1,-4,4,-3,-9,2,-2,9,2,-6,6,-5,-6,5,-10,7,5,4,1,-2,8,-4,10,5],[-2,10,-2,-8,-6,-1,-1,-9,7,5,1,-9,0,5,0,6,-4,3,7,8,7,-10,-3,-1,10,0,10,-5,-1,-2,-3,2,-6,7,3,-4,0,6,7,-4,-5,5,-3,8,1,10,1,-8,6,7],[4,10,-10,-4,-1,3,1,8,-6,6,10,0,6,5,5,4,5,-1,0,-2,-3,8,4,-9,-6,8,-3,-10,-5,-3,-2,-5,6,-3,-2,-9,-4,7,-6,-7,10,7,-1,10,6,-3,-8,10,6,-6],[1,1,4,8,2,8,-1,1,2,9,6,1,-10,-9,8,0,-9,-4,0,8,-7,-1,0,3,-3,5,-6,-6,-10,-1,-4,5,-1,-10,-6,-10,-9,9,-8,9,1,-1,5,-8,10,7,-2,-3,0,2],[-10,-8,-5,-8,-1,-9,0,-5,0,1,1,0,0,5,0,-8,3,-3,1,-2,4,3,7,3,-4,-8,-5,2,-7,-5,8,8,-1,-2,3,-9,-8,-3,-1,-4,0,-10,-3,9,9,-3,0,0,-4,-1],[-9,-10,-8,5,3,-7,-8,10,-8,-10,0,-5,3,5,6,-8,-7,8,-7,-3,-9,-3,4,-10,1,-9,3,-4,7,1,10,0,-8,-10,8,5,4,-6,-2,5,1,-9,-5,5,10,-7,2,-7,-10,-1],[8,-4,7,-5,10,-7,6,5,-4,5,-5,-10,9,-8,4,-10,9,3,-7,-3,6,7,7,-1,7,1,-3,4,8,7,10,-6,9,-8,4,-1,5,4,-8,-4,10,-3,-8,-8,10,-3,-3,-4,5,-9],[6,7,-2,-5,9,-7,-10,-6,8,-3,-6,-10,-5,-7,1,10,-9,7,8,-1,-5,2,-5,-5,-6,3,-3,7,-10,1,-3,-9,4,-7,9,1,-7,7,-7,2,-8,2,2,-8,4,5,-1,4,9,5],[8,1,-4,8,-5,-6,8,9,8,5,-6,9,-7,-7,2,3,1,4,-2,0,-7,-7,-7,1,2,10,-1,9,8,3,4,-8,1,6,9,-9,2,-3,9,3,7,2,10,-5,-7,-2,4,6,9,-8],[-8,6,-4,6,0,3,-6,8,-3,-7,2,2,-4,-3,-3,-2,-7,0,-4,-3,3,-8,-10,1,7,3,7,-3,3,6,-7,8,-9,6,0,10,-4,9,-10,6,2,-3,5,7,-7,-4,-6,-7,4,9],[-9,-2,8,-9,4,-5,4,6,-7,-8,4,-4,10,-2,4,10,-9,2,-9,9,-10,-9,10,-10,-9,-7,-4,-2,0,-6,3,5,8,-6,9,7,3,6,-7,10,2,-8,0,0,0,5,8,2,7,-9],[8,3,6,0,5,-8,7,9,-9,1,-7,-2,4,5,-1,8,3,-7,-6,5,-7,7,-8,0,-5,9,-2,-4,-10,-9,-7,-4,-4,2,-4,0,3,-8,5,7,0,-8,2,8,8,-5,6,-9,-4,-4],[-9,-5,-5,3,-6,7,3,7,7,-8,9,3,9,7,-2,8,-8,9,9,8,-10,-5,10,-10,4,-1,8,-7,0,-9,5,-9,3,10,-1,-9,6,-9,9,6,-1,5,-9,5,-1,-1,6,8,-4,-7],[-3,10,-3,-10,-2,2,8,-4,7,6,5,0,-5,3,2,9,-7,-4,-2,7,4,-10,-1,-10,-5,1,1,3,-2,5,-3,-9,-7,-4,10,-3,-4,7,-8,-6,-4,9,-9,0,4,8,-5,6,0,-9],[-10,-7,-7,-9,-1,2,8,2,3,-3,10,9,-3,9,-7,-3,3,-9,2,-10,4,-2,0,9,2,5,-3,10,-2,4,9,2,-9,3,-9,-7,3,-4,6,-10,0,-10,-1,-1,6,5,1,-10,1,0],[3,-8,3,2,-9,3,9,-1,-9,6,-3,-6,-1,2,0,-5,2,8,3,6,-1,-4,1,0,-1,7,4,5,-9,7,-1,-1,0,1,-5,-10,-10,-3,-8,-4,-8,0,8,-1,1,-6,8,2,-2,-1],[3,9,4,3,1,-6,-10,-9,3,8,10,5,0,-7,2,-7,2,-1,-8,9,6,-1,-3,10,-5,-10,-5,-6,-8,4,4,-1,-3,-6,8,2,9,-6,5,6,2,-1,-9,-8,0,0,6,-4,9,-7],[-3,4,-9,2,-9,-2,-5,-9,-9,5,4,-9,-7,-10,7,4,8,-7,0,6,3,8,-6,-3,2,-10,3,6,1,10,-6,6,0,0,0,8,-9,1,-7,-6,10,2,1,8,-10,7,6,7,-8,-10],[2,9,-3,9,-2,10,-4,8,3,5,8,0,8,3,2,3,-6,10,5,5,6,1,0,9,-6,2,9,-9,-2,2,-1,-7,-9,-10,8,9,6,-7,-2,2,5,4,-7,6,7,-1,-5,10,7,-9],[7,5,-2,3,3,-1,3,-5,-9,-4,-4,0,-4,-5,-4,-10,7,-1,-3,-5,4,4,3,-1,-3,-5,0,-1,4,7,-2,8,3,5,7,4,8,-3,-4,-1,-7,4,4,0,-1,-3,-8,-2,-5,-1],[3,-7,-2,-8,1,-4,9,6,5,-3,-1,-9,1,0,-9,2,-4,-5,-5,-2,-7,-6,-2,9,-7,-8,-3,-10,7,10,5,0,5,-7,-3,-4,-6,-8,5,5,-10,-2,4,9,-3,4,-3,1,2,0],[-4,5,5,-9,-1,-9,9,-4,8,-4,-3,-5,-4,1,-8,-6,4,-6,8,-8,-4,-10,-7,-3,-9,-3,-2,1,-3,-3,-9,10,9,-9,9,-8,1,8,3,2,-3,8,-10,-7,2,-7,7,3,-5,-4],[-6,-8,2,-9,4,0,4,2,7,7,-2,-4,-2,6,4,3,-6,9,-4,-2,-2,10,5,6,3,-8,2,-2,-2,-4,-6,1,0,-2,9,6,-1,-1,-7,1,-5,3,-5,-7,10,2,2,-7,3,-1],[10,-4,10,8,-3,2,2,10,8,5,9,2,-1,-2,-2,7,-3,0,7,1,3,0,4,4,-5,-6,1,1,-7,7,2,-3,-5,6,-1,-3,-8,-2,4,-5,1,-7,7,-4,-9,6,-7,4,-4,-8],[7,0,5,-5,6,-7,-5,-6,-3,-7,-3,-3,-8,2,-3,-1,0,7,-3,5,-4,7,-7,-6,-5,7,10,-6,1,5,3,-9,5,1,3,-6,4,4,0,-3,-5,2,5,8,3,6,8,-2,-4,-3],[-10,6,-4,-7,6,10,4,-8,-4,2,4,-4,9,-6,1,-2,8,-6,8,0,8,1,0,0,-6,-6,-7,9,-9,-3,0,-2,-2,-10,-10,-8,10,8,-7,-2,7,-2,-2,-6,10,10,4,-3,3,0],[1,-5,4,10,8,9,-3,0,-8,3,3,-2,-5,9,-1,2,-4,-7,7,-10,-5,2,8,2,8,4,-4,-5,10,-8,2,10,-5,-8,4,-2,3,0,-6,0,0,8,5,-5,-5,-1,1,-9,0,-5],[8,7,9,-8,-8,3,7,3,8,-2,9,8,6,10,-8,6,-5,-7,2,-2,-7,-3,7,-4,-9,10,2,9,0,0,-2,6,-6,-2,6,5,7,0,0,-5,-9,-3,-3,-1,8,-5,5,-2,-6,9],[-9,-7,-3,6,10,8,8,-5,10,8,6,1,-1,3,-3,-5,6,9,6,2,1,-2,7,0,5,-10,10,-5,-7,3,8,-9,6,9,-7,-10,0,-7,6,-6,2,-4,-7,-6,-1,-3,-2,-3,1,-1],[-5,8,-5,-3,4,5,10,-2,6,5,-10,5,10,-4,-9,-4,4,6,-3,-5,10,1,-10,3,-8,-9,5,5,9,-5,-9,4,9,8,-3,10,-1,-2,-1,5,-3,-9,-1,-5,9,8,5,-2,-3,-6]]
Re: pari matdet
chez moi c'est identique sur ton exemple
7>> matdet(A)
1046088311290524516321120623795924354054000507606011094835193721914521
// Time 0.02
8>> det(A)
1046088311290524516321120623795924354054000507606011094835193721914521
// Time 0.04
9>> det(A,bareiss)
1046088311290524516321120623795924354054000507606011094835193721914521
Quelle version de pari utilises-tu?
7>> matdet(A)
1046088311290524516321120623795924354054000507606011094835193721914521
// Time 0.02
8>> det(A)
1046088311290524516321120623795924354054000507606011094835193721914521
// Time 0.04
9>> det(A,bareiss)
1046088311290524516321120623795924354054000507606011094835193721914521
Quelle version de pari utilises-tu?
-
- Messages : 1137
- Inscription : dim. mai 20, 2007 7:09 am
- Localisation : Paris
- Contact :
Re: pari matdet
J'ai une debian stable, mais
xcas 0.9.0 est lie en statique sur pari
NB: sur une freebsd amd64 je n'ai pas reussi a avoir de PB, uniquement sur mes 32bits.
ca le fait avec le paquet suivant:
dpkg -s giac
Package: giac
Status: install ok installed
Priority: optional
Section: math
Installed-Size: 70819
Maintainer: Carleos Artime <carleos@pinon.ccu.uniovi.es>
Architecture: i386
Version: 0.9.0-2
ldd /usr/bin/xcas
linux-gate.so.1 => (0xb7f14000)
libgiac.so.0 => /usr/lib/libgiac.so.0 (0xb6fd3000)
libjpeg.so.62 => /usr/lib/libjpeg.so.62 (0xb6fb4000)
libz.so.1 => /usr/lib/libz.so.1 (0xb6f9e000)
libXext.so.6 => /usr/lib/libXext.so.6 (0xb6f90000)
libXft.so.2 => /usr/lib/libXft.so.2 (0xb6f7d000)
libfontconfig.so.1 => /usr/lib/libfontconfig.so.1 (0xb6f52000)
libXinerama.so.1 => /usr/lib/libXinerama.so.1 (0xb6f4f000)
libGL.so.1 => /usr/lib/libGL.so.1 (0xb6eec000)
librt.so.1 => /lib/i686/cmov/librt.so.1 (0xb6ee3000)
libpthread.so.0 => /lib/i686/cmov/libpthread.so.0 (0xb6ec9000)
libdl.so.2 => /lib/i686/cmov/libdl.so.2 (0xb6ec5000)
libpng12.so.0 => /usr/lib/libpng12.so.0 (0xb6ea1000)
libgmp.so.3 => /usr/lib/libgmp.so.3 (0xb6e5e000)
libstdc++.so.6 => /usr/lib/libstdc++.so.6 (0xb6d70000)
libm.so.6 => /lib/i686/cmov/libm.so.6 (0xb6d4a000)
libgcc_s.so.1 => /lib/libgcc_s.so.1 (0xb6d3c000)
libc.so.6 => /lib/i686/cmov/libc.so.6 (0xb6be1000)
libX11.so.6 => /usr/lib/libX11.so.6 (0xb6af2000)
libXau.so.6 => /usr/lib/libXau.so.6 (0xb6aef000)
libfreetype.so.6 => /usr/lib/libfreetype.so.6 (0xb6a7a000)
libXrender.so.1 => /usr/lib/libXrender.so.1 (0xb6a70000)
libexpat.so.1 => /usr/lib/libexpat.so.1 (0xb6a4a000)
libXxf86vm.so.1 => /usr/lib/libXxf86vm.so.1 (0xb6a45000)
libXdamage.so.1 => /usr/lib/libXdamage.so.1 (0xb6a42000)
libXfixes.so.3 => /usr/lib/libXfixes.so.3 (0xb6a3d000)
libdrm.so.2 => /usr/lib/libdrm.so.2 (0xb6a34000)
/lib/ld-linux.so.2 (0xb7f15000)
libxcb-xlib.so.0 => /usr/lib/libxcb-xlib.so.0 (0xb6a32000)
libxcb.so.1 => /usr/lib/libxcb.so.1 (0xb6a1a000)
libXdmcp.so.6 => /usr/lib/libXdmcp.so.6 (0xb6a15000)
Fred
xcas 0.9.0 est lie en statique sur pari
NB: sur une freebsd amd64 je n'ai pas reussi a avoir de PB, uniquement sur mes 32bits.
ca le fait avec le paquet suivant:
dpkg -s giac
Package: giac
Status: install ok installed
Priority: optional
Section: math
Installed-Size: 70819
Maintainer: Carleos Artime <carleos@pinon.ccu.uniovi.es>
Architecture: i386
Version: 0.9.0-2
ldd /usr/bin/xcas
linux-gate.so.1 => (0xb7f14000)
libgiac.so.0 => /usr/lib/libgiac.so.0 (0xb6fd3000)
libjpeg.so.62 => /usr/lib/libjpeg.so.62 (0xb6fb4000)
libz.so.1 => /usr/lib/libz.so.1 (0xb6f9e000)
libXext.so.6 => /usr/lib/libXext.so.6 (0xb6f90000)
libXft.so.2 => /usr/lib/libXft.so.2 (0xb6f7d000)
libfontconfig.so.1 => /usr/lib/libfontconfig.so.1 (0xb6f52000)
libXinerama.so.1 => /usr/lib/libXinerama.so.1 (0xb6f4f000)
libGL.so.1 => /usr/lib/libGL.so.1 (0xb6eec000)
librt.so.1 => /lib/i686/cmov/librt.so.1 (0xb6ee3000)
libpthread.so.0 => /lib/i686/cmov/libpthread.so.0 (0xb6ec9000)
libdl.so.2 => /lib/i686/cmov/libdl.so.2 (0xb6ec5000)
libpng12.so.0 => /usr/lib/libpng12.so.0 (0xb6ea1000)
libgmp.so.3 => /usr/lib/libgmp.so.3 (0xb6e5e000)
libstdc++.so.6 => /usr/lib/libstdc++.so.6 (0xb6d70000)
libm.so.6 => /lib/i686/cmov/libm.so.6 (0xb6d4a000)
libgcc_s.so.1 => /lib/libgcc_s.so.1 (0xb6d3c000)
libc.so.6 => /lib/i686/cmov/libc.so.6 (0xb6be1000)
libX11.so.6 => /usr/lib/libX11.so.6 (0xb6af2000)
libXau.so.6 => /usr/lib/libXau.so.6 (0xb6aef000)
libfreetype.so.6 => /usr/lib/libfreetype.so.6 (0xb6a7a000)
libXrender.so.1 => /usr/lib/libXrender.so.1 (0xb6a70000)
libexpat.so.1 => /usr/lib/libexpat.so.1 (0xb6a4a000)
libXxf86vm.so.1 => /usr/lib/libXxf86vm.so.1 (0xb6a45000)
libXdamage.so.1 => /usr/lib/libXdamage.so.1 (0xb6a42000)
libXfixes.so.3 => /usr/lib/libXfixes.so.3 (0xb6a3d000)
libdrm.so.2 => /usr/lib/libdrm.so.2 (0xb6a34000)
/lib/ld-linux.so.2 (0xb7f15000)
libxcb-xlib.so.0 => /usr/lib/libxcb-xlib.so.0 (0xb6a32000)
libxcb.so.1 => /usr/lib/libxcb.so.1 (0xb6a1a000)
libXdmcp.so.6 => /usr/lib/libXdmcp.so.6 (0xb6a15000)
Fred
Re: pari matdet
Oui, j'observe ca aussi avec la version du package, je fais un package 0.9.1 que je vais mettre dans testing, on verra si ca continue a poser probleme (je soupconne un probleme de conversion entre les GEN de pari et mes gen).
Re: pari matdet
bon, c'est pas mieux avec le 0.9.1. Et ca deconne uniquement avec la version dynamique, avec la version statique compilee a partir des memes fichiers objets ca marche, mystere!
-
- Messages : 1137
- Inscription : dim. mai 20, 2007 7:09 am
- Localisation : Paris
- Contact :
Re: pari matdet
Moi j'avais le bug avec des versions ou pari est en statique.
Du coup j'ai recompile pari en statique en ajoutant un -fPIC dans le Makefile
matdet a l'air correct. (0.9.1)
pour la 0.9.1 j'avais toujours le plantage exp(2*i*Pi/3) sous xcas
J'ai recompile avec
CXXFLAGS = -g -O2 -fno-strict-aliasing -DGIAC_GENERIC_CONSTANTS
et ce bug a disparu aussi
Fred
Du coup j'ai recompile pari en statique en ajoutant un -fPIC dans le Makefile
matdet a l'air correct. (0.9.1)
pour la 0.9.1 j'avais toujours le plantage exp(2*i*Pi/3) sous xcas
J'ai recompile avec
CXXFLAGS = -g -O2 -fno-strict-aliasing -DGIAC_GENERIC_CONSTANTS
et ce bug a disparu aussi
Fred
Re: pari matdet
salut,
chez moi, exp(2*i*pi/3) fonctionne normalement. Je vais recompiler pari avec -fPIC et voir si ca resout le probleme!
chez moi, exp(2*i*pi/3) fonctionne normalement. Je vais recompiler pari avec -fPIC et voir si ca resout le probleme!
Re: pari matdet
Apparamment, ce n'est pas ca, mais c'est parce que j'ai 2 libpari.a une dans /usr/lib et une dans /usr/local/lib, differentes, et le link ne se faisait pas avec la meme. J'ecrase la "mauvaise".
-
- Messages : 1137
- Inscription : dim. mai 20, 2007 7:09 am
- Localisation : Paris
- Contact :
Re: pari matdet
ce que je n'ai pas compris pour exp(2*i*Pi/3) c'est que avec giac ca ne plantait pas, c'est juste avec xcas...
fred
fred
Re: pari matdet
avec liaison dynamique dans les 2 cas?
-
- Messages : 1137
- Inscription : dim. mai 20, 2007 7:09 am
- Localisation : Paris
- Contact :
Re: pari matdet
Non, dans les 2 cas sont lies en statique sur pas mal de librairies dont: libgmp, pari, mpfr, ntl, gsl
Fred
Fred
-
- Messages : 1137
- Inscription : dim. mai 20, 2007 7:09 am
- Localisation : Paris
- Contact :
Re: pari matdet
avec la version du 29/11 je n'ai eu aucun de tous ces Pb. (sans mettre de FLAGS particuliers pour la compilation)
J'ai mis une version xcas_user 32bits avec pas mal de librairies statiques avec l'espoir d'avoir une grande compatibilite avec les divers linux.
http://people.math.jussieu.fr/~han/xcas ... inux32.tgz
(c'est un peu gros je ne sais pas trop comment le faire maigrir)
a+
Fred
J'ai mis une version xcas_user 32bits avec pas mal de librairies statiques avec l'espoir d'avoir une grande compatibilite avec les divers linux.
http://people.math.jussieu.fr/~han/xcas ... inux32.tgz
(c'est un peu gros je ne sais pas trop comment le faire maigrir)
a+
Fred
Re: pari matdet
J'ai rajouté automatiquement dans le configure -DGIAC_GENERIC_CONSTANTS, ce qui doit expliquer que ça marche. Super pour le binaire linux 32 bits, il faudra que je le teste, d'après mon expérience le problème principal c'était la compatibilité de la libc.
Sinon, pour info, dans la rubrique bugs corrigés en 0.9.1 mais pas dans la 0.9.0, bug repéré par mes étudiants: il y a le calcul de déterminant modulo p lorsque p>2^31. Mais je pense que d'ici les oraux de cette année, on aura une version 0.9.1 ou plus suffisamment stabilisée.
Sinon, pour info, dans la rubrique bugs corrigés en 0.9.1 mais pas dans la 0.9.0, bug repéré par mes étudiants: il y a le calcul de déterminant modulo p lorsque p>2^31. Mais je pense que d'ici les oraux de cette année, on aura une version 0.9.1 ou plus suffisamment stabilisée.
-
- Messages : 1137
- Inscription : dim. mai 20, 2007 7:09 am
- Localisation : Paris
- Contact :
Re: pari matdet
oui j'ai pris une centos car les versions des librairies sont un poil plus vieilles que sur debian stable du coup ton rpm ne marche pas sur ces redhat.
merci pour le bug du det, mais en fait j'ai remarque que c'etait plus rapide dans Z:
p:=nextprime(2^31);
n:=100;A:=matrix(n,n,(k,l)->rand(60)); B:=A % p;
b:=B*B; // 7 secondes
a:=(A*A) % p; // instantane
simplify(b-a) donne bien 0.
det(B); // 1.65s
det(A) % p ; // instantane et meme reponse en 0.9.1
Fred
merci pour le bug du det, mais en fait j'ai remarque que c'etait plus rapide dans Z:
p:=nextprime(2^31);
n:=100;A:=matrix(n,n,(k,l)->rand(60)); B:=A % p;
b:=B*B; // 7 secondes
a:=(A*A) % p; // instantane
simplify(b-a) donne bien 0.
det(B); // 1.65s
det(A) % p ; // instantane et meme reponse en 0.9.1
Fred
Re: pari matdet
En effet, les calculs modulaires sont faits avec des entiers longs, alors que sur Z c'est uniquement des entiers machines sauf à la reconstruction.