[ceph.git] / ceph / src / boost / libs / polygon / test / voronoi_ctypes_test.cpp

// Boost.Polygon library voronoi_ctypes_test.cpp file

//          Copyright Andrii Sydorchuk 2010-2012.
// Distributed under the Boost Software License, Version 1.0.
//    (See accompanying file LICENSE_1_0.txt or copy at
//          http://www.boost.org/LICENSE_1_0.txt)

// See http://www.boost.org for updates, documentation, and revision history.

#include <boost/core/lightweight_test.hpp>
#include <boost/polygon/detail/voronoi_ctypes.hpp>
#include <boost/random/mersenne_twister.hpp>
#include <vector>
#include <ctime>

using namespace boost::polygon::detail;

type_converter_fpt to_fpt;

void ulp_comparison_test1()
{
  ulp_comparison<double> ulp_cmp;
  uint64 a = 22;
  uint64 b = 27;
  fpt64 da, db;
  std::memcpy(&da, &a, sizeof(uint64));
  std::memcpy(&db, &b, sizeof(uint64));
  BOOST_TEST_EQ(ulp_cmp(da, db, 1), ulp_cmp.LESS);
  BOOST_TEST_EQ(ulp_cmp(db, da, 1), ulp_cmp.MORE);
  BOOST_TEST_EQ(ulp_cmp(da, db, 4), ulp_cmp.LESS);
  BOOST_TEST_EQ(ulp_cmp(da, db, 5), ulp_cmp.EQUAL);
  BOOST_TEST_EQ(ulp_cmp(da, db, 6), ulp_cmp.EQUAL);
}

void ulp_comparison_test2()
{
  ulp_comparison<fpt64> ulp_cmp;
  uint64 a = 0ULL;
  uint64 b = 0x8000000000000002ULL;
  fpt64 da, db;
  std::memcpy(&da, &a, sizeof(uint64));
  std::memcpy(&db, &b, sizeof(uint64));
  BOOST_TEST_EQ(ulp_cmp(da, db, 1), ulp_cmp.MORE);
  BOOST_TEST_EQ(ulp_cmp(db, da, 1), ulp_cmp.LESS);
  BOOST_TEST_EQ(ulp_cmp(da, db, 2), ulp_cmp.EQUAL);
  BOOST_TEST_EQ(ulp_cmp(da, db, 3), ulp_cmp.EQUAL);
}

void extended_exponent_fpt_test1()
{
  boost::mt19937_64 gen(static_cast<uint32>(time(NULL)));
  fpt64 b = 0.0;
  efpt64 eeb(b);
  for (int i = 0; i < 1000; ++i) {
    fpt64 a = to_fpt(static_cast<int64>(gen()));
    efpt64 eea(a);
    efpt64 neg = -eea;
    efpt64 sum = eea + eeb;
    efpt64 dif = eea - eeb;
    efpt64 mul = eea * eeb;
    BOOST_TEST_EQ(to_fpt(neg), -a);
    BOOST_TEST_EQ(to_fpt(sum), a + b);
    BOOST_TEST_EQ(to_fpt(dif), a - b);
    BOOST_TEST_EQ(to_fpt(mul), a * b);
  }
}

void extended_exponent_fpt_test2()
{
  boost::mt19937_64 gen(static_cast<uint32>(time(NULL)));
  fpt64 a = 0.0;
  efpt64 eea(a);
  for (int i = 0; i < 1000; ++i) {
    fpt64 b = to_fpt(static_cast<int64>(gen()));
    if (b == 0.0)
      continue;
    efpt64 eeb(b);
    efpt64 neg = -eea;
    efpt64 sum = eea + eeb;
    efpt64 dif = eea - eeb;
    efpt64 mul = eea * eeb;
    efpt64 div = eea / eeb;
    BOOST_TEST_EQ(to_fpt(neg), -a);
    BOOST_TEST_EQ(to_fpt(sum), a + b);
    BOOST_TEST_EQ(to_fpt(dif), a - b);
    BOOST_TEST_EQ(to_fpt(mul), a * b);
    BOOST_TEST_EQ(to_fpt(div), a / b);
  }
}

