Newform basis and dimension don't agree
ahorawa opened this issue · 4 comments
ahorawa commented
I tried to run:
F:=QuadraticField(10);
prec:=20;
M:=GradedRingOfHMFs(F, prec);
ZF:=Integers(F);
N:= 23*ZF;
H := HeckeCharacterGroup(N, [1,2]);
chi := H.2^11; // (H.1^11); // aka 11 mod 22
M1chi := HMFSpace(M, N, [1,1], chi);
SetVerbose("HilbertModularForms", 3);
Space := Weight1CuspBasis(M1chi : prove := false);
but got:
Computing basis of cusp forms of weight [2,2] and level
Principal Prime Ideal of ZF
Generator:
23
Decomposing HilbertCuspForms for N = 1 and weight = [ 2, 2 ]
Decomposing HilbertCuspForms for N = 23 and weight = [ 2, 2 ]
Weight1CuspBasis(
Mk: Space of Hilbert modular forms over Quadratic Field with def...
)
CuspFormBasis(
Mk: Space of Hilbert modular forms over Quadratic Field with def...
)
In file "/Users/alexh/Documents/Michigan/mathematics/Computations/HMS/hilbertmodularforms/ModFrmHilD/Basis.m", line 68, column 7:
>> assert CuspDimension(Mk) eq #Mk`CuspFormBasis;
^
Runtime error in assert: Assertion failed
The problem seems to come up when we compute the basis of weight [2,2] and level N = 23*ZF forms in this case. Is it the trace formula correction problem when weight is 2?
edgarcosta commented
Here is a minimal example:
F:=QuadraticField(10);
prec:=5;
M:=GradedRingOfHMFs(F, prec);
ZF:=Integers(F);
N:= 23*ZF;
SetVerbose("HilbertModularForms", 3);
M2 := HMFSpace(M, N, [2,2]);
C2 := CuspFormBasis(M2);
for me, this doesn't even finish.
@ahorawa, you have not yet pushed your code with the flag prove
edgarcosta commented
After 10 mins my minimal example finished, trying now yours with prec = 5
abhijit-mudigonda commented
I suspect this is fixed but I haven't checked. If someone can, please do!
ahorawa commented
I computed it and it works.