cyclic6 et giac 1.5.0-37
Publié : 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.
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));