void extended_exponent_fpt_test3()
{
  boost::mt19937_64 gen(static_cast<uint32>(time(NULL)));
  for (int i = 0; i < 1000; ++i) {
    fpt64 a = to_fpt(static_cast<int64>(gen()));
    fpt64 b = to_fpt(static_cast<int64>(gen()));
    if (b == 0.0)
      continue;
    efpt64 eea(a);
    efpt64 eeb(b);
    efpt64 neg = -eea;
    efpt64 sum = eea + eeb;
    efpt64 dif = eea - eeb;
    efpt64 mul = eea * eeb;
    efpt64 div = eea / eeb;
    BOOST_TEST_EQ(to_fpt(neg), -a);
    BOOST_TEST_EQ(to_fpt(sum), a + b);
    BOOST_TEST_EQ(to_fpt(dif), a - b);
    BOOST_TEST_EQ(to_fpt(mul), a * b);
    BOOST_TEST_EQ(to_fpt(div), a / b);
  }
}

void extended_exponent_fpt_test4()
{
  for (int exp = 0; exp < 64; ++exp)
  for (int i = 1; i < 100; ++i) {
    fpt64 a = i;
    fpt64 b = to_fpt(1LL << exp);
    efpt64 eea(a);
    efpt64 eeb(b);
    efpt64 neg = -eea;
    efpt64 sum = eea + eeb;
    efpt64 dif = eea - eeb;
    efpt64 mul = eea * eeb;
    efpt64 div = eea / eeb;
    BOOST_TEST_EQ(to_fpt(neg), -a);
    BOOST_TEST_EQ(to_fpt(sum), a + b);
    BOOST_TEST_EQ(to_fpt(dif), a - b);
    BOOST_TEST_EQ(to_fpt(mul), a * b);
    BOOST_TEST_EQ(to_fpt(div), a / b);
  }
}

void extended_exponent_fpt_test5()
{
  for (int i = 0; i < 100; ++i) {
    efpt64 a(to_fpt(i * i));
    efpt64 b = a.sqrt();
    BOOST_TEST_EQ(to_fpt(b), to_fpt(i));
  }
}

void extended_exponent_fpt_test6()
{
  for (int i = -10; i <= 10; ++i) {
    efpt64 a(to_fpt(i));
    BOOST_TEST_EQ(is_pos(a), i > 0);
    BOOST_TEST_EQ(is_neg(a), i < 0);
    BOOST_TEST_EQ(is_zero(a), !i);
  }
}

void extended_int_test1()
{
  typedef extended_int<1> eint32;
  eint32 e1(0), e2(32), e3(-32);
  BOOST_TEST_EQ(e1.count(), 0);
  BOOST_TEST_EQ(e1.size(), 0U);
  BOOST_TEST_EQ(e2.count(), 1);
  BOOST_TEST_EQ(e2.chunks()[0], 32U);
  BOOST_TEST_EQ(e2.size(), 1U);
  BOOST_TEST_EQ(e3.count(), -1);
  BOOST_TEST_EQ(e3.chunks()[0], 32U);
  BOOST_TEST_EQ(e3.size(), 1U);
}

void extended_int_test2()
{
  typedef extended_int<2> eint64;
  int64 val64 = 0x7fffffffffffffffLL;
  eint64 e1(0), e2(32), e3(-32), e4(val64), e5(-val64);
  BOOST_TEST_EQ(e1.count(), 0);
  BOOST_TEST_EQ(e2.count(), 1);
  BOOST_TEST_EQ(e2.chunks()[0], 32U);
  BOOST_TEST_EQ(e3.count(), -1);
  BOOST_TEST_EQ(e3.chunks()[0], 32U);
  BOOST_TEST_EQ(e4.count(), 2);
  BOOST_TEST_EQ(e4.chunks()[0], 0xffffffff);
  BOOST_TEST_EQ(e4.chunks()[1], val64 >> 32);
  BOOST_TEST_EQ(e5.count(), -2);
  BOOST_TEST_EQ(e5.chunks()[0], 0xffffffff);
  BOOST_TEST_EQ(e5.chunks()[1], val64 >> 32);
}

