12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273 |
- // Copyright John Maddock 2007.
- // Copyright Paul A. Bristow 2010
- // 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)
- // Note that this file contains quickbook mark-up as well as code
- // and comments, don't change any of the special comment mark-ups!
- #include <iostream>
- using std::cout; using std::endl;
- #include <cerrno> // for ::errno
- //[policy_eg_7
- #include <boost/math/distributions.hpp> // All distributions.
- // using boost::math::normal; // Would create an ambguity between
- // boost::math::normal_distribution<RealType> boost::math::normal and
- // 'anonymous-namespace'::normal'.
- namespace
- { // anonymous or unnnamed (rather than named as in policy_eg_6.cpp).
- using namespace boost::math::policies;
- // using boost::math::policies::errno_on_error; // etc.
- typedef policy<
- // return infinity and set errno rather than throw:
- overflow_error<errno_on_error>,
- // Don't promote double -> long double internally:
- promote_double<false>,
- // Return the closest integer result for discrete quantiles:
- discrete_quantile<integer_round_nearest>
- > my_policy;
- BOOST_MATH_DECLARE_DISTRIBUTIONS(double, my_policy)
- } // close namespace my_namespace
- int main()
- {
- // Construct distribution with something we know will overflow.
- normal norm(10, 2); // using 'anonymous-namespace'::normal
- errno = 0;
- cout << "Result of quantile(norm, 0) is: "
- << quantile(norm, 0) << endl;
- cout << "errno = " << errno << endl;
- errno = 0;
- cout << "Result of quantile(norm, 1) is: "
- << quantile(norm, 1) << endl;
- cout << "errno = " << errno << endl;
- //
- // Now try a discrete distribution:
- binomial binom(20, 0.25);
- cout << "Result of quantile(binom, 0.05) is: "
- << quantile(binom, 0.05) << endl;
- cout << "Result of quantile(complement(binom, 0.05)) is: "
- << quantile(complement(binom, 0.05)) << endl;
- }
- //] //[/policy_eg_7]
- /*
- Output:
- Result of quantile(norm, 0) is: -1.#INF
- errno = 34
- Result of quantile(norm, 1) is: 1.#INF
- errno = 34
- Result of quantile(binom, 0.05) is: 1
- Result of quantile(complement(binom, 0.05)) is: 8
- */
|