pari matdet

Xcas devel: interface utilisateur/user interface

Modérateur : xcasadmin

frederic han
Messages : 1137
Inscription : dim. mai 20, 2007 7:09 am
Localisation : Paris
Contact :

pari matdet

Message par frederic han » jeu. nov. 25, 2010 8:09 am

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]]

parisse
Messages : 5731
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: pari matdet

Message par parisse » jeu. nov. 25, 2010 10:30 am

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?

frederic han
Messages : 1137
Inscription : dim. mai 20, 2007 7:09 am
Localisation : Paris
Contact :

Re: pari matdet

Message par frederic han » jeu. nov. 25, 2010 12:04 pm

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

parisse
Messages : 5731
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: pari matdet

Message par parisse » jeu. nov. 25, 2010 12:47 pm

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).

parisse
Messages : 5731
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: pari matdet

Message par parisse » jeu. nov. 25, 2010 1:34 pm

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!

frederic han
Messages : 1137
Inscription : dim. mai 20, 2007 7:09 am
Localisation : Paris
Contact :

Re: pari matdet

Message par frederic han » jeu. nov. 25, 2010 9:57 pm

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

parisse
Messages : 5731
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: pari matdet

Message par parisse » ven. nov. 26, 2010 8:36 am

salut,

chez moi, exp(2*i*pi/3) fonctionne normalement. Je vais recompiler pari avec -fPIC et voir si ca resout le probleme!

parisse
Messages : 5731
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: pari matdet

Message par parisse » ven. nov. 26, 2010 8:45 am

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".

frederic han
Messages : 1137
Inscription : dim. mai 20, 2007 7:09 am
Localisation : Paris
Contact :

Re: pari matdet

Message par frederic han » ven. nov. 26, 2010 1:29 pm

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

parisse
Messages : 5731
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: pari matdet

Message par parisse » ven. nov. 26, 2010 3:11 pm

avec liaison dynamique dans les 2 cas?

frederic han
Messages : 1137
Inscription : dim. mai 20, 2007 7:09 am
Localisation : Paris
Contact :

Re: pari matdet

Message par frederic han » ven. nov. 26, 2010 3:16 pm

Non, dans les 2 cas sont lies en statique sur pas mal de librairies dont: libgmp, pari, mpfr, ntl, gsl

Fred

frederic han
Messages : 1137
Inscription : dim. mai 20, 2007 7:09 am
Localisation : Paris
Contact :

Re: pari matdet

Message par frederic han » lun. nov. 29, 2010 8:00 pm

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

parisse
Messages : 5731
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: pari matdet

Message par parisse » lun. nov. 29, 2010 8:46 pm

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.

frederic han
Messages : 1137
Inscription : dim. mai 20, 2007 7:09 am
Localisation : Paris
Contact :

Re: pari matdet

Message par frederic han » lun. nov. 29, 2010 9:09 pm

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

parisse
Messages : 5731
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: pari matdet

Message par parisse » mar. nov. 30, 2010 7:26 am

En effet, les calculs modulaires sont faits avec des entiers longs, alors que sur Z c'est uniquement des entiers machines sauf à la reconstruction.

Répondre