summaryrefslogtreecommitdiff log msg author committer range
path: root/main/mapc_poly_main.cc
blob: 9aa0f65d4afe8a53a0e5d5a62a75aacb88d98c3b (plain)
 ```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 #include #include #include #include #include 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; } ``````