ceph/src/boost/libs/math/tools/gauss_kronrod_constants.cpp

   1 //  Copyright (c) 2017 John Maddock
   2 //  Use, modification and distribution are subject to the
   3 //  Boost Software License, Version 1.0. (See accompanying file
   4 //  LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
   5
   6 #include <boost/math/quadrature/gauss_kronrod.hpp>
   7 #include <boost/multiprecision/cpp_bin_float.hpp>
   8
   9 template <class T>
  10 void print_gauss_constants(const char* suffix, int prec, int tag)
  11 {
  12    auto ab = T::abscissa();
  13    auto w = T::weights();
  14    std::cout << std::setprecision(prec) << std::scientific;
  15    std::size_t order = (ab[0] == 0) ? (ab.size() * 2) - 1 : ab.size() * 2;
  16    std::cout <<
  17       "template <class T>\n"
  18       "class gauss_detail<T, " << order << ", " << tag << ">\n"
  19       "   {\n"
  20       "   public:\n"
  21       "      static " << (prec > 40 ? " " : "constexpr ") << "std::array<T, " << ab.size() << "> const & abscissa()\n"
  22       "      {\n"
  23       "         static " << (prec > 40 ? " " : "constexpr ") << "std::array<T, " << ab.size() << "> data = {\n";
  24    for (unsigned i = 0; i < ab.size(); ++i)
  25       std::cout << "            " << (prec > 40 ? "BOOST_MATH_HUGE_CONSTANT(T, 0, " : "") << ab[i] << (prec > 40 ? ")" : suffix) << ",\n";
  26    std::cout <<
  27       "};\n"
  28       "         return data;\n"
  29       "      }\n"
  30       "      static " << (prec > 40 ? " " : "constexpr ") << "std::array<T, " << w.size() << "> const & weights()\n"
  31       "      {\n"
  32       "         static " << (prec > 40 ? " " : "constexpr ") << "std::array<T, " << w.size() << "> data = {\n";
  33    for (unsigned i = 0; i < w.size(); ++i)
  34       std::cout << "            " << (prec > 40 ? "BOOST_MATH_HUGE_CONSTANT(T, 0, " : "") << w[i] << (prec > 40 ? ")" : suffix) << ",\n";
  35
  36    std::cout << "         };\n"
  37       "         return data;\n"
  38       "      }\n"
  39       "   };\n\n";
  40 }
  41
  42 template <class T>
  43 void print_gauss_kronrod_constants(const char* suffix, int prec, int tag)
  44 {
  45    auto ab = T::abscissa();
  46    auto w = T::weights();
  47    std::cout << std::setprecision(prec) << std::scientific;
  48    std::size_t order = (ab.size() * 2) - 1;
  49    std::cout <<
  50       "   template <class T>\n"
  51       "   class gauss_kronrod_detail<T, " << order << ", " << tag << ">\n"
  52       "   {\n"
  53       "   public:\n"
  54       "      static " << (prec > 40 ? " " : "constexpr ") << "std::array<T, " << ab.size() << "> const & abscissa()\n"
  55       "      {\n"
  56       "         static " << (prec > 40 ? " " : "constexpr ") << "std::array<T, " << ab.size() << "> data = {\n";
  57
  58    for (unsigned i = 0; i < ab.size(); ++i)
  59       std::cout << "            " << (prec > 40 ? "BOOST_MATH_HUGE_CONSTANT(T, 0, " : "") << ab[i] << (prec > 40 ? ")" : suffix) << ",\n";
  60
  61    std::cout << "         };\n"
  62       "         return data;\n"
  63       "      }\n"
  64       "      static " << (prec > 40 ? " " : "constexpr ") << "std::array<T, " << w.size() << "> const & weights()\n"
  65       "      {\n"
  66       "         static " << (prec > 40 ? " " : "constexpr ") << "std::array<T, " << w.size() << "> data = {\n";
  67
  68    for (unsigned i = 0; i < w.size(); ++i)
  69       std::cout << "            " << (prec > 40 ? "BOOST_MATH_HUGE_CONSTANT(T, 0, " : "") << w[i] << (prec > 40 ? ")" : suffix) << ",\n";
  70
  71    std::cout << "         };\n"
  72       "         return data;\n"
  73       "      }\n"
  74       "   };\n\n";
  75 }
  76
  77
  78
  79 int main()
  80 {
  81    typedef boost::multiprecision::number<boost::multiprecision::cpp_bin_float<250> > mp_type;
  82
  83    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 7> >("f", 9, 0);
  84    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 7> >("", 17, 1);
  85    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 7> >("L", 35, 2);
  86    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 7> >("Q", 35, 3);
  87    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 7> >("", 115, 4);
  88
  89    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 10> >("f", 9, 0);
  90    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 10> >("", 17, 1);
  91    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 10> >("L", 35, 2);
  92    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 10> >("Q", 35, 3);
  93    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 10> >("", 115, 4);
  94
  95    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 15> >("f", 9, 0);
  96    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 15> >("", 17, 1);
  97    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 15> >("L", 35, 2);
  98    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 15> >("Q", 35, 3);
  99    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 15> >("", 115, 4);
 100
 101    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 20> >("f", 9, 0);
 102    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 20> >("", 17, 1);
 103    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 20> >("L", 35, 2);
 104    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 20> >("Q", 35, 3);
 105    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 20> >("", 115, 4);
 106
 107    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 25> >("f", 9, 0);
 108    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 25> >("", 17, 1);
 109    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 25> >("L", 35, 2);
 110    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 25> >("Q", 35, 3);
 111    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 25> >("", 115, 4);
 112
 113    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 30> >("f", 9, 0);
 114    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 30> >("", 17, 1);
 115    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 30> >("L", 35, 2);
 116    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 30> >("Q", 35, 3);
 117    print_gauss_constants<boost::math::quadrature::gauss<mp_type, 30> >("", 115, 4);
 118
 119
 120    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 15> >("f", 9, 0);
 121    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 15> >("", 17, 1);
 122    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 15> >("L", 35, 2);
 123    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 15> >("Q", 35, 3);
 124    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 15> >("", 115, 4);
 125
 126    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 21> >("f", 9, 0);
 127    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 21> >("", 17, 1);
 128    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 21> >("L", 35, 2);
 129    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 21> >("Q", 35, 3);
 130    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 21> >("", 115, 4);
 131
 132    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 31> >("f", 9, 0);
 133    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 31> >("", 17, 1);
 134    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 31> >("L", 35, 2);
 135    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 31> >("Q", 35, 3);
 136    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 31> >("", 115, 4);
 137
 138    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 41> >("f", 9, 0);
 139    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 41> >("", 17, 1);
 140    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 41> >("L", 35, 2);
 141    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 41> >("Q", 35, 3);
 142    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 41> >("", 115, 4);
 143
 144    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 51> >("f", 9, 0);
 145    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 51> >("", 17, 1);
 146    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 51> >("L", 35, 2);
 147    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 51> >("Q", 35, 3);
 148    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 51> >("", 115, 4);
 149
 150    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 61> >("f", 9, 0);
 151    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 61> >("", 17, 1);
 152    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 61> >("L", 35, 2);
 153    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 61> >("Q", 35, 3);
 154    print_gauss_kronrod_constants<boost::math::quadrature::gauss_kronrod<mp_type, 61> >("", 115, 4);
 155
 156    return 0;
 157 }
 158