summaryrefslogtreecommitdiffstats
path: root/divers
diff options
context:
space:
mode:
authorDavid A. Madore <david+git@madore.org>2011-11-27 08:21:49 +0100
committerDavid A. Madore <david+git@madore.org>2011-11-27 08:21:49 +0100
commit3488be4c25d9a8bc59952e84e2fb1179798a84bf (patch)
tree471e79faa4029e3c4e02f5faee18867bec99750f /divers
parente7b595293338f69f9ee7b7fe7391a7306a5ca551 (diff)
downloadgalois-3488be4c25d9a8bc59952e84e2fb1179798a84bf.zip
galois-3488be4c25d9a8bc59952e84e2fb1179798a84bf.tar.gz
galois-3488be4c25d9a8bc59952e84e2fb1179798a84bf.tar.bz2
sageries: une correction + je continue le calcul un peu plus loin.
Diffstat (limited to 'divers')
-rw-r--r--divers/sageries/resolvante-d5-dans-m2010
1 files changed, 8 insertions, 2 deletions
diff --git a/divers/sageries/resolvante-d5-dans-m20 b/divers/sageries/resolvante-d5-dans-m20
index 2277bea..be91ce2 100644
--- a/divers/sageries/resolvante-d5-dans-m20
+++ b/divers/sageries/resolvante-d5-dans-m20
@@ -61,7 +61,7 @@ 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)
+dexpr = (z1-z2)*(z1-z3)*(z1-z4)*(z1-z5)*(z2-z3)*(z2-z4)*(z2-z5)*(z3-z4)*(z3-z5)*(z4-z5)
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)]
@@ -79,4 +79,10 @@ 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()
+ results[i][j] = (-1)^(i+j)*(qsubmat[complements[j],complements[i]].determinant())
+
+qsubadj = matrix(R, 12, 12, results)
+tmp = [sum([pdsubmat[1][k]*qsubadj[k][j] for k in range(12)]) for j in range(12)]
+subfactor = -tmp[11]/1280
+qsubdetfact = R(qsubdet/subfactor)
+donq = [R(pol/subfactor)/qsubdetfact for pol in tmp]