Ariel Pacetti on Wed, 28 Mar 2012 16:52:39 +0200


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

forqfvec


Dear pari developers,

sometime ago Bill wrote a routime called forqfvec which is a mix of qfminim and somehow forvec (see http://pari.math.u-bordeaux.fr/archives/pari-dev-1005/msg00003.html).
The advantage of such routine is that instead of giving the whole output 
of qfminim (which might blow up the memory) it lets you manage each vector 
at a time and then discard it. The first application of it is the routine 
qfrep writen by Karim.
In such routine, one takes each vector of qfrep and compute its norm, and 
then puts a 1 in the right place of a vector. This gives the q-expansion 
of modular forms (of integral weight or half integral weight) which is 
quite useful. I was asked by Bill to make some lobby on this routine, so 
let me just mention some other applications of it so as to consider 
including it in the future releases.
1) Say you want to compute theta series of quadratic forms. Then the 
routine qfrep does the job. But if you want to compute the form with a 
character for example, then you need the forqfvec routine. I used such 
routine in the paper "Computing central values of twisted L-series, the 
case of composite levels" with Gonzalo TornarĂ­a. At the time, Gonzalo 
implemented the particular routine in C.
2) Say you want to compute modular forms (from quaternion algebras) of 
weight greater than 2. Then what you do is to add an harmonic polynomial 
to the theta function, and evaluate it in the elements for the theta 
series. Again, forqfvec will do the job (you can check Eichler's paper on 
the basis problem to see how this works).
3) Say you want to compute Hilbert modular forms (as in Lassina Dembele's 
works). Then you have forms over a totaly real field and need to evaluate 
the elements at some representation (if the class number is greater than 
one). Again, forqfvec does the job. Furthermore, I wrote a script for a 
particular case of a HMF over Q[\sqrt{5}] to compute many Fourier 
coefficients using forqfvec (and was the first time I asked Bill such 
routine).
There should be many other applications relating modular forms that I am 
missing, so I believe having it in the official release might be quite 
useful.
Cheers,

Ariel