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
105
106
107
108
109
110
111
112
|
// getcurve.cc: implementation of function getcurve() for curve input
//////////////////////////////////////////////////////////////////////////
//
// 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 <eclib/getcurve.h>
int getcurve(Curvedata& CD, int verb)
{
Curve C0;
if(verb) cout << "Enter curve: ";
cin>>ws; if(cin.eof()) return 0; // quit if EOF reached
cin >> C0;
if (verb) cout << endl;
if(C0.isnull()) return 0; // quit if null curve entered
CD = Curvedata(C0,0); // DON'T change coords
if(CD.isnull()) // input curve was singular, non-null
{
cout<<C0<<" is singular"<<endl;
return 0;
}
return 1;
}
//#define DEBUG_Q_INPUT
int getcurve(vector<bigrational>& ai, int verb)
{
// read the coefficients, either as "a1 a2 a3 a4 a6" or as
// "[a1,a2,a3,a4,a6]" using code essentially the same as in
// ../qcurves/curve.cc
ai.resize(5);
if(verb) cerr << "Enter curve: ";
cin>>ws; if(cin.eof()) return 0; // quit if EOF reached
char c=0;
cin.get(c);
#ifdef DEBUG_Q_INPUT
cout<<"First char read = "<<c<<"\n";
#endif
if(c=='[')
{
#ifdef DEBUG_Q_INPUT
cout<<"Reading [a1,a2,a3,a4,a6]...\n";
#endif
cin >> ai[0] >> c;
if(c!=',')
{
cerr << "syntax error on curve input" << endl;
return 0;
}
cin >> ai[1] >> c;
if(c!=',')
{
cerr << "syntax error on curve input" << endl;
return 0;
}
cin >> ai[2] >> c;
if(c!=',')
{
cerr << "syntax error on curve input" << endl;
return 0;
}
cin >> ai[3] >> c;
if(c!=',')
{
cerr << "syntax error on curve input" << endl;
return 0;
}
cin >> ai[4] >> c;
if(c!=']')
{
cerr << "syntax error on curve input" << endl;
return 0;
}
#ifdef DEBUG_Q_INPUT
cout<<ai<<endl;
#endif
}
else
{
#ifdef DEBUG_Q_INPUT
cout<<"Reading a1 a2 a3 a4 a6 ...\n";
#endif
cin.unget();
cin >> ai[0] >> ai[1] >> ai[2] >> ai[3] >> ai[4];
#ifdef DEBUG_Q_INPUT
cout<<ai<<endl;
#endif
}
// test for null curve input
if((num(ai[0])==0)&&(num(ai[1])==0)&&(num(ai[2])==0)&&(num(ai[3])==0)&&(num(ai[4])==0))
return 0; // quit if null curve entered
return 1;
}
|