void extended_int_test3()
{
  typedef extended_int<2> eint64;
  std::vector<uint32> chunks;
  chunks.push_back(1);
  chunks.push_back(2);
  eint64 e1(chunks, true), e2(chunks, false);
  BOOST_TEST_EQ(e1.count(), 2);
  BOOST_TEST_EQ(e1.chunks()[0], 2U);
  BOOST_TEST_EQ(e1.chunks()[1], 1U);
  BOOST_TEST_EQ(e2.count(), -2);
  BOOST_TEST_EQ(e2.chunks()[0], 2U);
  BOOST_TEST_EQ(e2.chunks()[1], 1U);
}

void extended_int_test4()
{
  typedef extended_int<2> eint64;
  std::vector<uint32> chunks;
  chunks.push_back(1);
  chunks.push_back(2);
  eint64 e1(chunks, true), e2(chunks, false);
  BOOST_TEST_EQ(e1 == e2, false);
  BOOST_TEST_EQ(e1 == -e2, true);
  BOOST_TEST_EQ(e1 != e2, true);
  BOOST_TEST_EQ(e1 != -e2, false);
  BOOST_TEST_EQ(e1 < e2, false);
  BOOST_TEST_EQ(e1 < -e2, false);
  BOOST_TEST_EQ(e1 <= e2, false);
  BOOST_TEST_EQ(e1 <= -e2, true);
  BOOST_TEST_EQ(e1 > e2, true);
  BOOST_TEST_EQ(e1 > -e2, false);
  BOOST_TEST_EQ(e1 >= e2, true);
  BOOST_TEST_EQ(e1 >= -e2, true);
}

void extended_int_test5()
{
  typedef extended_int<2> eint64;
  boost::mt19937_64 gen(static_cast<uint32>(time(NULL)));
  for (int i = 0; i < 1000; ++i) {
    int64 i1 = static_cast<int64>(gen());
    int64 i2 = static_cast<int64>(gen());
    eint64 e1(i1), e2(i2);
    BOOST_TEST_EQ(e1 == e2, i1 == i2);
    BOOST_TEST_EQ(e1 != e2, i1 != i2);
    BOOST_TEST_EQ(e1 > e2, i1 > i2);
    BOOST_TEST_EQ(e1 >= e2, i1 >= i2);
    BOOST_TEST_EQ(e1 < e2, i1 < i2);
    BOOST_TEST_EQ(e1 <= e2, i1 <= i2);
  }
}

void extended_int_test6()
{
  typedef extended_int<1> eint32;
  eint32 e1(32);
  eint32 e2 = -e1;
  BOOST_TEST_EQ(e2.count(), -1);
  BOOST_TEST_EQ(e2.size(), 1U);
  BOOST_TEST_EQ(e2.chunks()[0], 32U);
}

void extended_int_test7()
{
  typedef extended_int<2> eint64;
  boost::mt19937_64 gen(static_cast<uint32>(time(NULL)));
  for (int i = 0; i < 1000; ++i) {
    int64 i1 = static_cast<int64>(gen()) >> 2;
    int64 i2 = static_cast<int64>(gen()) >> 2;
    eint64 e1(i1), e2(i2), e3(i1 + i2), e4(i1 - i2);
    BOOST_TEST(e1 + e2 == e3);
    BOOST_TEST(e1 - e2 == e4);
  }
}

void extended_int_test8()
{
  typedef extended_int<2> eint64;
  boost::mt19937 gen(static_cast<uint32>(time(NULL)));
  for (int i = 0; i < 1000; ++i) {
    int64 i1 = static_cast<int32>(gen());
    int64 i2 = static_cast<int32>(gen());
    eint64 e1(i1), e2(i2), e3(i1 * i2);
    BOOST_TEST(e1 * e2 == e3);
  }
}

