Letting the magma function also compute the trace for p = 1

assaferan · Jun 23, 2023 · f7a1bf4 · f7a1bf4
1 parent 8686eab
commit f7a1bf4
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/magma_functions.m b/magma_functions.m
@@ -201,6 +201,7 @@ function TrivialContribution(N, k, p)
 
 // At the moment only works for Hecke operators at primes
 function TraceFormulaGamma0ALNew(p, N, k)
+    if (p eq 1) then return TraceFormulaGamma0ALTrivialNew(N, k); end if;
     assert IsPrime(p);
     ms := [d : d in Divisors(N) | (N mod d^2 eq 0) and (d mod p ne 0)];
     trace := &+[Integers() | MoebiusMu(m)*(TraceFormulaGamma0AL(p, N div m^2, k) - TrivialContribution(N div m^2, k, p)) : m in ms];