forked from JohnCremona/eclib
/
tperiods.cc
59 lines (53 loc) · 1.86 KB
/
tperiods.cc
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
// tperiods.cc -- input a curve, find its periods
//////////////////////////////////////////////////////////////////////////
//
// 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
//
//////////////////////////////////////////////////////////////////////////
//
//#define TEST
#include <eclib/compproc.h>
#include <eclib/cperiods.h>
#include <eclib/getcurve.h>
int main(){
#ifdef MPFP
set_precision("Enter precision in bits");
long original_output_precision = RR::OutputPrecision();
RR::SetOutputPrecision(original_output_precision-3);
#endif
initprimes("PRIMES",0);
int verb=1;
bigint v;
vector<bigrational> ai(5);
while (getcurve(ai,verb))
{
cout << "Input curve = ";
cout <<"["<<ai[0]<<","<<ai[1]<<","<<ai[2]<<","<<ai[3]<<","<<ai[4]<<"]:" << endl;
Curvedata CDin(ai,v); // v holds scaling factor
Curvedata CD(CDin,1); // 1 for minimise
cout << "Minimal model = "<<(Curve)CD<<endl;
Cperiods cp(CD);
cout << "Periods: " << cp << endl;
Curve EE = cp.trans_to_curve();
cout << "Curve from periods: " << EE << endl;
}
#ifdef MPFP
RR::SetOutputPrecision(original_output_precision);
#endif
}