Merge branch 'main' of https://github.com/kedlaya-s-algorithm-research/griffiths-dwork-construction

Author: MaximusMellberg

experiments/ZetaData.jl

@@ -17,3 +17,19 @@ function data_fermat_curves()
         end
     end
 end
+
+function data_fermat_deform_6()
+
+    #ps = [7,11,13,17,19,23,29,31,37,41]
+    ps = [11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97]
+
+    d = 6
+    for p in ps
+        f = DeRham.fermat_hypersurface(3,6,p); f += prod((gens(parent(f)).^2))
+
+        z = DeRham.zeta_function(f,algorithm=:naive,fastevaluation=true)
+        println("p: $p, d: $d, p mod d = $(p % d), p mod 12 = $(p % 12)")
+        println(z)
+        println()
+    end
+end

src/ControlledReduction.jl

@@ -1177,7 +1177,7 @@ function computeRuvS(V,S,f,pseudoInverseMat,Ruvs,explookup,params)
     if haskey(Ruvs, V)
         return get(Ruvs, V, 0)
     else
-        println("New key: $V")
+        (4 < params.verbose) && println("New key: $V")
     end
     ev1 = gen_exp_vec(n+1,n*d-n,termorder)
     ev2 = gen_exp_vec(n+1,n*d-n+d-length(S),termorder)

src/SmoothNondegenerate.jl

@@ -18,6 +18,11 @@ function find_Ssmooth_model(f, M, S_target, params)
     vars = gens(PR)
     l = d * n - n + d - length(S_target)
 
+    #TODO: fix the anti-pattern below and remove this band-aid
+    if !check_smoothness(f)
+        error("f is not smooth")
+    end
+
     bool = true 
     f_transformed = f
     while bool

src/ZetaFunction.jl

@@ -176,9 +176,9 @@ function precision_information(f,basis,verbose=0)
     (0 < verbose) && println("Hodge numbers = $hodge_numbers")
 
     # for K3 over F3 only 
-    if (n == 3) && (p == 3) && (d == 4)
-        return (hodge_polygon, [1, 2, 2], [4, 4, 3], 6)
-    end 
+    #if (n == 3) && (p == 3) && (d == 4)
+    #    return (hodge_polygon, [1, 2, 2], [4, 4, 3], 6)
+    #end 
 
     k = sum(hodge_numbers) # dimension of H^n
     (0 < verbose) && println("There are $k basis elements in H^$n")

test/varieties/vanishing.csv

@@ -0,0 +1,3 @@
+3;11;5*x1^6*x2 + 3*x1^6*x3 + 7*x1^5*x2^2 + 9*x1^5*x2*x3 + x1^5*x3^2 + 8*x1^4*x2^3 + x1^4*x2^2*x3 + 5*x1^4*x2*x3^2 + 6*x1^3*x2^4 + 6*x1^3*x2^3*x3 + 9*x1^3*x2^2*x3^2 + 8*x1^3*x2*x3^3 + 5*x1^3*x3^4 + 5*x1^2*x2^5 + 2*x1^2*x2^4*x3 + 7*x1^2*x2^3*x3^2 + 4*x1^2*x2*x3^4 + 10*x1^2*x3^5 + 7*x1*x2^6 + 3*x1*x2^5*x3 + 5*x1*x2^4*x3^2 + 8*x1*x2^3*x3^3 + 2*x1*x2^2*x3^4 + 3*x1*x2*x3^5 + 5*x1*x3^6 + 2*x2^7 + 4*x2^6*x3 + 2*x2^5*x3^2 + 2*x2^4*x3^3 + 9*x2^3*x3^4 + 7*x2^2*x3^5 + 4*x2*x3^6 + 2*x3^7;ZZRingElem[4177248169415651,-1518999334332964,655931530734689,-200859416110144,34237400473320,-8455600419926,858292959524,-267305524607,68887149485,-21921295814,5806049601,2472191414,-1995260839,889046048,-336993701,104438702,-30635791,7347488,-1499069,168854,36051,-12374,3535,-1247,364,-326,120,-64,19,-4,1]
+3;11;10*x1^5 + 6*x1^4*x2 + 9*x1^4*x3 + 6*x1^3*x2^2 + x1^3*x2*x3 + 3*x1^3*x3^2 + 8*x1^2*x2^3 + 10*x1^2*x2^2*x3 + x1^2*x2*x3^2 + 3*x1^2*x3^3 + x1*x2^4 + 5*x1*x2^3*x3 + x1*x2^2*x3^2 + 2*x1*x2*x3^3 + 6*x1*x3^4 + 6*x2^5 + 7*x2^4*x3 + 6*x2^3*x3^2 + 6*x2^2*x3^3 + 9*x2*x3^4 + 9*x3^5;ZZRingElem[1771561, 0, 0, 31944, 5445, -517, 642, -47, 45, 24, 0, 0, 1]
+3;11;4*x1^7 + 3*x1^6*x2 + 4*x1^6*x3 + 6*x1^5*x2^2 + 10*x1^5*x2*x3 + 9*x1^5*x3^2 + x1^4*x2^3 + 9*x1^4*x2^2*x3 + x1^4*x2*x3^2 + 5*x1^4*x3^3 + 4*x1^3*x2^4 + 4*x1^3*x2^3*x3 + 10*x1^3*x2^2*x3^2 + 6*x1^3*x2*x3^3 + 4*x1^3*x3^4 + 3*x1^2*x2^5 + 6*x1^2*x2^4*x3 + 5*x1^2*x2^3*x3^2 + 8*x1^2*x2^2*x3^3 + 5*x1^2*x2*x3^4 + 7*x1^2*x3^5 + 7*x1*x2^6 + 5*x1*x2^5*x3 + 5*x1*x2^4*x3^2 + 3*x1*x2^3*x3^3 + 4*x1*x2^2*x3^4 + 10*x1*x2*x3^5 + 9*x1*x3^6 + 2*x2^7 + 3*x2^6*x3 + 4*x2^5*x3^2 + 2*x2^4*x3^3 + 5*x2^3*x3^4 + 7*x2^2*x3^5 + 9*x2*x3^6 + 6*x3^7;ZZRingElem[4177248169415651, 759499667166482, 379749833583241, 0, 30528348755377, -6951229793068, 3560501013410, -405567002852, 410107513695, -89407140548, 30683114418, -8633709854, 3882243497, -724848806, 390799937, -51056182, 35527267, -5990486, 2916787, -589694, 190518, -50468, 21045, -1892, 1510, -268, 107, 0, 11, 2, 1]