]>
Commit | Line | Data |
---|---|---|
7c673cae FG |
1 | // Copyright Benjamin Sobotta 2012 |
2 | ||
3 | // Use, modification and distribution are subject to the | |
4 | // Boost Software License, Version 1.0. (See accompanying file | |
5 | // LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt) | |
6 | ||
7 | #ifdef _MSC_VER | |
8 | # pragma warning(disable : 4127) // conditional expression is constant. | |
9 | #endif | |
10 | ||
11 | #include <boost/math/special_functions/owens_t.hpp> | |
12 | #include <iostream> | |
13 | ||
14 | int main() | |
15 | { | |
16 | double h = 0.0,a; | |
17 | std::cout << std::setprecision(20); | |
18 | ||
19 | static const double a_vec[] = { | |
20 | 0.5000000000000000E+00, | |
21 | 0.1000000000000000E+01, | |
22 | 0.2000000000000000E+01, | |
23 | 0.3000000000000000E+01, | |
24 | 0.5000000000000000E+00, | |
25 | 0.1000000000000000E+01, | |
26 | 0.2000000000000000E+01, | |
27 | 0.3000000000000000E+01, | |
28 | 0.5000000000000000E+00, | |
29 | 0.1000000000000000E+01, | |
30 | 0.2000000000000000E+01, | |
31 | 0.3000000000000000E+01, | |
32 | 0.5000000000000000E+00, | |
33 | 0.1000000000000000E+01, | |
34 | 0.2000000000000000E+01, | |
35 | 0.3000000000000000E+01, | |
36 | 0.5000000000000000E+00, | |
37 | 0.1000000000000000E+01, | |
38 | 0.2000000000000000E+01, | |
39 | 0.3000000000000000E+01, | |
40 | 0.1000000000000000E+02, | |
41 | 0.1000000000000000E+03 }; | |
42 | ||
43 | static const double h_vec[] = { | |
44 | 0.1000000000000000E+01, | |
45 | 0.1000000000000000E+01, | |
46 | 0.1000000000000000E+01, | |
47 | 0.1000000000000000E+01, | |
48 | 0.5000000000000000E+00, | |
49 | 0.5000000000000000E+00, | |
50 | 0.5000000000000000E+00, | |
51 | 0.5000000000000000E+00, | |
52 | 0.2500000000000000E+00, | |
53 | 0.2500000000000000E+00, | |
54 | 0.2500000000000000E+00, | |
55 | 0.2500000000000000E+00, | |
56 | 0.1250000000000000E+00, | |
57 | 0.1250000000000000E+00, | |
58 | 0.1250000000000000E+00, | |
59 | 0.1250000000000000E+00, | |
60 | 0.7812500000000000E-02, | |
61 | 0.7812500000000000E-02, | |
62 | 0.7812500000000000E-02, | |
63 | 0.7812500000000000E-02, | |
64 | 0.7812500000000000E-02, | |
65 | 0.7812500000000000E-02 }; | |
66 | ||
67 | static const double t_vec[] = { | |
68 | 0.4306469112078537E-01, | |
69 | 0.6674188216570097E-01, | |
70 | 0.7846818699308410E-01, | |
71 | 0.7929950474887259E-01, | |
72 | 0.6448860284750376E-01, | |
73 | 0.1066710629614485E+00, | |
74 | 0.1415806036539784E+00, | |
75 | 0.1510840430760184E+00, | |
76 | 0.7134663382271778E-01, | |
77 | 0.1201285306350883E+00, | |
78 | 0.1666128410939293E+00, | |
79 | 0.1847501847929859E+00, | |
80 | 0.7317273327500385E-01, | |
81 | 0.1237630544953746E+00, | |
82 | 0.1737438887583106E+00, | |
83 | 0.1951190307092811E+00, | |
84 | 0.7378938035365546E-01, | |
85 | 0.1249951430754052E+00, | |
86 | 0.1761984774738108E+00, | |
87 | 0.1987772386442824E+00, | |
88 | 0.2340886964802671E+00, | |
89 | 0.2479460829231492E+00 }; | |
90 | ||
91 | for(unsigned i = 0; i != 22; ++i) | |
92 | { | |
93 | h = h_vec[i]; | |
94 | a = a_vec[i]; | |
95 | const double t = boost::math::owens_t(h, a); | |
96 | std::cout << "h=" << h << "\ta=" << a << "\tcomp=" | |
97 | << t << "\ttab=" << t_vec[i] | |
98 | << "\tdiff=" << std::fabs(t_vec[i]-t) << std::endl;; | |
99 | } | |
100 | ||
101 | return 0; | |
102 | } | |
103 | ||
104 | ||
105 | // EOF | |
106 | ||
107 | ||
108 | ||
109 | ||
110 | ||
111 | ||
112 | ||
113 |