void extended_int_test9()
{
  typedef extended_int<1> eint32;
  for (int i = -10; i <= 10; ++i) {
    for (int j = -10; j <= 10; ++j) {
      eint32 e1(i), e2(j), e3(i+j), e4(i-j), e5(i*j);
      BOOST_TEST(e1 + e2 == e3);
      BOOST_TEST(e1 - e2 == e4);
      BOOST_TEST(e1 * e2 == e5);
    }
  }
}

void extended_int_test10()
{
  typedef extended_int<2> eint64;
  boost::mt19937_64 gen(static_cast<uint32>(time(NULL)));
  for (int i = 0; i < 100; ++i) {
    int64 i1 = static_cast<int64>(gen()) >> 20;
    int64 i2 = i1 >> 32;
    eint64 e1(i1), e2(i2);
    BOOST_TEST(to_fpt(e1) == static_cast<fpt64>(i1));
    BOOST_TEST(to_fpt(e2) == static_cast<fpt64>(i2));
  }
}

void extened_int_test11()
{
  typedef extended_int<64> eint2048;
  eint2048 two(2), value(1);
  for (int i = 0; i < 1024; ++i)
    value = value * two;
  BOOST_TEST_EQ(value.count(), 33);
  for (std::size_t i = 1; i < value.size(); ++i)
    BOOST_TEST_EQ(value.chunks()[i-1], 0U);
  BOOST_TEST_EQ(value.chunks()[32], 1U);
}

int main()
{
    ulp_comparison_test1();
    ulp_comparison_test2();
    extended_exponent_fpt_test1();
    extended_exponent_fpt_test2();
    extended_exponent_fpt_test3();
    extended_exponent_fpt_test4();
    extended_exponent_fpt_test5();
    extended_exponent_fpt_test6();
    extended_int_test1();
    extended_int_test2();
    extended_int_test3();
    extended_int_test4();
    extended_int_test5();
    extended_int_test6();
    extended_int_test7();
    extended_int_test8();
    extended_int_test9();
    extended_int_test10();
    extened_int_test11();
    return boost::report_errors();
}
Commit	Line	Data
7c673cae FG	1	// Boost.Polygon library voronoi_ctypes_test.cpp file
	2
	3	// Copyright Andrii Sydorchuk 2010-2012.
	4	// Distributed under the Boost Software License, Version 1.0.
	5	// (See accompanying file LICENSE_1_0.txt or copy at
	6	// http://www.boost.org/LICENSE_1_0.txt)
	7
	8	// See http://www.boost.org for updates, documentation, and revision history.
	9
92f5a8d4 TL	10	#include <boost/core/lightweight_test.hpp>
92f5a8d4 TL	11	#include <boost/polygon/detail/voronoi_ctypes.hpp>
7c673cae	12	#include <boost/random/mersenne_twister.hpp>
92f5a8d4 TL	13	#include <vector>
92f5a8d4 TL	14	#include <ctime>
7c673cae	15
7c673cae FG	16	using namespace boost::polygon::detail;
	17
	18	type_converter_fpt to_fpt;
	19
