1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
|
// FILE moreap.cc: computes more ap for given level(s)
//////////////////////////////////////////////////////////////////////////
//
// Copyright 1990-2012 John Cremona
//
// This file is part of the eclib package.
//
// eclib is free software; you can redistribute it and/or modify it
// under the terms of the GNU General Public License as published by the
// Free Software Foundation; either version 2 of the License, or (at your
// option) any later version.
//
// eclib is distributed in the hope that it will be useful, but WITHOUT
// ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
// FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
// for more details.
//
// You should have received a copy of the GNU General Public License
// along with eclib; if not, write to the Free Software Foundation,
// Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA
//
//////////////////////////////////////////////////////////////////////////
#include <fstream>
#include <eclib/moddata.h>
#include <eclib/symb.h>
#include <eclib/cusp.h>
#include <eclib/homspace.h>
#include <eclib/oldforms.h>
#include <eclib/cperiods.h> //from qcurves, for computing conductors
#include <eclib/newforms.h>
//#define AUTOLOOP
int main(void)
{
int n=1;
int lastn, stopp, output, showeigs, showforms, findcurves;
int nnf, nap;
cout << "Program moreap\n";
cout << "---------------\n\n";
cout << "For each N, assumes that the file newforms/xN exists, and computes more\n";
cout << "Hecke eigenvalues.\n";
cout << "Output new eigs to file (1/0)? "; cin>>output;
cout << "Output new eigs to screen (1/0)? "; cin>>showeigs;
cout << "Display newforms (1/0)? "; cin>>showforms;
cout << "Attempt curve construction (1/0)? "; cin>>findcurves;
#ifdef MPFP
if(findcurves) set_precision(60);
#endif
#ifdef AUTOLOOP
cout << "How many primes for Hecke eigenvalues? ";
cin >> stopp; cout << endl;
cout<<"Enter first and last N: ";cin>>n>>lastn; n--;
while (n<lastn) { n++;
#else
while (cout<<"Enter level: ", cin>>n, n>0) {
cout << "How many primes for Hecke eigenvalues? ";
cin >> stopp; cout << endl;
#endif
cout << "\n>>>Level " << n << "\t";
newforms nf(n,showforms);
nf.createfromdata(1,0,0);
if (showforms) nf.newforms::display();
nnf = nf.n1ds;
if(nnf==0)
{
cout<<"No newforms."<<endl;
continue;
}
nap = nf.nflist[0].aplist.size();
if(nap>=stopp)
{
cout<<"Already have "<<nap<<" eigenvalues on file, no need to compute more."<<endl;
}
else
{
nf.makebases(1);
cout << "About to start computing ap..."<<flush;
nf.make_projcoord();
nf.find_jlist();
nf.addap(stopp);
cout << "...done."<<endl;
if(output) nf.output_to_file();
}
if(findcurves)
{
int inf, nnf = nf.n1ds;
// Now we compute the curves
cout<<"Computing "<<nnf<<" curves...\n";
vector<int> forms;
for(inf=0; inf<nnf; inf++) forms.push_back(inf);
vector<int> failures = nf.showcurves(forms,0,"no");
if(failures.size()>0)
{
cout<<"No curve found for "<<failures.size()<<" forms: "<<failures<<endl;
}
else
{
cout<<"All curves found OK"<<endl;
}
}
} // end of while()
} // end of main()
|
