cyclic6 et giac 1.5.0-37

Bugs

Modérateur : xcasadmin

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

cyclic6 et giac 1.5.0-37

Message par frederic han » dim. févr. 17, 2019 8:20 pm

Bonjour Bernard,
avec les sources 1.5.0-37 j'obtiens un résultat different d'avec une version de giac plus ancienne 1.4.9 pour une base de groebner de cyclic6 modulo
p=2**31-1
(les termes B[24] et B[40] sont differents.

Code : Tout sélectionner

p:=2**31-1;
I:=[x0 + x1 + x2 + x3 + x4 + x5, x0*x1 + x1*x2 + x2*x3 + x3*x4 + x0*x5 + x4*x5, x0*x1*x2 + x1*x2*x3 + x2*x3*x4 + x0*x1*x5 + x0*x4*x5 + x3*x4*x5, x0*x1*x2*x3 + x1*x2*x3*x4 + x0*x1*x2*x5 + x0*x1*x4*x5 + x0*x3*x4*x5 + x2*x3*x4*x5, x0*x1*x2*x3*x4 + x0*x1*x2*x3*x5 + x0*x1*x2*x4*x5 + x0*x1*x3*x4*x5 + x0*x2*x3*x4*x5 + x1*x2*x3*x4*x5, x0*x1*x2*x3*x4*x5 - 1];
B:=gbasis(I % p,[x0, x1, x2, x3, x4, x5]);
old40:=x1*x2*x3 + x2*x3^2 + 1073741821*x1*x2*x4 + 1073741822*x2^2*x4 + x1*x3*x4 + 1073741823*x2*x3*x4 + 2*x3^2*x4 + 1073741822*x2*x4^2 + x3*x4^2 + x1*x2*x5 + x2^2*x5 - 1073741823*x1*x3*x5 - x2*x3*x5 - x3^2*x5 - x1*x4*x5 + 1073741823*x2*x4*x5 + x3*x4*x5 + x2*x5^2 - x3*x5^2;
print(version())
print(normal(B[40]-old40 % p));

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

Re: cyclic6 et giac 1.5.0-37

Message par parisse » lun. févr. 18, 2019 4:09 pm

Un probleme d'overflow, voila le correctif
https://dev.geogebra.org/trac/changeset/66972/
Ca ne touche que les calculs modulo un premier compris entre 2^30 et 2^31. Comme ils ne sont pas utilises dans les calculs sur Q, ca ne devrait pas trop deranger les utilisateurs de Xcas. J'attends toujours de valider le nouveau code de bases de Groebner pour faire une mise a jour, j'aimerais bien calculer cyclic10 sur Q, mais ca necessite beaucoup de memoire (256G au moins) et du temps (probablement une dizaine de jours sur une architecture avec 36 CPU).

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

Re: cyclic6 et giac 1.5.0-37

Message par parisse » jeu. févr. 21, 2019 12:53 pm

Bon, j'ai fait une mise a jour (-43) sans attendre de savoir si cyclic10 fonctionne, les autres tests ont l'air bons.

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

Re: cyclic6 et giac 1.5.0-37

Message par frederic han » ven. févr. 22, 2019 10:00 pm

OK merci, en fait j'ai mis ce petit patch dans le package sage, il a deja eu une revue positive.

Fred

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

Re: cyclic6 et giac 1.5.0-37

Message par parisse » sam. févr. 23, 2019 7:54 am

Oui. Mais en voyant le nombre de messages echanges pour ce correctif, je me dis quand meme que les processus utilises par sage sont bien lourds: il m'a fallu environ 1 heure pour corriger (ce n'etait pas un bug trivial, heureusement j'ai pu le localiser avec le log du SVN de geogebra), il en a surement fallu plus pour rediger tous ces messages.

Répondre