Библиотека для работы с длинными числами (BigInteger) -> Форум на Исходниках.RU [Powered by Invision Power Board]

<{CODE_COLLAPSE_OFF}><{CODE_WRAP_OFF}>

#include <iostream>

#include <iterator>

#include <vector>

#include <algorithm>

#include <iomanip>

#include <boost/multiprecision/cpp_bin_float.hpp>

#include <chrono>

#include <cmath>

#include <typeinfo>

// вывод времени в микросекундах

template <typename Ch, typename Tr>

std::basic_ostream<Ch, Tr>& operator<<(std::basic_ostream<Ch, Tr>& os,

std::chrono::high_resolution_clock::duration time)

{

return os << std::chrono::duration_cast<std::chrono::microseconds>(time).count();

}

// точная тестовая функция

template <typename T>

T func(T x)

{

return pow(x-5, 6);

}

// тестовая функция с потерей точности

template <typename T>

T poly(T x)

{

return pow(x, 6) - 30 * pow(x, 5) + 375 * pow(x, 4) - 2500 * pow(x, 3) + 9375 * pow(x, 2) - 18750 * x + 15625;

}

// запуск тестовой функции на критическом интервале с малым шагом

template <typename T, T Fn(T), typename It>

void run(It it)

{

const T start = static_cast<T>(4.98),

end = static_cast<T>(5.02),

step = static_cast<T>(0.0001);

for (T f = start; f < end; f += step)

*it++ = Fn(f);

}

// запуск теста COUNT раз с замером времени

template <typename T, T Fn(T), typename Cont>

void test(Cont &cont)

{

std::cerr << typeid(T).name() << ":\n";

const int COUNT = 100;

auto start = std::chrono::high_resolution_clock::now();

for (int i = 0; i < COUNT; ++i)

{

cont.clear();

run<T, Fn>(std::back_inserter<Cont>(cont));

}

auto diff = std::chrono::high_resolution_clock::now() - start;

std::cerr << diff << " mcs" << std::endl;

}

// тестирование обоих тестовых функций с выводом их результатов

template <typename T> void boost_test()

{

std::vector<T> vT_1;

std::vector<T> vT_2;

test<T, func>(vT_1);

test<T, poly>(vT_2);

for (int i = 0; i < std::min(vT_1.size(), vT_2.size()); ++i)

std::cout << i << ":\t" << std::setw(20) << vT_1[i] << '\t'

<< std::setw(20) << vT_2[i] << '\n';

std::cout << std::endl;

}

int main()

{

using namespace boost::multiprecision;

typedef number<cpp_bin_float<50, digit_base_10>> cpp_dec_50;

typedef number<cpp_bin_float<1000, digit_base_2>> cpp_bin_float_1000;

boost_test<cpp_bin_float_double>();// IEEE's double precision

boost_test<cpp_bin_float_quad>(); // IEEE's quad precision

boost_test<cpp_dec_50>(); // 50-dec

boost_test<cpp_bin_float_1000>(); // 1000-bit

}