92f5a8d4 TL	20	void ulp_comparison_test1()
92f5a8d4 TL	21	{
7c673cae FG	22	ulp_comparison<double> ulp_cmp;
	23	uint64 a = 22;
	24	uint64 b = 27;
	25	fpt64 da, db;
	26	std::memcpy(&da, &a, sizeof(uint64));
	27	std::memcpy(&db, &b, sizeof(uint64));
92f5a8d4 TL	28	BOOST_TEST_EQ(ulp_cmp(da, db, 1), ulp_cmp.LESS);
	29	BOOST_TEST_EQ(ulp_cmp(db, da, 1), ulp_cmp.MORE);
	30	BOOST_TEST_EQ(ulp_cmp(da, db, 4), ulp_cmp.LESS);
	31	BOOST_TEST_EQ(ulp_cmp(da, db, 5), ulp_cmp.EQUAL);
	32	BOOST_TEST_EQ(ulp_cmp(da, db, 6), ulp_cmp.EQUAL);
7c673cae FG	33	}
7c673cae FG	34
92f5a8d4 TL	35	void ulp_comparison_test2()
92f5a8d4 TL	36	{
7c673cae FG	37	ulp_comparison<fpt64> ulp_cmp;
	38	uint64 a = 0ULL;
	39	uint64 b = 0x8000000000000002ULL;
	40	fpt64 da, db;
	41	std::memcpy(&da, &a, sizeof(uint64));
	42	std::memcpy(&db, &b, sizeof(uint64));
92f5a8d4 TL	43	BOOST_TEST_EQ(ulp_cmp(da, db, 1), ulp_cmp.MORE);
	44	BOOST_TEST_EQ(ulp_cmp(db, da, 1), ulp_cmp.LESS);
	45	BOOST_TEST_EQ(ulp_cmp(da, db, 2), ulp_cmp.EQUAL);
	46	BOOST_TEST_EQ(ulp_cmp(da, db, 3), ulp_cmp.EQUAL);
7c673cae FG	47	}
7c673cae FG	48
92f5a8d4 TL	49	void extended_exponent_fpt_test1()
92f5a8d4 TL	50	{
7c673cae FG	51	boost::mt19937_64 gen(static_cast<uint32>(time(NULL)));
	52	fpt64 b = 0.0;
	53	efpt64 eeb(b);
	54	for (int i = 0; i < 1000; ++i) {
	55	fpt64 a = to_fpt(static_cast<int64>(gen()));
	56	efpt64 eea(a);
	57	efpt64 neg = -eea;
	58	efpt64 sum = eea + eeb;
	59	efpt64 dif = eea - eeb;
	60	efpt64 mul = eea * eeb;
92f5a8d4 TL	61	BOOST_TEST_EQ(to_fpt(neg), -a);
	62	BOOST_TEST_EQ(to_fpt(sum), a + b);
	63	BOOST_TEST_EQ(to_fpt(dif), a - b);
	64	BOOST_TEST_EQ(to_fpt(mul), a * b);
7c673cae FG	65	}
	66	}
	67
92f5a8d4 TL	68	void extended_exponent_fpt_test2()
92f5a8d4 TL	69	{
7c673cae FG	70	boost::mt19937_64 gen(static_cast<uint32>(time(NULL)));
	71	fpt64 a = 0.0;
	72	efpt64 eea(a);
	73	for (int i = 0; i < 1000; ++i) {
	74	fpt64 b = to_fpt(static_cast<int64>(gen()));
	75	if (b == 0.0)
	76	continue;
	77	efpt64 eeb(b);
	78	efpt64 neg = -eea;
	79	efpt64 sum = eea + eeb;
	80	efpt64 dif = eea - eeb;
	81	efpt64 mul = eea * eeb;
	82	efpt64 div = eea / eeb;
92f5a8d4 TL	83	BOOST_TEST_EQ(to_fpt(neg), -a);
	84	BOOST_TEST_EQ(to_fpt(sum), a + b);
	85	BOOST_TEST_EQ(to_fpt(dif), a - b);
	86	BOOST_TEST_EQ(to_fpt(mul), a * b);
	87	BOOST_TEST_EQ(to_fpt(div), a / b);
7c673cae FG	88	}
	89	}
	90
92f5a8d4 TL	91	void extended_exponent_fpt_test3()
92f5a8d4 TL	92	{
7c673cae FG	93	boost::mt19937_64 gen(static_cast<uint32>(time(NULL)));
	94	for (int i = 0; i < 1000; ++i) {
	95	fpt64 a = to_fpt(static_cast<int64>(gen()));
	96	fpt64 b = to_fpt(static_cast<int64>(gen()));
	97	if (b == 0.0)
	98	continue;
	99	efpt64 eea(a);
	100	efpt64 eeb(b);
	101	efpt64 neg = -eea;
	102	efpt64 sum = eea + eeb;
	103	efpt64 dif = eea - eeb;
	104	efpt64 mul = eea * eeb;
	105	efpt64 div = eea / eeb;
92f5a8d4 TL	106	BOOST_TEST_EQ(to_fpt(neg), -a);
	107	BOOST_TEST_EQ(to_fpt(sum), a + b);
	108	BOOST_TEST_EQ(to_fpt(dif), a - b);
	109	BOOST_TEST_EQ(to_fpt(mul), a * b);
	110	BOOST_TEST_EQ(to_fpt(div), a / b);
7c673cae FG	111	}
	112	}
	113
