
(* lambda_14 *)
<<cox.m
Print["\n\n             SALEM 14  \n"]
DiStrat[salem14]


(*---------------------------------Salem factor---------------------------*)
(* Twist at 5, 13, 5x19, 5x41 *)
Print["\n\n             Obstruction over 5x19:  \n"]

s = salem14
qs = QEuler[s]
fs = Mat[s]

prime=5
ps = PrimeTwists[qs,fs,prime]
ps = Select[ps,CycloDim[GluePeriod[QTwist[qs,fs,#],fs]]<= 8&]
p5=ps[[1]]

prime=19
ps = PrimeTwists[qs,fs,prime]
ps = Select[ps,CycloDim[GluePeriod[QTwist[qs,fs,#],fs]]<= 8&]
p19 = ps[[1]]

qs = QTwist[qs,fs,p5 p19]
DiGlue[qs,fs]

c6 = Cyclotomic[6,x]
q6 = 5 QEuler[c6]
f6 = Mat[c6]

c9 = Cyclotomic[9,x]
q9 = QEuler[c9]
f9 = Mat[c9]
(* ps = PrimeTwists[q9,f9,19] *)
ps = {-1 - 2*x, -3 + 2*x^2, -3 + x}

poly9 = ps[[1]]
q9 = QTwist[q9,f9,poly9]

DiF[q6,f6]
DiF[q9,f9]

Print["Gluing is impossible because of different actions mod 3"]
