summaryrefslogtreecommitdiffstats
path: root/divers
diff options
context:
space:
mode:
authorDavid A. Madore <david+git@madore.org>2011-11-26 21:32:06 (GMT)
committerDavid A. Madore <david+git@madore.org>2011-11-26 21:32:06 (GMT)
commite7b595293338f69f9ee7b7fe7391a7306a5ca551 (patch)
tree807315d78444ee980531c3527aeca9e1d4d0b18d /divers
parentc0e24d5a23a98d0ce0961b2422465827e92501c2 (diff)
downloadgalois-e7b595293338f69f9ee7b7fe7391a7306a5ca551.zip
galois-e7b595293338f69f9ee7b7fe7391a7306a5ca551.tar.gz
galois-e7b595293338f69f9ee7b7fe7391a7306a5ca551.tar.bz2
Amélioration des sageries, mais toujours sans succès...
Diffstat (limited to 'divers')
-rw-r--r--divers/sageries/resolvante-d5-dans-m2016
1 files changed, 12 insertions, 4 deletions
diff --git a/divers/sageries/resolvante-d5-dans-m20 b/divers/sageries/resolvante-d5-dans-m20
index 2d15b9f..2277bea 100644
--- a/divers/sageries/resolvante-d5-dans-m20
+++ b/divers/sageries/resolvante-d5-dans-m20
@@ -58,17 +58,25 @@ asym3 = -e([3]).expand(5).subs(x0=z1,x1=z2,x2=z3,x3=z4,x4=z5)
asym4 = +e([4]).expand(5).subs(x0=z1,x1=z2,x2=z3,x3=z4,x4=z5)
asym5 = -e([5]).expand(5).subs(x0=z1,x1=z2,x2=z3,x3=z4,x4=z5)
Isym = R.ideal([a1-asym1, a2-asym2, a3-asym3, a4-asym4, a5-asym5])
+B = Isym.groebner_basis()
pexpr = z1^2*(z2*z5 + z3*z4) + z2^2*(z1*z3 + z4*z5) + z3^2*(z1*z5 + z2*z4) + z4^2*(z1*z2 + z3*z5) + z5^2*(z1*z4 + z2*z3)
qexpr = z1*z2 + z2*z3 + z3*z4 + z4*z5 + z1*z5
dexpr = (z1-z2)*(z1-z3)*(z1-z4)*(z2-z3)*(z2-z4)*(z3-z4)
-B = Isym.groebner_basis()
-def debase(r):
+def debase(poly):
+ r = poly.reduce(B)
return [r.coefficient({z1:0,z2:i2,z3:i3,z4:i4,z5:i5}) for i2 in range(2) for i3 in range(3) for i4 in range(4) for i5 in range(5)]
-qtab = [debase((qexpr^i).reduce(B)) for i in range(12)]
-pdtab = [debase((pexpr^i*dexpr^j).reduce(B)) for i in range(6) for j in range(2)]
+qtab = [debase(qexpr^i) for i in range(12)]
+pdtab = [debase(pexpr^i*dexpr^j) for i in range(6) for j in range(2)]
qmat = matrix(R, 12, 120, qtab)
pdmat = matrix(R, 12, 120, pdtab)
cols = [0, 10, 105, 84, 109, 119, 57, 91, 111, 90, 110, 15]
qsubmat = qmat[range(12),cols]
pdsubmat = pdmat[range(12),cols]
+qsubdet = qsubmat.determinant()
+complements = [[i for i in range(12) if i!=j] for j in range(12)]
+results = [[None for j in range(12)] for i in range(12)]
+for i in range(12):
+ for j in range(12):
+ print (i,j)
+ results[i][j] = qsubmat[complements[i],complements[j]].determinant()