QEtaModularCategory¶
qetamod.spad line 162 [edit on github]
QEtaModularCategory provides functions to check modularity conditions.
- candidateLevelsCoEtaQuotient: (QEtaSpecification, PositiveInteger, NonNegativeInteger) -> List PositiveInteger
candidateLevelsCoEtaQuotient(rspec,m,t)returns an ascendingly sorted list of thennup to 24*m*level(rspec) such that conditionCoEtaQuotient?(nn,rspec,m,t) istrue.
- conditionCoEtaQuotient?: (PositiveInteger, QEtaSpecification, PositiveInteger, NonNegativeInteger) -> Boolean
conditionCoEtaQuotient?(nn,rspec,m,t)returnstrueiff a cofactor eta-quotient can be found such that the product of the cofactor with the orbit product corresponding to the dissection (m,t) gives a modular function in levelnn.
- cuspMatrices: PositiveInteger -> List Matrix Integer
cuspMatrices(nn)returns [cuspToMatrix(nn,c) forcin cusps(nn)].
- cusps: PositiveInteger -> List Cusp
cusps(nn)returns representatives for all the (inequivalent) cusps for the congruence subgroup that the respective package represents. They are sorted by their size as rational numbers with infinity being the biggest cusp. Note that cusps()=[cusp(x(1,1),x(2,1)) forxin doubleCosetRepresentatives()] where doupleCosetRepresentatives() is from the respective congruence subgroup that is connected to the package.
- cuspToMatrix: (PositiveInteger, Cusp) -> Matrix Integer
For cusp=(a:c), cuspToMatrix(
nn,cusp) returns a matrix gamma=[[a,b],[c,d]] corresponding to the cusp (a:c) of this domain. We assume that a/c is a normalized cusp, i.e. cusp=normalizeCusp(cusp).
- etaCofactorSpace: (PositiveInteger, List List Integer, QEtaSpecification, PositiveInteger, NonNegativeInteger) -> Record(level: PositiveInteger, indices: List List Integer, particular: Union(Vector Integer, failed), basis: List Vector Integer)
etaCofactorSpace(nn,idxs,rspec,m,t)returns a vectorvand the basis of a space such that #v=#idxs and modular?(nn,members(s),rspec,m,t) istruefor anys=v+ reduce(_+, [z.i* basis.iforiin 1..#basis]) and any sufficiently long listzof integers. The function fails, if there is no such solution. The indices part of the result is equal toidxsand corresponds to the entries of the solution vectors.
- etaCofactorSpaceSystem: (PositiveInteger, List List Integer, QEtaSpecification, PositiveInteger, NonNegativeInteger) -> Record(qmat: Matrix Fraction Integer, zrels: Vector Integer, qrhs: Vector Fraction Integer, qsgn: List Integer)
etaCofactorSpaceSystem(nn,idxs,rspec,m,t)returns a matrix mat and a vectorvsuch that for the integer solutionssof the equation mat*s=v it holds modular?(sspec,rspec,m,t) where sspec results from the list of the first #idxs entries ofs. Withidxsone can specify what indexes should be allowed in the solution space.
- etaQuotientMonoidInfinitySystem: (PositiveInteger, List List Integer) -> Record(qmat: Matrix Fraction Integer, zrels: Vector Integer, qrhs: Vector Fraction Integer, qsgn: List Integer)
etaQuotientMonoidInfinitySystem(nn,idxs)returns etaQuotientMonoidSystem(nn,idxs,[infinity()$Cusp]).
- etaQuotientMonoidNoPolesSystem: (PositiveInteger, List List Integer, List Cusp) -> Record(qmat: Matrix Fraction Integer, zrels: Vector Integer, qrhs: Vector Fraction Integer, qsgn: List Integer)
etaQuotientMonoidNoPolesSystem(nn,idxs,nopolesat)returns etaQuotientMonoid(nn,idxs,[],nopolesat).
- etaQuotientMonoidSystem: (PositiveInteger, List List Integer, List Cusp) -> Record(qmat: Matrix Fraction Integer, zrels: Vector Integer, qrhs: Vector Fraction Integer, qsgn: List Integer)
etaQuotientMonoidSystem(nn,idxs,polesat)returns etaQuotientMonoidNoPolesSystem(nn,idxs,nopolesat) = etaQuotientMonoidSystem(nn,idxs,[],nopolesat) for nopolesat:=[cforcin cuspsnn| not member?(c,polesat)].
- etaQuotientMonoidSystem: (PositiveInteger, List List Integer, List Cusp, List Cusp) -> Record(qmat: Matrix Fraction Integer, zrels: Vector Integer, qrhs: Vector Fraction Integer, qsgn: List Integer)
etaQuotientMonoidSystem(nn,idxs,polesat,nopolesat)returns a record that can be given to zsolve$X4ti2in order to find the generators of the monoid of (generalized) eta-quotients that are modularwrt. QMOD and have no poles at the cusps given bynopolesat. Furthermore, setup the system in such a way that at all cusps inpolesatthe expansion order of a solution of the system will be non-positive. The indices of the eta-functions that may appear is given byidxs. Note that the union of poleat andnopolesatcan be a subset of all cusps given by cusps(nn).
- etaQuotientMonoidSystem: (PositiveInteger, List List Integer, List Cusp, List Integer, List Integer, List Cusp) -> Record(qmat: Matrix Fraction Integer, zrels: Vector Integer, qrhs: Vector Fraction Integer, qsgn: List Integer)
etaQuotientMonoidSystem(nn,idxs,polesat,rels,ords,nopolesat)returns a record that can be given to zsolve$X4ti2in order to find the generators of the monoid of (generalized) eta-quotients that are modularwrt. QMOD and have no poles at the cusps given bynopolesatand for all cusps inpolesatthe expansion order fulfills the relation at the respective cusp. #polesat>=#rels and #polesat>=#ords. Missing values are assumed to be 0. Entries inrelscan be 1,0,-1and stand for>=, =, and<=, respectively. The indices of the eta-functions that may appear is given byidxs.
- genus: PositiveInteger -> NonNegativeInteger
genus(nn)returns the genus of the congruence subgroup connected to the package (eitherGamma_0(nn) orGamma_1(nn).
- matrixEtaOrder: (PositiveInteger, List Cusp, List List Integer) -> Matrix Fraction Integer
matrixEtaOrder(nn,spitzen,idxs)returns the matrix given by [matrixEtaOrderRow(nn,cusp,idxs) for cusp inspitzen] with duplicate rows removed.
- matrixEtaOrderElement: (PositiveInteger, Cusp, List Integer) -> Fraction Integer
z:=matrixEtaOrderElement(
nn,cusp,idx). The parameter corresponds to the index of a (generalized) eta-quotient, i.e. it can be of the form idx=[delta,g] or of the form idx=[delta] withgin the range 0..delta. Thenzis the coefficient ofr_{delta,g} in cite[Thm.~4]{Robins_GeneralizedDedekindEtaProducts_1994}, except for the casenn=4,gcd(denom(cusp),nn)=2(orGamma1(4)-equivalent) where we return half of this coefficient. If idx=[delta] then z=1/2*matrixEtaOrderElement(nn,cusp,[delta,0])/2. and thus corresponds to the coefficient of $r_delta$ in eqref{eq:order-rbar-non-adjusted} in Theorem~ref{thm:matrixEtaOrderRobins}. Note that Robins refers to an expansion in the uniformizing variableq^(gcd(c,nn)/nn) (where cusp=a/c). That is equal toq^(1/width(c)$CongruenceSubgroupGamma1(nn)) except for the casenn=4andgcd(nn,c)=2. In that case the width of the cusp 1/2 is 1 and not 2=4/gcd(4.2). For this case, we deviate from Robins’ coefficients and adapt to an expansion inqrather than an expansion inq^(1/2).
- matrixEtaOrderFull: (PositiveInteger, List Cusp, List List Integer) -> Matrix Fraction Integer
matrixEtaOrderFull(nn,spitzen,idxs)returns a matrix given by [matrixEtaOrderRow(nn,cusp,idxs) for cusp inspitzen].
- matrixEtaOrderRow: (PositiveInteger, Cusp, List List Integer) -> List Fraction Integer
matrixEtaOrderRow(nn, cusp, idxs)returns [matrixEtaOrderElement(nn,cusp,idx) for idx in idxs].
- matrixModular: (PositiveInteger, List List Integer) -> Matrix Fraction Integer
matrixModular(nn, idxs)returns a matrix matmat*v=0for an (extended) exponent vector where the exponents correspond to idxs, i.e. #v=#idxs+k wherekcorresponds to additional (slack) variables due to the modularity conditions. The rows corresponding to slack variables have been normalized so that the entry of the slack variable is 1. The first row corresponds to 1/2 of the sum of exponents of pure eta-functions, so that when multiplied by the exponent vectors corresponding to idxs it gives the weight of the eta-quotient. The second row (when multiplied by the exponent vectors corresponding to idxs) gives the expansion order at infinity. The third row (when multiplied by the exponent vectors corresponding to idxs) gives the expansion order at the cusp 0.
- minimalLevelCoEtaQuotient: (QEtaSpecification, PositiveInteger, NonNegativeInteger) -> PositiveInteger
minimalLevelCoEtaQuotient(rspec, m, t)returns the smallestnnof candidateLevalsCoEtaQuotient(rspec,m,t).
- modular?: (QEtaSpecification, QEtaSpecification, PositiveInteger, NonNegativeInteger) -> Boolean
modular?(sspec,rspec,m,t)returnstrueiff (sspec,rspec,m,t) specifies a modular function in level(sspec).
- modular?: QEtaSpecification -> Boolean
modular?(rspec)returnstrueiff spec specifies a modular function in level(rspec).
- modularOrbit: (QEtaSpecification, PositiveInteger, NonNegativeInteger) -> List NonNegativeInteger
modularOrbit(rspec, m, t)(for $Gamma_0(N) whereNis level(rspec)) computes the elements of $modularOrbit{r,m,t)$ as defined in qeta.tex, cite[Def.~42]{Radu_RamanujanKolberg_2015} and cite[Lemma 4.35]{Radu_PhD_2010} wereris pureExponents(rspec) andrspecis expected to specify a pure eta-quotient. For $Gamma_1(N)$, modularOrbit(rspec,m,t) returns [t].
- numberOfGapsForSamba: PositiveInteger -> Integer
numberOfGapsForSamba(nn)returns the number of gaps (counted from the elements in the algorithm SAMBA that is sufficient to ensure a module basis, i.e. it is a criterion to avoid needless computation. ForCongruenceSubgroupGamma0this is equal to genus(nn). If no theorem for another group is known it will return-1, with the meaning that SAMBA must really check all critical element for zero reduction.
- rightCosetRepresentatives: (PositiveInteger, PositiveInteger) -> List Matrix Integer
rightCosetRepresentatives(nn, mm)returns representatives of $G(nn) diagdownG(mm)$ for $G$ being the group constructor represented by this domain, i.e. either $Gamma_0$ or $Gamma_1$.
- weaklyModular?: QEtaSpecification -> Boolean
weaklyModular?(spec)returnstrueiffspecspecifies a weakly modular function of some integer weightk.
- width: (PositiveInteger, Cusp) -> PositiveInteger
width(nn, cusp)returns the width of the cusp=(a:c) forG(nn) whereGis the congruence subgroup ofSL_2(ZZ) connected to this package.
- width: (PositiveInteger, Matrix Integer) -> PositiveInteger
width(nn, gamma)returns the width of the cusp=(a:c) forG(nn) whereGis the congruence subgroup ofSL_2(ZZ) connected to this package and gamma = matrix [[a,b],[c,d]].