text: ///////////////////////////////////////////////////////////////
  //  Copyright Christopher Kormanyos 2002 - 2011.
  //  Copyright 2011 John Maddock. Distributed under the Boost
  //  Software License, Version 1.0. (See accompanying file
  //  LICENSE_1_0.txt or copy at https://www.boost.org/LICENSE_1_0.txt
  //
  // This work is based on an earlier work:
  // "Algorithm 910: A Portable C++ Multiple-Precision System for Special-Function Calculations",
  // in ACM TOMS, {VOL 37, ISSUE 4, (February 2011)} (C) ACM, 2011. http://doi.acm.org/10.1145/1916461.1916469
  
  #include "setup.hpp"
  #include "table_type.hpp"
  
  #include <boost/math/special_functions/math_fwd.hpp>
  #include "libs/math/test/test_bessel_k.hpp"
  
  void expected_results() {
      //
      // Define the max and mean errors expected for
      // various compilers and platforms.
      //
      add_expected_result(".*",                 // compiler
                          ".*",                 // stdlib
                          ".*",                 // platform
                          ".*gmp.*",            // test type(s)
                          ".*",                 // test data group
                          ".*", 2000, 1500);    // test function
  #ifdef BOOST_INTEL
      add_expected_result(".*",               // compiler
                          ".*",               // stdlib
                          ".*",               // platform
                          ".*float128.*",     // test type(s)
                          ".*",               // test data group
                          ".*", 300, 100);    // test function
  #endif
      add_expected_result(".*",                            // compiler
                          ".*",                            // stdlib
                          ".*",                            // platform
                          ".*mpfr_float_backend<18>.*",    // test type(s)
                          ".*",                            // test data group
                          ".*", 3000, 1000);               // test function
      add_expected_result(".*",                            // compiler
                          ".*",                            // stdlib
                          ".*",                            // platform
                          ".*mpfr_float_backend<0>.*",     // test type(s)
                          ".*",                            // test data group
                          ".*", 300, 150);                 // test function
      add_expected_result(".*",                            // compiler
                          ".*",                            // stdlib
                          ".*",                            // platform
                          ".*",                            // test type(s)
                          ".*large.*",                     // test data group
                          ".*", 80, 50);                   // test function
      add_expected_result(".*",                            // compiler
                          ".*",                            // stdlib
                          ".*",                            // platform
                          ".*cpp_bin_float.*",             // test type(s)
                          ".*",                            // test data group
                          ".*", 300, 150);                 // test function
      add_expected_result(".*",                            // compiler
                          ".*",                            // stdlib
                          ".*",                            // platform
                          ".*",                            // test type(s)
                          ".*",                            // test data group
                          ".*", 50, 15);                   // test function
      //
      // Finish off by printing out the compiler/stdlib/platform names,
      // we do this to make it easier to mark up expected error rates.
      //
      std::cout << "Tests run with " << BOOST_COMPILER << ", " << BOOST_STDLIB << ", " << BOOST_PLATFORM << std::endl;
  }
  
  template<class T>
  void test(T t, const char* p) {
      test_bessel(t, p);
  }
  
  BOOST_AUTO_TEST_CASE(test_main) {
      using namespace nil::crypto3::multiprecision;
      expected_results();
      //
      // Test at:
      // 18 decimal digits: tests 80-bit long double approximations
      // 30 decimal digits: tests 128-bit long double approximations
      // 35 decimal digits: tests arbitrary precision code
      //
      ALL_TESTS
  }
  