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
|
#include <cassert>
#include <cstdlib>
#include <fstream>
#include <iostream>
#include <kratpoly.h>
#include <kfloatpoly.h>
using namespace std;
//int main(const int argc, const char* argv[])
int main()
{
// ofstream of(argv[1]);
unsigned long i, j;
i = 2;
j = 1;
if (i > j)
cerr << " i > j " << endl << flush;
// long xp[] = { 1, 1 };
// bigrational xc[] = { 1, 1, 1, 1 };
long xp[] = { 1, 2 };
bigrational xc[] = { 1, 1, 1, 0, 1, 1 };
K_RATPOLY X(2, xp, 6, xc);
cerr << " X = " << endl << flush;
cerr << X << endl << flush;
// long yp[] = { 3 };
// bigrational yc[] = { bigrational(1, 2), 0, 0, bigrational(1, 2) };
// bigrational yc[] = { bigrational(1, 2), 0, 0, 0 };
long yp[] = { 4 };
bigrational yc[] = { bigrational(1, 2), - bigrational(1, 2), - bigrational(3, 2), bigrational(5, 2), - 1 };
K_RATPOLY Y(1, yp, 5, yc);
cerr << " Y = " << endl << flush;
cerr << Y << endl << flush;
long zp[] = { 1 };
bigrational zc[] = { bigrational(1, 2), bigrational(1, 3) };
// long zp[] = { 0 };
// bigrational zc[] = { bigrational(1, 2) };
K_RATPOLY Z(1, zp, 2, zc);
cerr << " Z = " << endl << flush;
cerr << Z << endl << flush;
// K_RATPOLY W = X.subst_expr(0, Y, Z);
K_RATPOLY W = X.subst_expr(0, Z);
cerr << " W = " << endl << flush;
cerr << W << endl << flush;
// K_RATPOLY Q, R;
// Q = div(Y, Z, R);
// cerr << " Q = " << endl << flush;
// cerr << Q << endl << flush;
// cerr << " R = " << endl << flush;
// cerr << R << endl << flush;
// K_RATPOLY G = gcd(Z, Y);
// cerr << " G = " << endl << flush;
// cerr << G << endl << flush;
bigrational v = - 1;
bigrational V[] = { - 1 };
cerr << " Y(" << v << ") = " << Y.evaluate(v) << endl << flush;
cerr << " Y(" << V[0] << ") = " << Y.evaluate(V) << endl << flush;
cerr << " sgn(Y(" << V[0] << ")) = " << Y.sgn_at(V) << endl << flush;
Y.set_Sturm_seq();
cerr << Y.num_Sturm_seq_perm(- 1) << endl << flush;
cerr << Y.num_Sturm_seq_perm(0) << endl << flush;
cerr << Y.num_Sturm_seq_perm(1) << endl << flush;
// long fp[] = { 3 };
// double fc[] = { 1.0, 0.0, - 1.0, 0.0 };
// K_FLOATPOLY F(1, fp, fc);
// cerr << " F = " << endl << flush;
// cerr << F << endl << flush;
//
// unsigned long l_FR;
// double* FR;
// l_FR = F.gen_fp_roots(-2.0, 2.0, &FR);
// cerr << " F has " << l_FR << " real roots " << endl << flush;
//
// for (i = 0; i < l_FR; i++)
// cerr << FR[i] << endl << flush;
//
// if (l_FR)
// delete [] FR;
return 0;
}
|
