#include "math/polynomial/algebraic_numbers.h"
#include "util/rlimit.h"
#include <iostream>
namespace polynomial {
void test_algebraic_basic_operations() {
std::cout << "test_algebraic_basic_operations\n";
reslimit rl;
unsynch_mpq_manager qm;
anum_manager am(rl, qm);
scoped_anum a(am), b(am), c(am);
scoped_mpq q2(qm), q3(qm);
qm.set(q2, 2);
qm.set(q3, 3);
am.set(a, q2);
am.set(b, q3);
am.add(a, b, c);
scoped_mpq q5(qm);
qm.set(q5, 5);
VERIFY(am.eq(c, q5));
}
void test_algebraic_arithmetic() {
std::cout << "test_algebraic_arithmetic\n";
reslimit rl;
unsynch_mpq_manager qm;
anum_manager am(rl, qm);
scoped_anum a(am), b(am), sum(am), diff(am), prod(am);
scoped_mpq qa(qm), qb(qm);
qm.set(qa, 5, 2);
qm.set(qb, 3, 4);
am.set(a, qa);
am.set(b, qb);
am.add(a, b, sum); am.sub(a, b, diff); am.mul(a, b, prod);
scoped_mpq qsum(qm), qdiff(qm), qprod(qm);
qm.set(qsum, 13, 4);
qm.set(qdiff, 7, 4);
qm.set(qprod, 15, 8);
VERIFY(am.eq(sum, qsum));
VERIFY(am.eq(diff, qdiff));
VERIFY(am.eq(prod, qprod));
}
void test_algebraic_comparison() {
std::cout << "test_algebraic_comparison\n";
reslimit rl;
unsynch_mpq_manager qm;
anum_manager am(rl, qm);
scoped_anum a(am), b(am), c(am);
scoped_mpq q2(qm), q3(qm);
qm.set(q2, 2);
qm.set(q3, 3);
am.set(a, q2);
am.set(b, q3);
am.set(c, q2);
VERIFY(am.lt(a, b)); VERIFY(am.gt(b, a)); VERIFY(am.eq(a, c)); VERIFY(!am.eq(a, b)); }
void test_algebraic_degree() {
std::cout << "test_algebraic_degree\n";
reslimit rl;
unsynch_mpq_manager qm;
anum_manager am(rl, qm);
scoped_anum rational_num(am);
scoped_mpq q(qm);
qm.set(q, 5, 3);
am.set(rational_num, q);
VERIFY(am.is_rational(rational_num));
scoped_anum zero(am);
am.reset(zero);
VERIFY(am.is_zero(zero));
VERIFY(am.is_rational(zero));
}
void test_algebraic_signs() {
std::cout << "test_algebraic_signs\n";
reslimit rl;
unsynch_mpq_manager qm;
anum_manager am(rl, qm);
scoped_anum pos(am), neg(am), zero(am);
scoped_mpq qpos(qm), qneg(qm);
qm.set(qpos, 5);
qm.set(qneg, -3);
am.set(pos, qpos);
am.set(neg, qneg);
am.reset(zero);
VERIFY(am.is_pos(pos));
VERIFY(am.is_neg(neg));
VERIFY(am.is_zero(zero));
VERIFY(!am.is_pos(neg));
VERIFY(!am.is_neg(pos));
VERIFY(!am.is_zero(pos));
}
void test_algebraic_numbers() {
test_algebraic_basic_operations();
test_algebraic_arithmetic();
test_algebraic_comparison();
test_algebraic_degree();
test_algebraic_signs();
}
}
void tst_algebraic_numbers() {
polynomial::test_algebraic_numbers();
}