92f5a8d4 TL	114	void extended_exponent_fpt_test4()
92f5a8d4 TL	115	{
7c673cae FG	116	for (int exp = 0; exp < 64; ++exp)
	117	for (int i = 1; i < 100; ++i) {
	118	fpt64 a = i;
	119	fpt64 b = to_fpt(1LL << exp);
	120	efpt64 eea(a);
	121	efpt64 eeb(b);
	122	efpt64 neg = -eea;
	123	efpt64 sum = eea + eeb;
	124	efpt64 dif = eea - eeb;
	125	efpt64 mul = eea * eeb;
	126	efpt64 div = eea / eeb;
92f5a8d4 TL	127	BOOST_TEST_EQ(to_fpt(neg), -a);
	128	BOOST_TEST_EQ(to_fpt(sum), a + b);
	129	BOOST_TEST_EQ(to_fpt(dif), a - b);
	130	BOOST_TEST_EQ(to_fpt(mul), a * b);
	131	BOOST_TEST_EQ(to_fpt(div), a / b);
7c673cae FG	132	}
	133	}
	134
92f5a8d4 TL	135	void extended_exponent_fpt_test5()
92f5a8d4 TL	136	{
7c673cae FG	137	for (int i = 0; i < 100; ++i) {
	138	efpt64 a(to_fpt(i * i));
	139	efpt64 b = a.sqrt();
92f5a8d4	140	BOOST_TEST_EQ(to_fpt(b), to_fpt(i));
7c673cae FG	141	}
	142	}
	143
92f5a8d4 TL	144	void extended_exponent_fpt_test6()
92f5a8d4 TL	145	{
7c673cae FG	146	for (int i = -10; i <= 10; ++i) {
7c673cae FG	147	efpt64 a(to_fpt(i));
92f5a8d4 TL	148	BOOST_TEST_EQ(is_pos(a), i > 0);
	149	BOOST_TEST_EQ(is_neg(a), i < 0);
	150	BOOST_TEST_EQ(is_zero(a), !i);
7c673cae FG	151	}
	152	}
	153
