Blame view

3rdparty/boost_1_81_0/libs/multiprecision/tools/sincos.cpp 1.75 KB
73ef4ff3   Hu Chunming   提交三方库
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
  // Copyright John Maddock 2013.
  // Use, modification and distribution are subject to 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)
  
  //[special_data_example
  
  #include <boost/multiprecision/mpfr.hpp>
  #include <boost/math/tools/test_data.hpp>
  #include <boost/test/included/prg_exec_monitor.hpp>
  #include <boost/math/tools/tuple.hpp>
  #include <fstream>
  
  using namespace boost::math::tools;
  using namespace boost::math;
  using namespace std;
  using namespace boost::multiprecision;
  
  typedef number<mpfr_float_backend<1000> > mp_type;
  
  
  boost::math::tuple<mp_type, mp_type, mp_type> generate(mp_type a)
  {
     mp_type ts, tc;
     ts = sin(a);
     tc = cos(a);
     return boost::math::make_tuple(a, ts, tc);
  }
  
  int cpp_main(int argc, char*argv [])
  {
     parameter_info<mp_type> arg1, arg2;
     test_data<mp_type> data;
  
     bool cont;
     std::string line;
  
     if(argc < 1)
        return 1;
  
     do{
        //
        // User interface which prompts for 
        // range of input parameters:
        //
        if(0 == get_user_parameter_info(arg1, "a"))
           return 1;
        arg1.type |= dummy_param;
  
        data.insert(&generate, arg1);
  
        std::cout << "Any more data [y/n]?";
        std::getline(std::cin, line);
        boost::algorithm::trim(line);
        cont = (line == "y");
     }while(cont);
     //
     // Just need to write the results to a file:
     //
     std::cout << "Enter name of test data file [default=sincos.ipp]";
     std::getline(std::cin, line);
     boost::algorithm::trim(line);
     if(line == "")
        line = "sincos.ipp";
     std::ofstream ofs(line.c_str());
     line.erase(line.find('.'));
     ofs << std::scientific << std::setprecision(500);
     write_code(ofs, data, line.c_str());
  
     return 0;
  }
  
  //]