1 // (C) Copyright John Maddock 2007.
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/special_functions/expint.hpp>
7 #include <boost/math/constants/constants.hpp>
8 #include <fstream>
9 #include <boost/math/tools/test_data.hpp>
10 #include "mp_t.hpp"
11
12 using namespace boost::math::tools;
13
14 struct expint_data_generator
15 {
operator ()expint_data_generator16 mp_t operator()(mp_t a, mp_t b)
17 {
18 unsigned n = boost::math::tools::real_cast<unsigned>(a);
19 std::cout << n << " " << b << " ";
20 mp_t result = boost::math::expint(n, b);
21 std::cout << result << std::endl;
22 return result;
23 }
24 };
25
26
main()27 int main()
28 {
29 boost::math::expint(1, 0.06227754056453704833984375);
30 std::cout << boost::math::expint(1, mp_t(0.5)) << std::endl;
31
32 parameter_info<mp_t> arg1, arg2;
33 test_data<mp_t> data;
34
35 std::cout << "Welcome.\n"
36 "This program will generate spot tests for the expint function:\n"
37 " expint(a, b)\n\n";
38
39 bool cont;
40 std::string line;
41
42 do{
43 get_user_parameter_info(arg1, "a");
44 get_user_parameter_info(arg2, "b");
45 data.insert(expint_data_generator(), arg1, arg2);
46
47 std::cout << "Any more data [y/n]?";
48 std::getline(std::cin, line);
49 boost::algorithm::trim(line);
50 cont = (line == "y");
51 }while(cont);
52
53 std::cout << "Enter name of test data file [default=expint_data.ipp]";
54 std::getline(std::cin, line);
55 boost::algorithm::trim(line);
56 if(line == "")
57 line = "expint_data.ipp";
58 std::ofstream ofs(line.c_str());
59 ofs << std::scientific << std::setprecision(40);
60 write_code(ofs, data, "expint_data");
61
62 return 0;
63 }
64
65