92f5a8d4 TL	154	void extended_int_test1()
92f5a8d4 TL	155	{
7c673cae FG	156	typedef extended_int<1> eint32;
7c673cae FG	157	eint32 e1(0), e2(32), e3(-32);
92f5a8d4 TL	158	BOOST_TEST_EQ(e1.count(), 0);
	159	BOOST_TEST_EQ(e1.size(), 0U);
	160	BOOST_TEST_EQ(e2.count(), 1);
	161	BOOST_TEST_EQ(e2.chunks()[0], 32U);
	162	BOOST_TEST_EQ(e2.size(), 1U);
	163	BOOST_TEST_EQ(e3.count(), -1);
	164	BOOST_TEST_EQ(e3.chunks()[0], 32U);
	165	BOOST_TEST_EQ(e3.size(), 1U);
7c673cae FG	166	}
7c673cae FG	167
92f5a8d4 TL	168	void extended_int_test2()
92f5a8d4 TL	169	{
7c673cae FG	170	typedef extended_int<2> eint64;
	171	int64 val64 = 0x7fffffffffffffffLL;
	172	eint64 e1(0), e2(32), e3(-32), e4(val64), e5(-val64);
92f5a8d4 TL	173	BOOST_TEST_EQ(e1.count(), 0);
	174	BOOST_TEST_EQ(e2.count(), 1);
	175	BOOST_TEST_EQ(e2.chunks()[0], 32U);
	176	BOOST_TEST_EQ(e3.count(), -1);
	177	BOOST_TEST_EQ(e3.chunks()[0], 32U);
	178	BOOST_TEST_EQ(e4.count(), 2);
	179	BOOST_TEST_EQ(e4.chunks()[0], 0xffffffff);
	180	BOOST_TEST_EQ(e4.chunks()[1], val64 >> 32);
	181	BOOST_TEST_EQ(e5.count(), -2);
	182	BOOST_TEST_EQ(e5.chunks()[0], 0xffffffff);
	183	BOOST_TEST_EQ(e5.chunks()[1], val64 >> 32);
7c673cae FG	184	}
7c673cae FG	185
92f5a8d4 TL	186	void extended_int_test3()
92f5a8d4 TL	187	{
7c673cae FG	188	typedef extended_int<2> eint64;
	189	std::vector<uint32> chunks;
	190	chunks.push_back(1);
	191	chunks.push_back(2);
	192	eint64 e1(chunks, true), e2(chunks, false);
92f5a8d4 TL	193	BOOST_TEST_EQ(e1.count(), 2);
	194	BOOST_TEST_EQ(e1.chunks()[0], 2U);
	195	BOOST_TEST_EQ(e1.chunks()[1], 1U);
	196	BOOST_TEST_EQ(e2.count(), -2);
	197	BOOST_TEST_EQ(e2.chunks()[0], 2U);
	198	BOOST_TEST_EQ(e2.chunks()[1], 1U);
7c673cae FG	199	}
7c673cae FG	200
92f5a8d4 TL	201	void extended_int_test4()
92f5a8d4 TL	202	{
7c673cae FG	203	typedef extended_int<2> eint64;
	204	std::vector<uint32> chunks;
	205	chunks.push_back(1);
	206	chunks.push_back(2);
	207	eint64 e1(chunks, true), e2(chunks, false);
92f5a8d4 TL	208	BOOST_TEST_EQ(e1 == e2, false);
	209	BOOST_TEST_EQ(e1 == -e2, true);
	210	BOOST_TEST_EQ(e1 != e2, true);
	211	BOOST_TEST_EQ(e1 != -e2, false);
	212	BOOST_TEST_EQ(e1 < e2, false);
	213	BOOST_TEST_EQ(e1 < -e2, false);
	214	BOOST_TEST_EQ(e1 <= e2, false);
	215	BOOST_TEST_EQ(e1 <= -e2, true);
	216	BOOST_TEST_EQ(e1 > e2, true);
	217	BOOST_TEST_EQ(e1 > -e2, false);
	218	BOOST_TEST_EQ(e1 >= e2, true);
	219	BOOST_TEST_EQ(e1 >= -e2, true);
7c673cae FG	220	}
7c673cae FG	221
92f5a8d4 TL	222	void extended_int_test5()
92f5a8d4 TL	223	{
7c673cae FG	224	typedef extended_int<2> eint64;
	225	boost::mt19937_64 gen(static_cast<uint32>(time(NULL)));
	226	for (int i = 0; i < 1000; ++i) {
	227	int64 i1 = static_cast<int64>(gen());
	228	int64 i2 = static_cast<int64>(gen());
	229	eint64 e1(i1), e2(i2);
92f5a8d4 TL	230	BOOST_TEST_EQ(e1 == e2, i1 == i2);
	231	BOOST_TEST_EQ(e1 != e2, i1 != i2);
	232	BOOST_TEST_EQ(e1 > e2, i1 > i2);
	233	BOOST_TEST_EQ(e1 >= e2, i1 >= i2);
	234	BOOST_TEST_EQ(e1 < e2, i1 < i2);
	235	BOOST_TEST_EQ(e1 <= e2, i1 <= i2);
7c673cae FG	236	}
	237	}
	238
92f5a8d4 TL	239	void extended_int_test6()
92f5a8d4 TL	240	{
7c673cae FG	241	typedef extended_int<1> eint32;
	242	eint32 e1(32);
	243	eint32 e2 = -e1;
92f5a8d4 TL	244	BOOST_TEST_EQ(e2.count(), -1);
	245	BOOST_TEST_EQ(e2.size(), 1U);
	246	BOOST_TEST_EQ(e2.chunks()[0], 32U);
7c673cae FG	247	}
7c673cae FG	248
92f5a8d4 TL	249	void extended_int_test7()
92f5a8d4 TL	250	{
7c673cae FG	251	typedef extended_int<2> eint64;
	252	boost::mt19937_64 gen(static_cast<uint32>(time(NULL)));
	253	for (int i = 0; i < 1000; ++i) {
	254	int64 i1 = static_cast<int64>(gen()) >> 2;
	255	int64 i2 = static_cast<int64>(gen()) >> 2;
	256	eint64 e1(i1), e2(i2), e3(i1 + i2), e4(i1 - i2);
92f5a8d4 TL	257	BOOST_TEST(e1 + e2 == e3);
92f5a8d4 TL	258	BOOST_TEST(e1 - e2 == e4);
7c673cae FG	259	}
	260	}
	261
92f5a8d4 TL	262	void extended_int_test8()
92f5a8d4 TL	263	{
7c673cae FG	264	typedef extended_int<2> eint64;
	265	boost::mt19937 gen(static_cast<uint32>(time(NULL)));
	266	for (int i = 0; i < 1000; ++i) {
	267	int64 i1 = static_cast<int32>(gen());
	268	int64 i2 = static_cast<int32>(gen());
	269	eint64 e1(i1), e2(i2), e3(i1 * i2);
92f5a8d4	270	BOOST_TEST(e1 * e2 == e3);
7c673cae FG	271	}
	272	}
	273
92f5a8d4 TL	274	void extended_int_test9()
92f5a8d4 TL	275	{
7c673cae FG	276	typedef extended_int<1> eint32;
	277	for (int i = -10; i <= 10; ++i) {
	278	for (int j = -10; j <= 10; ++j) {
	279	eint32 e1(i), e2(j), e3(i+j), e4(i-j), e5(i*j);
92f5a8d4 TL	280	BOOST_TEST(e1 + e2 == e3);
	281	BOOST_TEST(e1 - e2 == e4);
	282	BOOST_TEST(e1 * e2 == e5);
7c673cae FG	283	}
	284	}
	285	}
	286
92f5a8d4 TL	287	void extended_int_test10()
92f5a8d4 TL	288	{
7c673cae FG	289	typedef extended_int<2> eint64;
	290	boost::mt19937_64 gen(static_cast<uint32>(time(NULL)));
	291	for (int i = 0; i < 100; ++i) {
	292	int64 i1 = static_cast<int64>(gen()) >> 20;
	293	int64 i2 = i1 >> 32;
	294	eint64 e1(i1), e2(i2);
92f5a8d4 TL	295	BOOST_TEST(to_fpt(e1) == static_cast<fpt64>(i1));
92f5a8d4 TL	296	BOOST_TEST(to_fpt(e2) == static_cast<fpt64>(i2));
7c673cae FG	297	}
	298	}
	299
92f5a8d4 TL	300	void extened_int_test11()
92f5a8d4 TL	301	{
7c673cae FG	302	typedef extended_int<64> eint2048;
	303	eint2048 two(2), value(1);
	304	for (int i = 0; i < 1024; ++i)
	305	value = value * two;
92f5a8d4	306	BOOST_TEST_EQ(value.count(), 33);
7c673cae	307	for (std::size_t i = 1; i < value.size(); ++i)
92f5a8d4 TL	308	BOOST_TEST_EQ(value.chunks()[i-1], 0U);
	309	BOOST_TEST_EQ(value.chunks()[32], 1U);
	310	}
	311
	312	int main()
	313	{
	314	ulp_comparison_test1();
	315	ulp_comparison_test2();
	316	extended_exponent_fpt_test1();
	317	extended_exponent_fpt_test2();
	318	extended_exponent_fpt_test3();
	319	extended_exponent_fpt_test4();
	320	extended_exponent_fpt_test5();
	321	extended_exponent_fpt_test6();
	322	extended_int_test1();
	323	extended_int_test2();
	324	extended_int_test3();
	325	extended_int_test4();
	326	extended_int_test5();
	327	extended_int_test6();
	328	extended_int_test7();
	329	extended_int_test8();
	330	extended_int_test9();
	331	extended_int_test10();
	332	extened_int_test11();
	333	return boost::report_errors();
7c673cae	334	}