Boost-Commit :

Date view	Thread view	Subject view	Author view

Subject: [Boost-commit] svn:boost r77513 - in trunk: boost/math/special_functions/detail libs/math/doc/sf_and_dist libs/math/doc/sf_and_dist/html libs/math/doc/sf_and_dist/html/index libs/math/doc/sf_and_dist/html/math_toolkit/constants libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/binom_eg libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/find_eg libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/nccs_eg libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/neg_binom_eg libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/normal_example libs/math/doc/sf_and_dist/html/math_toolkit/main_overview libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial libs/math/doc/sf_and_dist/html/math_toolkit/status libs/math/doc/sf_and_dist/html/math_toolkit/utils/fp_facets libs/math/test
From: john_at_[hidden]
Date: 2012-03-24 07:59:26

Next message: jewillco_at_[hidden]: "[Boost-commit] svn:boost r77514 - in trunk: boost/property_map libs/property_map/doc"
Previous message: joel_at_[hidden]: "[Boost-commit] svn:boost r77512 - trunk/boost/spirit/include"

Author: johnmaddock
Date: 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
New Revision: 77513
URL: http://svn.boost.org/trac/boost/changeset/77513

Log:
Fix numeric error in bessel_jy.hpp.
Remove some dead code from bessel_jy.hpp.
Add additional test cases for above bug.
Update and regenerate docs to match.
Text files modified:
   trunk/boost/math/special_functions/detail/bessel_jy.hpp | 16 +++++---
   trunk/libs/math/doc/sf_and_dist/html/index.html | 2
   trunk/libs/math/doc/sf_and_dist/html/index/s13.html | 4 +-
   trunk/libs/math/doc/sf_and_dist/html/index/s14.html | 4 +-
   trunk/libs/math/doc/sf_and_dist/html/index/s15.html | 4 +-
   trunk/libs/math/doc/sf_and_dist/html/index/s16.html | 4 +-
   trunk/libs/math/doc/sf_and_dist/html/index/s17.html | 4 +-
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/constants/new_const.html | 1
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/binom_eg/binomial_coinflip_example.html | 39 ++++++++++++++-------
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/binom_eg/binomial_quiz_example.html | 71 +++++++++++++++++++++++++-------------
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/dist_construct_eg.html | 45 ++++++++++++++++--------
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/error_eg.html | 6 ++-
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/find_eg/find_location_eg.html | 22 +++++++----
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/find_eg/find_mean_and_sd_eg.html | 66 ++++++++++++++++++++++++------------
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/find_eg/find_scale_eg.html | 24 ++++++++----
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/geometric_eg.html | 66 ++++++++++++++++++++++++------------
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/inverse_chi_squared_eg.html | 6 +-
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/nccs_eg/nccs_power_eg.html | 18 ++++++---
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/neg_binom_eg/neg_binom_conf.html | 13 ++++--
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/neg_binom_eg/neg_binom_size_eg.html | 13 +++---
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/neg_binom_eg/negative_binomial_example1.html | 72 ++++++++++++++++++++++++++-------------
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/normal_example/normal_misc.html | 64 +++++++++++++++++++++++------------
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/main_overview/history1.html | 40 +++++++++++++--------
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/discrete_quant_ref.html | 4 --
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/error_handling_policies.html | 1
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/internal_promotion.html | 2 -
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/namespace_pol.html | 2 -
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/precision_pol.html | 2 -
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/ad_hoc_dist_policies.html | 1
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/ad_hoc_sf_policies.html | 2 -
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/changing_policy_defaults.html | 6 ++-
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/namespace_policies.html | 29 ++++++++++------
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/understand_dis_quant.html | 9 +++-
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/user_def_err_pol.html | 51 ++++++++++++++++++---------
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/status/history1.html | 40 +++++++++++++--------
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/utils/fp_facets/examples.html | 1
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/utils/fp_facets/intro.html | 1
   trunk/libs/math/doc/sf_and_dist/roadmap.qbk | 4 ++
   trunk/libs/math/test/pch.hpp | 7 +--
   trunk/libs/math/test/test_bessel_j.cpp | 10 +++++
   trunk/libs/math/test/test_bessel_j.hpp | 7 ++-
   trunk/libs/math/test/test_bessel_y.cpp | 10 +++++
   trunk/libs/math/test/test_bessel_y.hpp | 5 ++
   43 files changed, 510 insertions(+), 288 deletions(-)

Modified: trunk/boost/math/special_functions/detail/bessel_jy.hpp
==============================================================================
--- trunk/boost/math/special_functions/detail/bessel_jy.hpp (original)
+++ trunk/boost/math/special_functions/detail/bessel_jy.hpp 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -215,8 +215,6 @@
 Cr = br + a / temp;
 Dr = br;
 Di = bi;
- //std::cout << "C = " << Cr << " " << Ci << std::endl;
- //std::cout << "D = " << Dr << " " << Di << std::endl;
 if (fabs(Cr) + fabs(Ci) < tiny) { Cr = tiny; }
 if (fabs(Dr) + fabs(Di) < tiny) { Dr = tiny; }
 temp = Dr * Dr + Di * Di;
@@ -227,7 +225,6 @@
 temp = fr;
 fr = temp * delta_r - fi * delta_i;
 fi = temp * delta_i + fi * delta_r;
- //std::cout << fr << " " << fi << std::endl;
 for (k = 2; k < policies::get_max_series_iterations<Policy>(); k++)
 {
 a = k - 0.5f;
@@ -239,8 +236,6 @@
 Ci = bi - a * Ci / temp;
 Dr = br + a * Dr;
 Di = bi + a * Di;
- //std::cout << "C = " << Cr << " " << Ci << std::endl;
- //std::cout << "D = " << Dr << " " << Di << std::endl;
 if (fabs(Cr) + fabs(Ci) < tiny) { Cr = tiny; }
 if (fabs(Dr) + fabs(Di) < tiny) { Dr = tiny; }
 temp = Dr * Dr + Di * Di;
@@ -253,7 +248,6 @@
 fi = temp * delta_i + fi * delta_r;
 if (fabs(delta_r - 1) + fabs(delta_i) < tolerance)
 break;
- //std::cout << fr << " " << fi << std::endl;
 }
 policies::check_series_iterations<T>("boost::math::bessel_jy<%1%>(%1%,%1%) in CF2_jy", k, pol);
 *p = fr;
@@ -491,6 +485,16 @@
 CF2_jy(u, x, &p, &q, pol); // continued fraction CF2_jy
 T t = u / x - fu; // t = J'/J
 gamma = (p - t) / q;
+ //
+ // We can't allow gamma to cancel out to zero competely as it messes up
+ // the subsequent logic. So pretend that one bit didn't cancel out
+ // and set to a suitably small value. The only test case we've been able to
+ // find for this, is when v = 8.5 and x = 4*PI.
+ //
+ if(gamma == 0)
+ {
+ gamma = u * tools::epsilon<T>() / x;
+ }
 Ju = sign(current) * sqrt(W / (q + gamma * (p - t)));

 Jv = Ju * ratio; // normalization

Modified: trunk/libs/math/doc/sf_and_dist/html/index.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/index.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/index.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -527,7 +527,7 @@

</div>
<table xmlns:rev="http://www.cs.rpi.edu/~gregod/boost/tools/doc/revision" width="100%"><tr>
-<td align="left">Last revised: January 27, 2012 at 18:59:48 GMT</td>
+<td align="left">Last revised: March 24, 2012 at 11:48:31 GMT</td>
<td align="right"><div class="copyright-footer"></div></td>
</tr></table>
<hr>

Modified: trunk/libs/math/doc/sf_and_dist/html/index/s13.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/index/s13.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/index/s13.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -22,9 +22,9 @@
<div class="spirit-nav">
<a accesskey="p" href="../math_toolkit/status/credits.html"><img src="../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../index.html"><img src="../../../../../../doc/src/images/up.png" alt="Up"></a><a accesskey="h" href="../index.html"><img src="../../../../../../doc/src/images/home.png" alt="Home"></a><a accesskey="n" href="s14.html"><img src="../../../../../../doc/src/images/next.png" alt="Next"></a>
</div>
-<div class="section id1486091">
+<div class="section id1487367">
<div class="titlepage"><div><div><h2 class="title" style="clear: both">
-<a name="id1486091"></a>Function Index</h2></div></div></div>
+<a name="id1487367"></a>Function Index</h2></div></div></div>
<a class="link" href="s13.html#idx_id_0">A</a> <a class="link" href="s13.html#idx_id_1">B</a> <a class="link" href="s13.html#idx_id_2">C</a> <a class="link" href="s13.html#idx_id_3">D</a> <a class="link" href="s13.html#idx_id_4">E</a> <a class="link" href="s13.html#idx_id_5">F</a> <a class="link" href="s13.html#idx_id_6">G</a> <a class="link" href="s13.html#idx_id_7">H</a> <a class="link" href="s13.html#idx_id_8">I</a> <a class="link" href="s13.html#idx_id_9">K</a> <a class="link" href="s13.html#idx_id_10">L</a> <a class="link" href="s13.html#idx_id_11">M</a> <a class="link" href="s13.html#idx_id_12">N</a> <a class="link" href="s13.html#idx_id_14">P</a> <a class="link" href="s13.html#idx_id_15">Q</a> <a class="link" href="s13.html#idx_id_16">R</a> <a class="link" href="s13.html#idx_id_17">S</a> <a class="link" href="s13.html#idx_id_18">T</a> <a class="link" href="s13.html#idx_id_19">U</a> <a class="link" href="s13.html#idx_id_20">V</a> <a class="link" href="s13.html#idx_id_21">W</a> <a class="link" href=
"s13.html#idx_id_22">Z</a>
<div class="variablelist"><dl>
<dt>

Modified: trunk/libs/math/doc/sf_and_dist/html/index/s14.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/index/s14.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/index/s14.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -22,9 +22,9 @@
<div class="spirit-nav">
<a accesskey="p" href="s13.html"><img src="../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../index.html"><img src="../../../../../../doc/src/images/up.png" alt="Up"></a><a accesskey="h" href="../index.html"><img src="../../../../../../doc/src/images/home.png" alt="Home"></a><a accesskey="n" href="s15.html"><img src="../../../../../../doc/src/images/next.png" alt="Next"></a>
</div>
-<div class="section id1501215">
+<div class="section id1502492">
<div class="titlepage"><div><div><h2 class="title" style="clear: both">
-<a name="id1501215"></a>Class Index</h2></div></div></div>
+<a name="id1502492"></a>Class Index</h2></div></div></div>
<a class="link" href="s14.html#idx_id_24">B</a> <a class="link" href="s14.html#idx_id_25">C</a> <a class="link" href="s14.html#idx_id_26">D</a> <a class="link" href="s14.html#idx_id_27">E</a> <a class="link" href="s14.html#idx_id_28">F</a> <a class="link" href="s14.html#idx_id_29">G</a> <a class="link" href="s14.html#idx_id_30">H</a> <a class="link" href="s14.html#idx_id_31">I</a> <a class="link" href="s14.html#idx_id_33">L</a> <a class="link" href="s14.html#idx_id_34">M</a> <a class="link" href="s14.html#idx_id_35">N</a> <a class="link" href="s14.html#idx_id_37">P</a> <a class="link" href="s14.html#idx_id_39">R</a> <a class="link" href="s14.html#idx_id_40">S</a> <a class="link" href="s14.html#idx_id_41">T</a> <a class="link" href="s14.html#idx_id_42">U</a> <a class="link" href="s14.html#idx_id_44">W</a>
<div class="variablelist"><dl>
<dt>

Modified: trunk/libs/math/doc/sf_and_dist/html/index/s15.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/index/s15.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/index/s15.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -22,9 +22,9 @@
<div class="spirit-nav">
<a accesskey="p" href="s14.html"><img src="../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../index.html"><img src="../../../../../../doc/src/images/up.png" alt="Up"></a><a accesskey="h" href="../index.html"><img src="../../../../../../doc/src/images/home.png" alt="Home"></a><a accesskey="n" href="s16.html"><img src="../../../../../../doc/src/images/next.png" alt="Next"></a>
</div>
-<div class="section id1502148">
+<div class="section id1503424">
<div class="titlepage"><div><div><h2 class="title" style="clear: both">
-<a name="id1502148"></a>Typedef Index</h2></div></div></div>
+<a name="id1503424"></a>Typedef Index</h2></div></div></div>
<a class="link" href="s15.html#idx_id_46">A</a> <a class="link" href="s15.html#idx_id_47">B</a> <a class="link" href="s15.html#idx_id_48">C</a> <a class="link" href="s15.html#idx_id_49">D</a> <a class="link" href="s15.html#idx_id_50">E</a> <a class="link" href="s15.html#idx_id_51">F</a> <a class="link" href="s15.html#idx_id_52">G</a> <a class="link" href="s15.html#idx_id_53">H</a> <a class="link" href="s15.html#idx_id_54">I</a> <a class="link" href="s15.html#idx_id_56">L</a> <a class="link" href="s15.html#idx_id_58">N</a> <a class="link" href="s15.html#idx_id_59">O</a> <a class="link" href="s15.html#idx_id_60">P</a> <a class="link" href="s15.html#idx_id_62">R</a> <a class="link" href="s15.html#idx_id_63">S</a> <a class="link" href="s15.html#idx_id_64">T</a> <a class="link" href="s15.html#idx_id_65">U</a> <a class="link" href="s15.html#idx_id_66">V</a> <a class="link" href="s15.html#idx_id_67">W</a>
<div class="variablelist"><dl>
<dt>

Modified: trunk/libs/math/doc/sf_and_dist/html/index/s16.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/index/s16.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/index/s16.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -22,9 +22,9 @@
<div class="spirit-nav">
<a accesskey="p" href="s15.html"><img src="../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../index.html"><img src="../../../../../../doc/src/images/up.png" alt="Up"></a><a accesskey="h" href="../index.html"><img src="../../../../../../doc/src/images/home.png" alt="Home"></a><a accesskey="n" href="s17.html"><img src="../../../../../../doc/src/images/next.png" alt="Next"></a>
</div>
-<div class="section id1506551">
+<div class="section id1506991">
<div class="titlepage"><div><div><h2 class="title" style="clear: both">
-<a name="id1506551"></a>Macro Index</h2></div></div></div>
+<a name="id1506991"></a>Macro Index</h2></div></div></div>
<a class="link" href="s16.html#idx_id_70">B</a> <a class="link" href="s16.html#idx_id_74">F</a>
<div class="variablelist"><dl>
<dt>

Modified: trunk/libs/math/doc/sf_and_dist/html/index/s17.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/index/s17.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/index/s17.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -21,9 +21,9 @@
<div class="spirit-nav">
<a accesskey="p" href="s16.html"><img src="../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../index.html"><img src="../../../../../../doc/src/images/up.png" alt="Up"></a><a accesskey="h" href="../index.html"><img src="../../../../../../doc/src/images/home.png" alt="Home"></a>
</div>
-<div class="section id1507777">
+<div class="section id1508218">
<div class="titlepage"><div><div><h2 class="title" style="clear: both">
-<a name="id1507777"></a>Index</h2></div></div></div>
+<a name="id1508218"></a>Index</h2></div></div></div>
<a class="link" href="s17.html#idx_id_92">A</a> <a class="link" href="s17.html#idx_id_93">B</a> <a class="link" href="s17.html#idx_id_94">C</a> <a class="link" href="s17.html#idx_id_95">D</a> <a class="link" href="s17.html#idx_id_96">E</a> <a class="link" href="s17.html#idx_id_97">F</a> <a class="link" href="s17.html#idx_id_98">G</a> <a class="link" href="s17.html#idx_id_99">H</a> <a class="link" href="s17.html#idx_id_100">I</a> <a class="link" href="s17.html#idx_id_101">K</a> <a class="link" href="s17.html#idx_id_102">L</a> <a class="link" href="s17.html#idx_id_103">M</a> <a class="link" href="s17.html#idx_id_104">N</a> <a class="link" href="s17.html#idx_id_105">O</a> <a class="link" href="s17.html#idx_id_106">P</a> <a class="link" href="s17.html#idx_id_107">Q</a> <a class="link" href="s17.html#idx_id_108">R</a> <a class="link" href="s17.html#idx_id_109">S</a> <a class="link" href="s17.html#idx_id_110">T</a> <a class="link" href="s17.html#idx_id_111">U</a> <a class="link" href="s17.html#idx_id_112">V</a
> <a class="link" href="s17.html#idx_id_113">W</a> <a class="link" href="s17.html#idx_id_114">Z</a>
<div class="variablelist"><dl>
<dt>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/constants/new_const.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/constants/new_const.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/constants/new_const.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -232,7 +232,6 @@
 } // namespace constants
 } // namespace math
} // namespace boost
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/binom_eg/binomial_coinflip_example.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/binom_eg/binomial_coinflip_example.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/binom_eg/binomial_coinflip_example.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -66,7 +66,8 @@
{
 cout << "Using Binomial distribution to predict how many heads and tails." << endl;
 try
- {</pre>
+ {
+</pre>

 

@@ -83,7 +84,8 @@
int flips = 10;
binomial flip(flips, success_fraction);

-cout.precision(4);</pre>
+cout.precision(4);
+</pre>

 

@@ -103,7 +105,8 @@
// if success_fraction is exactly one half,
// for example, when flipping 'fair' coins.
cout << "Skewness if success_fraction is " << flip.success_fraction()
- << " is " << skewness(flip) << endl << endl; // Expect zero for a 'fair' coin.</pre>
+ << " is " << skewness(flip) << endl << endl; // Expect zero for a 'fair' coin.
+</pre>

 

@@ -113,7 +116,8 @@

<pre class="programlisting">cout << "For " << flip.trials() << " coin flips: " << endl;
cout << "Probability of getting no heads is " << pdf(flip, 0) << endl;
-cout << "Probability of getting at least one head is " << 1. - pdf(flip, 0) << endl;</pre>
+cout << "Probability of getting at least one head is " << 1. - pdf(flip, 0) << endl;
+</pre>

 

@@ -123,7 +127,8 @@


<pre class="programlisting">cout << "Probability of getting 0 or 1 heads is "
- << pdf(flip, 0) + pdf(flip, 1) << endl; // sum of exactly == probabilities</pre>
+ << pdf(flip, 0) + pdf(flip, 1) << endl; // sum of exactly == probabilities
+</pre>

 

@@ -132,7 +137,8 @@


<pre class="programlisting">cout << "Probability of getting 0 or 1 (<= 1) heads is " << cdf(flip, 1) << endl;
-cout << "Probability of getting 9 or 10 heads is " << pdf(flip, 9) + pdf(flip, 10) << endl;</pre>
+cout << "Probability of getting 9 or 10 heads is " << pdf(flip, 9) + pdf(flip, 10) << endl;
+</pre>

 

@@ -140,7 +146,8 @@
 


-<pre class="programlisting">cout << "Probability of getting 9 or 10 heads is " << 1. - cdf(flip, 8) << endl;</pre>
+<pre class="programlisting">cout << "Probability of getting 9 or 10 heads is " << 1. - cdf(flip, 8) << endl;
+</pre>

 

@@ -148,7 +155,8 @@
 


-<pre class="programlisting">cout << "Probability of getting 9 or 10 heads is " << cdf(complement(flip, 8)) << endl;</pre>
+<pre class="programlisting">cout << "Probability of getting 9 or 10 heads is " << cdf(complement(flip, 8)) << endl;
+</pre>

 

@@ -165,7 +173,8 @@

<pre class="programlisting">cout << "Probability of between 4 and 6 heads (4 or 5 or 6) is "
 // P(X == 4) + P(X == 5) + P(X == 6)
- << pdf(flip, 4) + pdf(flip, 5) + pdf(flip, 6) << endl;</pre>
+ << pdf(flip, 4) + pdf(flip, 5) + pdf(flip, 6) << endl;
+</pre>

 

@@ -175,7 +184,8 @@

<pre class="programlisting">cout << "Probability of between 4 and 6 heads (4 or 5 or 6) is "
 // P(X <= 6) - P(X <= 3) == P(X < 4)
- << cdf(flip, 6) - cdf(flip, 3) << endl;</pre>
+ << cdf(flip, 6) - cdf(flip, 3) << endl;
+</pre>

 

@@ -186,7 +196,8 @@
<pre class="programlisting">cout << "Probability of between 3 and 7 heads (3, 4, 5, 6 or 7) is "
 // P(X <= 7) - P(X <= 2) == P(X < 3)
 << cdf(flip, 7) - cdf(flip, 2) << endl;
-cout << endl;</pre>
+cout << endl;
+</pre>

 

@@ -215,7 +226,8 @@
 cout << setw(2) << successes << " " << setw(10) << left
 << probability << " or 1 in " << 1. / probability << ", or "
 << probability * 100. << "%"<< endl;
-} // for i</pre>
+} // for i
+</pre>

 

@@ -226,7 +238,8 @@
<pre class="programlisting">}
catch(const std::exception& e)
{
- //</pre>
+ //
+</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/binom_eg/binomial_quiz_example.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/binom_eg/binomial_quiz_example.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/binom_eg/binomial_quiz_example.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -55,8 +55,6 @@
 using std::setw; using std::setprecision;
#include <exception>
 using std::exception;
-
-
</pre>

 
@@ -70,7 +68,8 @@
<pre class="programlisting">int questions = 16; // All the questions in the quiz.
int answers = 4; // Possible answers to each question.
double success_fraction = 1. / answers; // If a random guess, p = 1/4 = 0.25.
-binomial quiz(questions, success_fraction);</pre>
+binomial quiz(questions, success_fraction);
+</pre>

 

@@ -81,7 +80,8 @@
<pre class="programlisting">cout << "In a quiz with " << quiz.trials()
 << " questions and with a probability of guessing right of "
 << quiz.success_fraction() * 100 << " %"
- << " or 1 in " << static_cast<int>(1. / quiz.success_fraction()) << endl;</pre>
+ << " or 1 in " << static_cast<int>(1. / quiz.success_fraction()) << endl;
+</pre>

 

@@ -96,7 +96,8 @@
cout << "Probability of getting exactly " << pass_score << " answers right by chance is "
 << pdf(quiz, pass_score) << endl;
cout << "Probability of getting all " << questions << " answers right by chance is "
- << pdf(quiz, questions) << endl;</pre>
+ << pdf(quiz, questions) << endl;
+</pre>

 
<pre class="programlisting">Probability of getting none right is 0.0100226
@@ -119,7 +120,8 @@
 double probability = pdf(quiz, successes);
 cout << setw(2) << successes << " " << probability << endl;
}
-cout << endl;</pre>
+cout << endl;
+</pre>

 
<pre class="programlisting">Guessed Probability
@@ -146,7 +148,8 @@
 


-<pre class="programlisting">cout << "Probability of getting none or one right is " << pdf(quiz, 0) + pdf(quiz, 1) << endl;</pre>
+<pre class="programlisting">cout << "Probability of getting none or one right is " << pdf(quiz, 0) + pdf(quiz, 1) << endl;
+</pre>

 
<pre class="programlisting">Probability of getting none or one right is 0.0634764
@@ -158,7 +161,8 @@
 


-<pre class="programlisting">cout << "Probability of getting none or one right is " << cdf(quiz, 1) << endl;</pre>
+<pre class="programlisting">cout << "Probability of getting none or one right is " << cdf(quiz, 1) << endl;
+</pre>

 
<pre class="programlisting">Probability of getting none or one right is 0.0634764
@@ -169,7 +173,8 @@
 


-<pre class="programlisting">cout << "Probability of getting <= 10 right (to fail) is " << cdf(quiz, 10) << endl;</pre>
+<pre class="programlisting">cout << "Probability of getting <= 10 right (to fail) is " << cdf(quiz, 10) << endl;
+</pre>

 
<pre class="programlisting">Probability of getting <= 10 right (to fail) is 0.999715
@@ -184,7 +189,8 @@
 


-<pre class="programlisting">cout << "Probability of getting > 10 right (to pass) is " << 1 - cdf(quiz, 10) << endl;</pre>
+<pre class="programlisting">cout << "Probability of getting > 10 right (to pass) is " << 1 - cdf(quiz, 10) << endl;
+</pre>

 
<pre class="programlisting">Probability of getting > 10 right (to pass) is 0.000285239
@@ -195,7 +201,8 @@
 


-<pre class="programlisting">cout << "Probability of getting > 10 right (to pass) is " << cdf(complement(quiz, 10)) << endl;</pre>
+<pre class="programlisting">cout << "Probability of getting > 10 right (to pass) is " << cdf(complement(quiz, 10)) << endl;
+</pre>

 
<pre class="programlisting">Probability of getting > 10 right (to pass) is 0.000285239
@@ -205,7 +212,8 @@
 


-<pre class="programlisting">BOOST_ASSERT((cdf(quiz, 10) + cdf(complement(quiz, 10))) == 1.);</pre>
+<pre class="programlisting">BOOST_ASSERT((cdf(quiz, 10) + cdf(complement(quiz, 10))) == 1.);
+</pre>

 

@@ -216,7 +224,8 @@

<pre class="programlisting">cout << "Probability of getting less than " << pass_score
 << " (< " << pass_score << ") answers right by guessing is "
- << cdf(quiz, pass_score -1) << endl;</pre>
+ << cdf(quiz, pass_score -1) << endl;
+</pre>

 
<pre class="programlisting">Probability of getting less than 11 (< 11) answers right by guessing is 0.999715
@@ -235,7 +244,8 @@
 << cdf(complement(quiz, pass_score-1))
 << ", only 1 in " << 1/cdf(complement(quiz, pass_score-1)) << endl;

-BOOST_ASSERT((cdf(quiz, pass_score -1) + cdf(complement(quiz, pass_score-1))) == 1);</pre>
+BOOST_ASSERT((cdf(quiz, pass_score -1) + cdf(complement(quiz, pass_score-1))) == 1);
+</pre>

 
<pre class="programlisting">Probability of getting at least 11 (>= 11) answers right by guessing is 0.000285239, only 1 in 3505.83
@@ -251,7 +261,8 @@
 cout << setw(2) << score << " " << setprecision(10)
 << cdf(quiz, score) << endl;
}
-cout << endl;</pre>
+cout << endl;
+</pre>

 
<pre class="programlisting">At most (<=)
@@ -281,7 +292,8 @@
{
 cout << setw(2) << score << " " << setprecision(10)
 << cdf(complement(quiz, score)) << endl;
-}</pre>
+}
+</pre>

 
<pre class="programlisting">At least (>)
@@ -324,7 +336,8 @@
cout.precision(4);
cout << "Probability of getting between "
 << low << " and " << high << " answers right by guessing is "
- << sum << endl; // 0.61323</pre>
+ << sum << endl; // 0.61323
+</pre>

 
<pre class="programlisting">Probability of getting between 3 and 5 answers right by guessing is 0.6132
@@ -335,7 +348,8 @@


<pre class="programlisting">cout << "Probability of getting between " << low << " and " << high << " answers right by guessing is "
- << cdf(quiz, high) - cdf(quiz, low - 1) << endl; // 0.61323</pre>
+ << cdf(quiz, high) - cdf(quiz, low - 1) << endl; // 0.61323
+</pre>

 
<pre class="programlisting">Probability of getting between 3 and 5 answers right by guessing is 0.6132
@@ -353,7 +367,8 @@
 << cdf(quiz, high) - cdf(quiz, low - 1) << endl; // 1 <= x 8 P = 0.9825
low = 4; high = 4;
cout << "Probability of getting between " << low << " and " << high << " answers right by guessing is "
- << cdf(quiz, high) - cdf(quiz, low - 1) << endl; // 4 <= x 4 P = 0.22520</pre>
+ << cdf(quiz, high) - cdf(quiz, low - 1) << endl; // 4 <= x 4 P = 0.22520
+</pre>

 
<pre class="programlisting">Probability of getting between 1 and 6 answers right by guessing is 0.9104
@@ -377,7 +392,8 @@
 << ceil(mean(quiz) - standard_deviation(quiz)) << " and "
 << floor(mean(quiz) + standard_deviation(quiz)) << " correct." << endl;
cout << "Mode (the most frequent) is " << mode(quiz) << endl;
-cout << "Skewness is " << skewness(quiz) << endl;</pre>
+cout << "Skewness is " << skewness(quiz) << endl;
+</pre>

 
<pre class="programlisting">By guessing, on average, one can expect to get 4 correct answers.
@@ -410,7 +426,8 @@
 << quantile(complement(quiz, 0.02)) << endl; // 2 to 98%

cout << "If guessing then percentiles 1 to 99% will get " << quantile(quiz, 0.01)
- << " to " << quantile(complement(quiz, 0.01)) << " right." << endl;</pre>
+ << " to " << quantile(complement(quiz, 0.01)) << " right." << endl;
+</pre>

 

@@ -457,7 +474,8 @@
using boost::math::policies::discrete_quantile;
using boost::math::policies::real;
using boost::math::policies::integer_round_outwards; // Default.
-typedef boost::math::policies::policy<discrete_quantile<real> > real_quantile_policy;</pre>
+typedef boost::math::policies::policy<discrete_quantile<real> > real_quantile_policy;
+</pre>

 

@@ -473,7 +491,8 @@


<pre class="programlisting">using boost::math::binomial_distribution;
-typedef binomial_distribution<double, real_quantile_policy> real_quantile_binomial;</pre>
+typedef binomial_distribution<double, real_quantile_policy> real_quantile_binomial;
+</pre>

 

@@ -481,7 +500,8 @@
 


-<pre class="programlisting">real_quantile_binomial quiz_real(questions, success_fraction);</pre>
+<pre class="programlisting">real_quantile_binomial quiz_real(questions, success_fraction);
+</pre>

 

@@ -504,7 +524,8 @@
 << quantile(complement(quiz_real, 0.02)) << endl; // 2 to 98% 0.31311 7.7880

cout << "If guessing, then percentiles 1 to 99% will get " << quantile(quiz_real, 0.01)
- << " to " << quantile(complement(quiz_real, 0.01)) << " right." << endl;</pre>
+ << " to " << quantile(complement(quiz_real, 0.01)) << " right." << endl;
+</pre>

 
<pre class="programlisting">Real Quantiles

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/dist_construct_eg.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/dist_construct_eg.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/dist_construct_eg.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -52,7 +52,8 @@
#include <boost/math/distributions/binomial.hpp> // for binomial_distribution.
#include <boost/math/distributions/beta.hpp> // for beta_distribution.
#include <boost/math/distributions/gamma.hpp> // for gamma_distribution.
-#include <boost/math/distributions/normal.hpp> // for normal_distribution.</pre>
+#include <boost/math/distributions/normal.hpp> // for normal_distribution.
+</pre>

 

@@ -64,7 +65,8 @@
 


-<pre class="programlisting">boost::math::negative_binomial_distribution<double> mydist0(8., 0.25);</pre>
+<pre class="programlisting">boost::math::negative_binomial_distribution<double> mydist0(8., 0.25);
+</pre>

 

@@ -72,7 +74,8 @@
 


-<pre class="programlisting">using namespace boost::math;</pre>
+<pre class="programlisting">using namespace boost::math;
+</pre>

 

@@ -81,7 +84,8 @@
 


-<pre class="programlisting">using boost::math::negative_binomial_distribution;</pre>
+<pre class="programlisting">using boost::math::negative_binomial_distribution;
+</pre>

 

@@ -94,7 +98,8 @@
 


-<pre class="programlisting">negative_binomial_distribution<> mydist9(8., 0.25); // Uses default RealType = double.</pre>
+<pre class="programlisting">negative_binomial_distribution<> mydist9(8., 0.25); // Uses default RealType = double.
+</pre>

 

@@ -120,7 +125,8 @@
 


-<pre class="programlisting">using boost::math::negative_binomial;</pre>
+<pre class="programlisting">using boost::math::negative_binomial;
+</pre>

 

@@ -128,7 +134,8 @@
 


-<pre class="programlisting">negative_binomial mydist(8., 0.25);</pre>
+<pre class="programlisting">negative_binomial mydist(8., 0.25);
+</pre>

 

@@ -136,7 +143,8 @@
 


-<pre class="programlisting">negative_binomial mydist10(5., 0.4); // Both arguments double.</pre>
+<pre class="programlisting">negative_binomial mydist10(5., 0.4); // Both arguments double.
+</pre>

 

@@ -144,7 +152,8 @@
 


-<pre class="programlisting">negative_binomial mydist11(5, 0.4); // Using provided typedef double, int and double arguments.</pre>
+<pre class="programlisting">negative_binomial mydist11(5, 0.4); // Using provided typedef double, int and double arguments.
+</pre>

 

@@ -153,7 +162,8 @@


<pre class="programlisting">negative_binomial mydist12(5., 0.4F); // Double and float arguments.
-negative_binomial mydist13(5, 1); // Both arguments integer.</pre>
+negative_binomial mydist13(5, 1); // Both arguments integer.
+</pre>

 

@@ -162,7 +172,8 @@


<pre class="programlisting">binomial mybinomial(1, 0.5); // is more concise than
-binomial_distribution<> mybinomd1(1, 0.5);</pre>
+binomial_distribution<> mybinomd1(1, 0.5);
+</pre>

 

@@ -187,7 +198,8 @@


<pre class="programlisting">using boost::math::beta_distribution;
-beta_distribution<> mybetad1(1, 0.5);</pre>
+beta_distribution<> mybetad1(1, 0.5);
+</pre>

 

@@ -195,7 +207,8 @@
 


-<pre class="programlisting">gamma_distribution<> mygammad1(1, 0.5);</pre>
+<pre class="programlisting">gamma_distribution<> mygammad1(1, 0.5);
+</pre>

 

@@ -222,7 +235,8 @@
negative_binomial_distribution<double> mydist6(8., 0.25);

// Explicit long double precision:
-negative_binomial_distribution<long double> mydist7(8., 0.25);</pre>
+negative_binomial_distribution<long double> mydist7(8., 0.25);
+</pre>

 

@@ -265,7 +279,8 @@
 normal norm3(2, 3); // Mean = 2, std deviation = 3.

 return 0;
-} // int main()</pre>
+} // int main()
+</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/error_eg.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/error_eg.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/error_eg.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -93,7 +93,8 @@
 using std::exception;

#include <cstddef>
- // using ::errno</pre>
+ // using ::errno
+</pre>

 

@@ -135,7 +136,8 @@
 "\n""Message from thrown exception was:\n " << e.what() << std::endl;
 }
 return 0;
-} // int main()</pre>
+} // int main()
+</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/find_eg/find_location_eg.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/find_eg/find_location_eg.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/find_eg/find_location_eg.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -50,7 +50,6 @@
 using std::setw; using std::setprecision;
#include <limits>
 using std::numeric_limits;
-
</pre>

 
@@ -62,7 +61,8 @@


<pre class="programlisting">normal N01; // Default 'standard' normal distribution with zero mean and 
-double sd = 1.; // normal default standard deviation is 1.</pre>
+double sd = 1.; // normal default standard deviation is 1.
+</pre>

 

@@ -82,7 +82,8 @@
cout << "Normal distribution with mean = " << N01.location()
 << ", standard deviation " << N01.scale()
 << ", has " << "fraction > " << z
- << ", p = " << cdf(complement(N01, z)) << endl; // Note: uses complement.</pre>
+ << ", p = " << cdf(complement(N01, z)) << endl; // Note: uses complement.
+</pre>

 
<pre class="programlisting">Normal distribution with mean = 0, standard deviation 1, has fraction <= -2, p = 0.0227501
@@ -94,7 +95,8 @@


<pre class="programlisting">double l = find_location<normal>(z, p, sd);
-cout << "offset location (mean) = " << l << endl;</pre>
+cout << "offset location (mean) = " << l << endl;
+</pre>

 

@@ -112,7 +114,8 @@
 


-<pre class="programlisting">normal np001pc(l, sd); // Same standard_deviation (scale) but with mean (location) shifted.</pre>
+<pre class="programlisting">normal np001pc(l, sd); // Same standard_deviation (scale) but with mean (location) shifted.
+</pre>

 

@@ -125,7 +128,8 @@
 << ", p = " << cdf(np001pc, z) << endl;
 cout << "Normal distribution with mean = " << l
 << " has " << "fraction > " << z
- << ", p = " << cdf(complement(np001pc, z)) << endl;</pre>
+ << ", p = " << cdf(complement(np001pc, z)) << endl;
+</pre>

 
<pre class="programlisting">Normal distribution with mean = 1.09023 has fraction <= -2, p = 0.001
@@ -148,7 +152,8 @@

<pre class="programlisting">using boost::math::policies::policy;
using boost::math::policies::domain_error;
-using boost::math::policies::ignore_error;</pre>
+using boost::math::policies::ignore_error;
+</pre>

 

@@ -165,7 +170,8 @@
// Default policy, fully specified.
l = find_location<normal>(z, p, sd, boost::math::policies::policy<>());
// A new policy, ignoring domain errors, without using a typedef.
-l = find_location<normal>(z, p, sd, policy<domain_error<ignore_error> >());</pre>
+l = find_location<normal>(z, p, sd, policy<domain_error<ignore_error> >());
+</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/find_eg/find_mean_and_sd_eg.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/find_eg/find_mean_and_sd_eg.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/find_eg/find_mean_and_sd_eg.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -77,7 +77,8 @@

<pre class="programlisting">double mean = 3.; // kg
double standard_deviation = 0.1; // kg
-normal packs(mean, standard_deviation);</pre>
+normal packs(mean, standard_deviation);
+</pre>

 

@@ -88,7 +89,8 @@

<pre class="programlisting">double max_weight = 3.1; // kg
cout << "Percentage of packs > " << max_weight << " is "
-<< cdf(complement(packs, max_weight)) * 100. << endl; // P(X > 3.1)</pre>
+<< cdf(complement(packs, max_weight)) * 100. << endl; // P(X > 3.1)
+</pre>

 

@@ -105,7 +107,8 @@
<pre class="programlisting">double minimum_weight = 2.9;
cout <<"Fraction of packs <= " << minimum_weight << " with a mean of " << mean
 << " is " << cdf(complement(packs, minimum_weight)) << endl;
-// fraction of packs <= 2.9 with a mean of 3 is 0.841345</pre>
+// fraction of packs <= 2.9 with a mean of 3 is 0.841345
+</pre>

 

@@ -128,7 +131,8 @@
cout << "Fraction of packs >= " << minimum_weight
<< " with a mean of " << xpacks.mean()
 << " is " << cdf(complement(xpacks, minimum_weight)) << endl;
-// fraction of packs >= 2.9 with a mean of 3.06449 is 0.950005</pre>
+// fraction of packs >= 2.9 with a mean of 3.06449 is 0.950005
+</pre>

 

@@ -147,7 +151,8 @@
cout << "Setting the packer to " << nominal_mean << " will mean that "
 << "fraction of packs >= " << minimum_weight
 << " is " << cdf(complement(nominal_packs, minimum_weight)) << endl;
-// Setting the packer to 3.06449 will mean that fraction of packs >= 2.9 is 0.95</pre>
+// Setting the packer to 3.06449 will mean that fraction of packs >= 2.9 is 0.95
+</pre>

 

@@ -159,7 +164,8 @@


<pre class="programlisting">#include <boost/math/distributions/find_location.hpp>
- using boost::math::find_location;</pre>
+ using boost::math::find_location;
+</pre>

 

@@ -169,7 +175,8 @@


<pre class="programlisting">double safe_mean = find_location<normal>(minimum_weight, under_fraction, standard_deviation);
-normal good_packs(safe_mean, standard_deviation);</pre>
+normal good_packs(safe_mean, standard_deviation);
+</pre>

 

@@ -180,7 +187,8 @@
<pre class="programlisting">cout << "Setting the packer to " << nominal_mean << " will mean that "
 << "fraction of packs >= " << minimum_weight
 << " is " << cdf(complement(good_packs, minimum_weight)) << endl;
-// Setting the packer to 3.06449 will mean that fraction of packs >= 2.9 is 0.95</pre>
+// Setting the packer to 3.06449 will mean that fraction of packs >= 2.9 is 0.95
+</pre>

 
<h5>
@@ -213,7 +221,8 @@
cout << "Cauchy Setting the packer to " << mean << " will mean that "
 << "fraction of packs >= " << minimum_weight
 << " is " << cdf(complement(cpacks, minimum_weight)) << endl;
-// Cauchy Setting the packer to 3 will mean that fraction of packs >= 2.9 is 0.75</pre>
+// Cauchy Setting the packer to 3 will mean that fraction of packs >= 2.9 is 0.75
+</pre>

 

@@ -225,7 +234,8 @@

<pre class="programlisting">double lc = find_location<cauchy>(minimum_weight, under_fraction, standard_deviation);
cout << "find_location<cauchy>(minimum_weight, over fraction, standard_deviation); " << lc << endl;
-// find_location<cauchy>(minimum_weight, over fraction, packs.standard_deviation()); 3.53138</pre>
+// find_location<cauchy>(minimum_weight, over fraction, packs.standard_deviation()); 3.53138
+</pre>

 

@@ -241,7 +251,8 @@
cout << "Cauchy Setting the packer to " << lc << " will mean that "
 << "fraction of packs >= " << minimum_weight
 << " is " << cdf(complement(goodcpacks, minimum_weight)) << endl;
-// Cauchy Setting the packer to 3.53138 will mean that fraction of packs >= 2.9 is 0.95</pre>
+// Cauchy Setting the packer to 3.53138 will mean that fraction of packs >= 2.9 is 0.95
+</pre>

 

@@ -254,7 +265,8 @@
 << find_location<cauchy>(minimum_weight, 0.99, standard_deviation)
 << " will mean that "
 << "fraction of packs >= " << minimum_weight
- << " is " << cdf(complement(goodcpacks, minimum_weight)) << endl;</pre>
+ << " is " << cdf(complement(goodcpacks, minimum_weight)) << endl;
+</pre>

 

@@ -313,7 +325,8 @@

<pre class="programlisting">double p = 0.05; // wanted p th quantile.
cout << "Quantile of " << p << " = " << quantile(packs, p)
- << ", mean = " << packs.mean() << ", sd = " << packs.standard_deviation() << endl;</pre>
+ << ", mean = " << packs.mean() << ", sd = " << packs.standard_deviation() << endl;
+</pre>

 

@@ -338,7 +351,8 @@
cout <<"Fraction of packs >= " << minimum_weight << " with a mean of " << mean
 << " and standard deviation of " << pack05.standard_deviation()
 << " is " << cdf(complement(pack05, minimum_weight)) << endl;
-// Fraction of packs >= 2.9 with a mean of 3 and standard deviation of 0.05 is 0.97725</pre>
+// Fraction of packs >= 2.9 with a mean of 3 and standard deviation of 0.05 is 0.97725
+</pre>

 

@@ -357,7 +371,8 @@
cout <<"Fraction of packs >= " << minimum_weight << " with a mean of " << mean
 << " and standard deviation of " << pack06.standard_deviation()
 << " is " << cdf(complement(pack06, minimum_weight)) << endl;
-// Fraction of packs >= 2.9 with a mean of 3 and standard deviation of 0.06 is 0.95221</pre>
+// Fraction of packs >= 2.9 with a mean of 3 and standard deviation of 0.06 is 0.95221
+</pre>

 

@@ -394,7 +409,8 @@

cout << "For the "<< p << "th quantile to be located at "
 << minimum_weight << ", would need a standard deviation of " << sd95 << endl;
-// For the 0.05th quantile to be located at 2.9, would need a standard deviation of 0.0607957</pre>
+// For the 0.05th quantile to be located at 2.9, would need a standard deviation of 0.0607957
+</pre>

 

@@ -407,7 +423,8 @@
cout <<"Fraction of packs >= " << minimum_weight << " with a mean of " << mean
 << " and standard deviation of " << pack95.standard_deviation()
 << " is " << cdf(complement(pack95, minimum_weight)) << endl;
-// Fraction of packs >= 2.9 with a mean of 3 and standard deviation of 0.0607957 is 0.95</pre>
+// Fraction of packs >= 2.9 with a mean of 3 and standard deviation of 0.0607957 is 0.95
+</pre>

 

@@ -418,7 +435,8 @@

<pre class="programlisting">double ss = find_scale<normal>(minimum_weight, under_fraction, packs.mean());
cout << "find_scale<normal>(minimum_weight, under_fraction, packs.mean()); " << ss << endl;
-// find_scale<normal>(minimum_weight, under_fraction, packs.mean()); 0.0607957</pre>
+// find_scale<normal>(minimum_weight, under_fraction, packs.mean()); 0.0607957
+</pre>

 

@@ -426,7 +444,8 @@
 


-<pre class="programlisting">double over_fraction = 0.95;</pre>
+<pre class="programlisting">double over_fraction = 0.95;
+</pre>

 

@@ -450,7 +469,8 @@

<pre class="programlisting">double ss1o = find_scale<normal>(minimum_weight, 1 - over_fraction, packs.mean());
cout << "find_scale<normal>(minimum_weight, under_fraction, packs.mean()); " << ss1o << endl;
-// find_scale<normal>(minimum_weight, under_fraction, packs.mean()); 0.0607957</pre>
+// find_scale<normal>(minimum_weight, under_fraction, packs.mean()); 0.0607957
+</pre>

 

@@ -463,7 +483,8 @@

<pre class="programlisting">double ssc = find_scale<normal>(complement(minimum_weight, over_fraction, packs.mean()));
cout << "find_scale<normal>(complement(minimum_weight, over_fraction, packs.mean())); " << ssc << endl;
-// find_scale<normal>(complement(minimum_weight, over_fraction, packs.mean())); 0.0607957</pre>
+// find_scale<normal>(complement(minimum_weight, over_fraction, packs.mean())); 0.0607957
+</pre>

 

@@ -478,7 +499,8 @@
cout <<"Fraction of packs >= " << minimum_weight << " with a mean of " << mean
 << " and standard deviation of " << pack95c.standard_deviation()
 << " is " << cdf(complement(pack95c, minimum_weight)) << endl;
-// Fraction of packs >= 2.9 with a mean of 3 and standard deviation of 0.0607957 is 0.95</pre>
+// Fraction of packs >= 2.9 with a mean of 3 and standard deviation of 0.0607957 is 0.95
+</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/find_eg/find_scale_eg.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/find_eg/find_scale_eg.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/find_eg/find_scale_eg.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -59,7 +59,8 @@


<pre class="programlisting">normal N01; // Default 'standard' normal distribution with zero mean
-double sd = 1.; // and standard deviation is 1.</pre>
+double sd = 1.; // and standard deviation is 1.
+</pre>

 

@@ -79,7 +80,8 @@
cout << "Normal distribution with mean = " << N01.location()
 << ", standard deviation " << N01.scale()
 << ", has " << "fraction > " << z
- << ", p = " << cdf(complement(N01, z)) << endl; // Note: uses complement.</pre>
+ << ", p = " << cdf(complement(N01, z)) << endl; // Note: uses complement.
+</pre>

 
<pre class="programlisting">Normal distribution with mean = 0 has fraction <= -2, p = 0.0227501
@@ -94,7 +96,8 @@

<pre class="programlisting">double l = N01.location();
double s = find_scale<normal>(z, p, l);
-cout << "scale (standard deviation) = " << s << endl;</pre>
+cout << "scale (standard deviation) = " << s << endl;
+</pre>

 

@@ -111,7 +114,8 @@
 


-<pre class="programlisting">normal np001pc(N01.location(), s);</pre>
+<pre class="programlisting">normal np001pc(N01.location(), s);
+</pre>

 

@@ -124,7 +128,8 @@
 << ", p = " << cdf(np001pc, z) << endl;
cout << "Normal distribution with mean = " << l
 << " has " << "fraction > " << z
- << ", p = " << cdf(complement(np001pc, z)) << endl;</pre>
+ << ", p = " << cdf(complement(np001pc, z)) << endl;
+</pre>

 
<pre class="programlisting">Normal distribution with mean = 0 has fraction <= -2, p = 0.001
@@ -147,7 +152,8 @@

<pre class="programlisting">using boost::math::policies::policy;
using boost::math::policies::domain_error;
-using boost::math::policies::ignore_error;</pre>
+using boost::math::policies::ignore_error;
+</pre>

 

@@ -165,7 +171,8 @@
// Default policy, fully specified.
l = find_scale<normal>(z, p, l, boost::math::policies::policy<>());
// New policy, without typedef.
-l = find_scale<normal>(z, p, l, policy<domain_error<ignore_error> >());</pre>
+l = find_scale<normal>(z, p, l, policy<domain_error<ignore_error> >());
+</pre>

 

@@ -185,7 +192,8 @@
cout << "Normal distribution with mean = " << l << " has "
 << "fraction <= " << z << " = " << cdf(np95pc, z) << endl;
cout << "Normal distribution with mean = " << l << " has "
- << "fraction > " << z << " = " << cdf(complement(np95pc, z)) << endl;</pre>
+ << "fraction > " << z << " = " << cdf(complement(np95pc, z)) << endl;
+</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/geometric_eg.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/geometric_eg.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/geometric_eg.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -37,7 +37,8 @@


<pre class="programlisting">#define BOOST_MATH_OVERFLOW_ERROR_POLICY ignore_error
-#define BOOST_MATH_DISCRETE_QUANTILE_POLICY real</pre>
+#define BOOST_MATH_DISCRETE_QUANTILE_POLICY real
+</pre>

 
<div class="caution"><table border="0" summary="Caution">
@@ -109,7 +110,8 @@


<pre class="programlisting">double success_fraction = 1./6; // success_fraction (p) = 0.1666
-// (so failure_fraction is 1 - success_fraction = 5./6 = 1- 0.1666 = 0.8333)</pre>
+// (so failure_fraction is 1 - success_fraction = 5./6 = 1- 0.1666 = 0.8333)
+</pre>

 

@@ -126,7 +128,8 @@
 


-<pre class="programlisting">geometric g6(success_fraction); // type double by default.</pre>
+<pre class="programlisting">geometric g6(success_fraction); // type double by default.
+</pre>

 

@@ -134,7 +137,8 @@
 


-<pre class="programlisting">cout << "success fraction of a six-sided dice is " << g6.success_fraction() << endl;</pre>
+<pre class="programlisting">cout << "success fraction of a six-sided dice is " << g6.success_fraction() << endl;
+</pre>

 

@@ -144,7 +148,8 @@


<pre class="programlisting">cout << pdf(g6, 0) << endl; // 0.1667
-cout << cdf(g6, 0) << endl; // 0.1667</pre>
+cout << cdf(g6, 0) << endl; // 0.1667
+</pre>

 

@@ -154,7 +159,8 @@
 


-<pre class="programlisting">cout << pdf(g6, 1) << endl; // 0.1389</pre>
+<pre class="programlisting">cout << pdf(g6, 1) << endl; // 0.1389
+</pre>

 

@@ -163,7 +169,8 @@
 


-<pre class="programlisting">cout << "pdf(g6, 0) + pdf(g6, 1) = " << pdf(g6, 0) + pdf(g6, 1) << endl;</pre>
+<pre class="programlisting">cout << "pdf(g6, 0) + pdf(g6, 1) = " << pdf(g6, 0) + pdf(g6, 1) << endl;
+</pre>

 

@@ -172,7 +179,8 @@
 


-<pre class="programlisting">cout << "cdf(g6, 1) = " << cdf(g6, 1) << endl; // 0.3056</pre>
+<pre class="programlisting">cout << "cdf(g6, 1) = " << cdf(g6, 1) << endl; // 0.3056
+</pre>

 

@@ -181,7 +189,8 @@
 


-<pre class="programlisting">cout << "cdf(g6, 12) = " << cdf(g6, 12) << endl; // 0.9065 or 90% probability.</pre>
+<pre class="programlisting">cout << "cdf(g6, 12) = " << cdf(g6, 12) << endl; // 0.9065 or 90% probability.
+</pre>

 

@@ -190,7 +199,8 @@
 


-<pre class="programlisting">cout << "quantile(g6, 0.99) = " << quantile(g6, 0.99) << endl; // 24.26</pre>
+<pre class="programlisting">cout << "quantile(g6, 0.99) = " << quantile(g6, 0.99) << endl; // 24.26
+</pre>

 

@@ -210,7 +220,8 @@

<pre class="programlisting">negative_binomial nb(1, success_fraction);
cout << pdf(nb, 1) << endl; // 0.1389
-cout << cdf(nb, 1) << endl; // 0.3056</pre>
+cout << cdf(nb, 1) << endl; // 0.3056
+</pre>

 

@@ -219,7 +230,8 @@
 


-<pre class="programlisting">cout << "quantile(complement(g6, 1 - p)) " << quantile(complement(g6, 0.01)) << endl; // 24.26</pre>
+<pre class="programlisting">cout << "quantile(complement(g6, 1 - p)) " << quantile(complement(g6, 0.01)) << endl; // 24.26
+</pre>

 

@@ -251,7 +263,8 @@
 << max_digits10 << endl;
 cout.precision(max_digits10); // 

- cout << cdf(g05, 0.0001) << endl; // returns 0.5000346561579232, not exact 0.5.</pre>
+ cout << cdf(g05, 0.0001) << endl; // returns 0.5000346561579232, not exact 0.5.
+</pre>

 

@@ -260,7 +273,8 @@
 


-<pre class="programlisting">cout << cdf(g05, floor(0.0001)) << endl; // returns exactly 0.5</pre>
+<pre class="programlisting">cout << cdf(g05, floor(0.0001)) << endl; // returns exactly 0.5
+</pre>

 
<pre class="programlisting"><code class="computeroutput">> formatC(pgeom(0.9999999,0.5, FALSE), digits=17) [1] " 0.25"</code>
@@ -313,7 +327,8 @@
c = 0.80; // Only reasonably confident.
cout << "To be " << c * 100
 << "% confident of finding we customer with a fault, need to survey "
- << ceil(quantile(g, c)) << " customers." << endl; // 79</pre>
+ << ceil(quantile(g, c)) << " customers." << endl; // 79
+</pre>

 
<h6>
@@ -332,7 +347,8 @@


<pre class="programlisting">cout.precision(2);
-geometric gav(0.75); // Shooter averages 7.5 out of 10 in the basket.</pre>
+geometric gav(0.75); // Shooter averages 7.5 out of 10 in the basket.
+</pre>

 

@@ -341,7 +357,8 @@
 


-<pre class="programlisting">cout << "Probability of score on 1st try = " << pdf(gav, 0) << endl; // 0.75</pre>
+<pre class="programlisting">cout << "Probability of score on 1st try = " << pdf(gav, 0) << endl; // 0.75
+</pre>

 

@@ -351,7 +368,8 @@
 


-<pre class="programlisting">cout << "Probability of score on 5th try = " << pdf(gav, 4) << endl; // 0.0029</pre>
+<pre class="programlisting">cout << "Probability of score on 5th try = " << pdf(gav, 4) << endl; // 0.0029
+</pre>

 

@@ -365,7 +383,8 @@
<pre class="programlisting">geometric gbest(0.95);
cout << "Probability of score on 5th try = " << pdf(gbest, 4) << endl; // 5.9e-6
geometric gmediocre(0.50);
-cout << "Probability of score on 5th try = " << pdf(gmediocre, 4) << endl; // 0.031</pre>
+cout << "Probability of score on 5th try = " << pdf(gmediocre, 4) << endl; // 0.031
+</pre>

 

@@ -408,7 +427,8 @@
 << t << endl; // 0.00025
t = geometric::find_upper_bound_on_p(k, alpha/2);
cout << "geometric::find_upper_bound_on_p(" << int(k) << ", " << alpha/2 << ") = "
- << t << endl; // 0.037</pre>
+ << t << endl; // 0.037
+</pre>

 

@@ -433,7 +453,8 @@
 << t << endl; // 5e-005
t = geometric::find_upper_bound_on_p(k, alpha/2);
cout << "geometric::find_upper_bound_on_p(" << int(k) << ", " << alpha/2 << ") = "
- << t << endl; // 0.052</pre>
+ << t << endl; // 0.052
+</pre>

 

@@ -451,7 +472,8 @@
catch(const std::exception& e)
{ // Since we have set an overflow policy of ignore_error,
 // an overflow exception should never be thrown.
- std::cout << "\nMessage from thrown exception was:\n " << e.what() << std::endl;</pre>
+ std::cout << "\nMessage from thrown exception was:\n " << e.what() << std::endl;
+</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/inverse_chi_squared_eg.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/inverse_chi_squared_eg.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/inverse_chi_squared_eg.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -130,7 +130,6 @@
cout << " 2.5% quantile: " << quantile(prior, 0.025) << endl;
cout << " 50% quantile: " << quantile(prior, 0.5) << endl;
cout << " 97.5% quantile: " << quantile(prior, 0.975) << endl << endl;
-
</pre>

 
@@ -244,7 +243,8 @@
double posteriorScale =
 (priorDF * priorScale + (ballsSampleVariance * (ballsSampleSize - 1))) / posteriorDF;
cout << "prior scale: " << priorScale << endl;
-cout << "posterior scale: " << posteriorScale << endl;</pre>
+cout << "posterior scale: " << posteriorScale << endl;
+</pre>

 

@@ -262,6 +262,7 @@
posterior degrees-of-freedom: 120
prior scale: 25
posterior scale: 49.5
+
</pre>

 To compare the generic machinery performance with our suspect machine,
@@ -284,7 +285,6 @@
 cout << " probability variance <= 25: " << boost::math::cdf(posterior, 25.0) << endl;
 cout << " probability variance > 50: "
 << boost::math::cdf(boost::math::complement(posterior, 50.0)) << endl;
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/nccs_eg/nccs_power_eg.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/nccs_eg/nccs_power_eg.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/nccs_eg/nccs_power_eg.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -50,7 +50,8 @@
using std::setprecision;

int main()
-{</pre>
+{
+</pre>

 

@@ -64,7 +65,8 @@
{
 cout << "[[lambda]=" << lam << "]";
}
-cout << "]\n";</pre>
+cout << "]\n";
+</pre>

 

@@ -80,7 +82,8 @@
{
 cout << "[[" << n << "]";
 for(int lam = 2; lam <= 20; lam += 2)
- {</pre>
+ {
+</pre>

 

@@ -88,7 +91,8 @@
 


-<pre class="programlisting">double cs = quantile(complement(chi_squared(n), 0.05));</pre>
+<pre class="programlisting">double cs = quantile(complement(chi_squared(n), 0.05));
+</pre>

 

@@ -97,7 +101,8 @@
 


-<pre class="programlisting">double beta = cdf(complement(non_central_chi_squared(n, lam), cs));</pre>
+<pre class="programlisting">double beta = cdf(complement(non_central_chi_squared(n, lam), cs));
+</pre>

 

@@ -110,7 +115,8 @@
 cout << "]" << endl;
 }
 cout << "]" << endl;
-}</pre>
+}
+</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/neg_binom_eg/neg_binom_conf.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/neg_binom_eg/neg_binom_conf.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/neg_binom_eg/neg_binom_conf.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -67,7 +67,8 @@
using std::cout; using std::endl;
#include <iomanip>
using std::setprecision;
-using std::setw; using std::left; using std::fixed; using std::right;</pre>
+using std::setw; using std::left; using std::fixed; using std::right;
+</pre>

 

@@ -76,7 +77,8 @@
 


-<pre class="programlisting">double alpha[] = { 0.5, 0.25, 0.1, 0.05, 0.01, 0.001, 0.0001, 0.00001 };</pre>
+<pre class="programlisting">double alpha[] = { 0.5, 0.25, 0.1, 0.05, 0.01, 0.001, 0.0001, 0.00001 };
+</pre>

 

@@ -112,7 +114,8 @@
 "___________________________________________\n"
 "Confidence Lower Upper\n"
 " Value (%) Limit Limit\n"
- "___________________________________________\n";</pre>
+ "___________________________________________\n";
+</pre>

 

@@ -142,7 +145,8 @@
 cout << fixed << setprecision(5) << setw(15) << right << upper << endl;
 }
 cout << endl;
-} // void confidence_limits_on_frequency(unsigned trials, unsigned successes)</pre>
+} // void confidence_limits_on_frequency(unsigned trials, unsigned successes)
+</pre>

 

@@ -159,7 +163,6 @@

 return 0;
} // int main()
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/neg_binom_eg/neg_binom_size_eg.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/neg_binom_eg/neg_binom_size_eg.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/neg_binom_eg/neg_binom_size_eg.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -48,7 +48,8 @@
 


-<pre class="programlisting">void find_number_of_trials(double failures, double p);</pre>
+<pre class="programlisting">void find_number_of_trials(double failures, double p);
+</pre>

 

@@ -58,7 +59,8 @@
 


-<pre class="programlisting">double alpha[] = { 0.5, 0.25, 0.1, 0.05, 0.01, 0.001, 0.0001, 0.00001 };</pre>
+<pre class="programlisting">double alpha[] = { 0.5, 0.25, 0.1, 0.05, 0.01, 0.001, 0.0001, 0.00001 };
+</pre>

 

@@ -80,8 +82,7 @@
 


-<pre class="programlisting">
-void find_number_of_trials(double failures, double p)
+<pre class="programlisting">void find_number_of_trials(double failures, double p)
{
 // trials = number of trials
 // failures = number of failures before achieving required success(es).
@@ -107,7 +108,8 @@
 << endl;
 }
 cout << endl;
-} // void find_number_of_trials(double failures, double p)</pre>
+} // void find_number_of_trials(double failures, double p)
+</pre>

 

@@ -127,7 +129,6 @@

 return 0;
} // int main()
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/neg_binom_eg/negative_binomial_example1.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/neg_binom_eg/negative_binomial_example1.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/neg_binom_eg/negative_binomial_example1.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -73,7 +73,8 @@


<pre class="programlisting">#define BOOST_MATH_OVERFLOW_ERROR_POLICY ignore_error
-#define BOOST_MATH_DISCRETE_QUANTILE_POLICY real</pre>
+#define BOOST_MATH_DISCRETE_QUANTILE_POLICY real
+</pre>

 

@@ -124,7 +125,8 @@


<pre class="programlisting">try
-{</pre>
+{
+</pre>

 

@@ -134,7 +136,8 @@
 


-<pre class="programlisting">double sales_quota = 5; // Pat's sales quota - successes (r).</pre>
+<pre class="programlisting">double sales_quota = 5; // Pat's sales quota - successes (r).
+</pre>

 

@@ -143,7 +146,8 @@
 


-<pre class="programlisting">double success_fraction = 0.4; // success_fraction (p) - so failure_fraction is 0.6.</pre>
+<pre class="programlisting">double success_fraction = 0.4; // success_fraction (p) - so failure_fraction is 0.6.
+</pre>

 

@@ -160,7 +164,8 @@
 


-<pre class="programlisting">negative_binomial nb(sales_quota, success_fraction); // type double by default.</pre>
+<pre class="programlisting">negative_binomial nb(sales_quota, success_fraction); // type double by default.
+</pre>

 

@@ -178,7 +183,8 @@
cout << "With a success rate of " << nb.success_fraction()
 << ", he might expect, on average,\n"
 "to need to visit about " << success_fraction * all_houses
- << " houses in order to sell all " << nb.successes() << " bars. " << endl;</pre>
+ << " houses in order to sell all " << nb.successes() << " bars. " << endl;
+</pre>

 
<pre class="programlisting">Pat has a sales per house success rate of 0.4.
@@ -200,7 +206,8 @@


<pre class="programlisting">cout << "Probability that Pat finishes on the " << sales_quota << "th house is "
- << pdf(nb, 5 - sales_quota) << endl; // == pdf(nb, 0)</pre>
+ << pdf(nb, 5 - sales_quota) << endl; // == pdf(nb, 0)
+</pre>

 

@@ -220,7 +227,8 @@
cout << "Probability that Pat finishes on the 7th house is "
 << pdf(nb, 7 - sales_quota) << endl;
cout << "Probability that Pat finishes on the 8th house is "
- << pdf(nb, 8 - sales_quota) << endl;</pre>
+ << pdf(nb, 8 - sales_quota) << endl;
+</pre>

 
<pre class="programlisting">Probability that Pat finishes on the 6th house is 0.03072
@@ -240,7 +248,8 @@
 + pdf(nb, 6 - sales_quota) // 1 failure.
 + pdf(nb, 7 - sales_quota) // 2 failures.
 + pdf(nb, 8 - sales_quota) // 3 failures.
- << endl;</pre>
+ << endl;
+</pre>

 
<pre class="programlisting">pdf(sales_quota) + pdf(6) + pdf(7) + pdf(8) = 0.17367
@@ -253,7 +262,8 @@

<pre class="programlisting">cout << "\nProbability of selling his quota of " << sales_quota
 << " bars\non or before the " << 8 << "th house is "
- << cdf(nb, 8 - sales_quota) << endl;</pre>
+ << cdf(nb, 8 - sales_quota) << endl;
+</pre>

 
<pre class="programlisting">Probability of selling his quota of 5 bars on or before the 8th house is 0.17367
@@ -264,7 +274,8 @@
 << pdf(nb, 10 - sales_quota) << endl;
cout << "\nProbability of selling his quota of " << sales_quota
 << " bars\non or before the " << 10 << "th house is "
- << cdf(nb, 10 - sales_quota) << endl;</pre>
+ << cdf(nb, 10 - sales_quota) << endl;
+</pre>

 
<pre class="programlisting">Probability that Pat finishes exactly on the 10th house is 0.10033
@@ -276,7 +287,8 @@
 << pdf(nb, 11 - sales_quota) << endl;
cout << "\nProbability of selling his quota of " << sales_quota
 << " bars\non or before the " << 11 << "th house is "
- << cdf(nb, 11 - sales_quota) << endl;</pre>
+ << cdf(nb, 11 - sales_quota) << endl;
+</pre>

 
<pre class="programlisting">Probability that Pat finishes on the 11th house is 0.10033
@@ -290,7 +302,8 @@

cout << "\nProbability of selling his quota of " << sales_quota
 << " bars\non or before the " << 12 << "th house is "
- << cdf(nb, 12 - sales_quota) << endl;</pre>
+ << cdf(nb, 12 - sales_quota) << endl;
+</pre>

 
<pre class="programlisting">Probability that Pat finishes on the 12th house is 0.094596
@@ -306,7 +319,8 @@


<pre class="programlisting">cout << "Probability that Pat finishes on the " << all_houses
- << " house is " << pdf(nb, all_houses - sales_quota) << endl;</pre>
+ << " house is " << pdf(nb, all_houses - sales_quota) << endl;
+</pre>

 

@@ -326,7 +340,8 @@

<pre class="programlisting">cout << "\nProbability of selling his quota of " << sales_quota
 << " bars\non or before the " << all_houses << "th house is "
- << cdf(nb, all_houses - sales_quota) << endl;</pre>
+ << cdf(nb, all_houses - sales_quota) << endl;
+</pre>

 
<pre class="programlisting">Probability of selling his quota of 5 bars
@@ -367,7 +382,8 @@

<pre class="programlisting">cout << "\nProbability of failing to sell his quota of " << sales_quota
 << " bars\neven after visiting all " << all_houses << " houses is "
- << cdf(complement(nb, all_houses - sales_quota)) << endl;</pre>
+ << cdf(complement(nb, all_houses - sales_quota)) << endl;
+</pre>

 
<pre class="programlisting">Probability of failing to sell his quota of 5 bars
@@ -380,7 +396,8 @@


<pre class="programlisting">double p = cdf(nb, (8 - sales_quota));
-cout << "Probability of meeting sales quota on or before 8th house is "<< p << endl;</pre>
+cout << "Probability of meeting sales quota on or before 8th house is "<< p << endl;
+</pre>

 
<pre class="programlisting">Probability of meeting sales quota on or before 8th house is 0.174
@@ -390,7 +407,8 @@
<pre class="programlisting">cout << "If the confidence of meeting sales quota is " << p
 << ", then the finishing house is " << quantile(nb, p) + sales_quota << endl;

-cout<< " quantile(nb, p) = " << quantile(nb, p) << endl;</pre>
+cout<< " quantile(nb, p) = " << quantile(nb, p) << endl;
+</pre>

 
<pre class="programlisting">If the confidence of meeting sales quota is 0.17367, then the finishing house is 8
@@ -405,7 +423,8 @@

<pre class="programlisting">cout << "If the confidence of meeting sales quota is " << 1.
 << ", then the finishing house is " << quantile(nb, 1) + sales_quota << endl;
-// 1.#INF == infinity.</pre>
+// 1.#INF == infinity.
+</pre>

 
<pre class="programlisting">If the confidence of meeting sales quota is 1, then the finishing house is 1.#INF
@@ -422,7 +441,8 @@
 << ", then the finishing house is " << quantile(nb, 0.5) + sales_quota << endl;

cout << "If the confidence of meeting sales quota is " << 1 - 0.00151 // 30 th
- << ", then the finishing house is " << quantile(nb, 1 - 0.00151) + sales_quota << endl;</pre>
+ << ", then the finishing house is " << quantile(nb, 1 - 0.00151) + sales_quota << endl;
+</pre>

 
<pre class="programlisting">If the confidence of meeting sales quota is 0, then the finishing house is 5
@@ -441,7 +461,8 @@


<pre class="programlisting">cout << "If confidence of meeting quota is zero\n(we assume all houses are successful sales)"
- ", then finishing house is " << sales_quota << endl;</pre>
+ ", then finishing house is " << sales_quota << endl;
+</pre>

 
<pre class="programlisting">If confidence of meeting quota is zero (we assume all houses are successful sales), then finishing house is 5
@@ -460,7 +481,8 @@
 cout << "If confidence of meeting quota is " << ps[i]
 << ", then finishing house is " << quantile(nb, ps[i]) + sales_quota
 << endl;
-}</pre>
+}
+</pre>

 
<pre class="programlisting">If confidence of meeting quota is 0, then finishing house is 5
@@ -514,7 +536,8 @@
{
 cout << left << setw(3) << i << " " << setw(8) << cdf(nb, i - sales_quota) << endl;</spa
n>
}
-cout << endl;</pre>
+cout << endl;
+</pre>

 
<pre class="programlisting">House for 5 th (last) sale. Probability (%)
@@ -555,7 +578,8 @@
catch(const std::exception& e)
{ // Since we have set an overflow policy of ignore_error,
 // an overflow exception should never be thrown.
- std::cout << "\nMessage from thrown exception was:\n " << e.what() << std::endl;</pre>
+ std::cout << "\nMessage from thrown exception was:\n " << e.what() << std::endl;
+</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/normal_example/normal_misc.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/normal_example/normal_misc.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/normal_example/normal_misc.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -58,7 +58,8 @@

 try
 {
- { // Traditional tables and values.</pre>
+ { // Traditional tables and values.
+</pre>

 

@@ -77,7 +78,8 @@
// Construct a standard normal distribution s
 normal s; // (default mean = zero, and standard deviation = unity)
 cout << "Standard normal distribution, mean = "<< s.mean()
- << ", standard deviation = " << s.standard_deviation() << endl;</pre>
+ << ", standard deviation = " << s.standard_deviation() << endl;
+</pre>

 

@@ -93,7 +95,8 @@
 cout << left << setprecision(3) << setw(6) << z << " "
 << setprecision(precision) << setw(12) << pdf(s, z) << endl;
}
-cout.precision(6); // default</pre>
+cout.precision(6); // default
+</pre>

 

@@ -112,7 +115,8 @@
 cout << left << setprecision(3) << setw(6) << z << " "
 << setprecision(precision) << setw(12) << cdf(s, z) << endl;
}
-cout.precision(6); // default</pre>
+cout.precision(6); // default
+</pre>

 

@@ -134,7 +138,8 @@


<pre class="programlisting">double z = 2.;
-cout << "Area for z = " << z << " is " << cdf(s, z) << endl; // to get the area for z.</pre>
+cout << "Area for z = " << z << " is " << cdf(s, z) << endl; // to get the area for z.
+</pre>

 

@@ -146,7 +151,8 @@


<pre class="programlisting"> cout << "95% of area has a z below " << quantile(s, 0.95) << endl;
-// 95% of area has a z below 1.64485</pre>
+// 95% of area has a z below 1.64485
+</pre>

 

@@ -157,7 +163,8 @@

<pre class="programlisting"> cout << "95% of area has a z between " << quantile(s, 0.975)
 << " and " << -quantile(s, 0.975) << endl;
-// 95% of area has a z between 1.95996 and -1.95996</pre>
+// 95% of area has a z between 1.95996 and -1.95996
+</pre>

 

@@ -172,7 +179,8 @@


<pre class="programlisting">double alpha1 = cdf(s, -1) * 2; // 0.3173105078629142
-cout << setprecision(17) << "Significance level for z == 1 is " << alpha1 << endl;</pre>
+cout << setprecision(17) << "Significance level for z == 1 is " << alpha1 << endl;
+</pre>

 

@@ -181,7 +189,8 @@


<pre class="programlisting">double alpha[] = {0.3173105078629142, // z for 1 standard deviation.
- 0.20, 0.1, 0.05, 0.01, 0.001, 0.0001, 0.00001 };</pre>
+ 0.20, 0.1, 0.05, 0.01, 0.001, 0.0001, 0.00001 };
+</pre>

 

@@ -198,7 +207,8 @@
 cout << setw(15) << alpha[i] << setw(15) << alpha[i] /2 << setw(10) << quantile(complement(s, alpha[i]/2)) << endl;
 // Use quantile(complement(s, alpha[i]/2)) to avoid potential loss of accuracy from quantile(s, 1 - alpha[i]/2) 
}
-cout << endl;</pre>
+cout << endl;
+</pre>

 

@@ -241,7 +251,8 @@
cout << "Fraction 2 standard deviations within either side of mean is "
 << 1 - cdf(complement(s, 2 * s.standard_deviation())) * 2 << endl;
cout << "Fraction 3 standard deviations within either side of mean is "
- << 1 - cdf(complement(s, 3 * s.standard_deviation())) * 2 << endl;</pre>
+ << 1 - cdf(complement(s, 3 * s.standard_deviation())) * 2 << endl;
+</pre>

 

@@ -289,7 +300,8 @@


<pre class="programlisting">// K. Krishnamoorthy, Handbook of Statistical Distributions with Applications,
- // ISBN 1 58488 635 8, page 125, example 10.3.5</pre>
+ // ISBN 1 58488 635 8, page 125, example 10.3.5
+</pre>

 

@@ -303,7 +315,8 @@
<pre class="programlisting">double mean_life = 1100.;
double life_standard_deviation = 100.;
normal bulbs(mean_life, life_standard_deviation);
-double expected_life = 1000.;</pre>
+double expected_life = 1000.;
+</pre>

 

@@ -321,7 +334,8 @@
cout << "Fraction of bulbs that will last between "
 << min_life << " and " << max_life << " is "
 << cdf(bulbs, max_life) // P(X <= 1200)
- - cdf(bulbs, min_life) << endl; // P(X <= 900)</pre>
+ - cdf(bulbs, min_life) << endl; // P(X <= 900)
+</pre>

 
<div class="note"><table border="0" summary="Note">
@@ -354,7 +368,8 @@
double stock = 160.; // per week.
cout << "Percentage of weeks overstocked "
 << cdf(sacks, stock) * 100. << endl; // P(X <=160)
-// Percentage of weeks overstocked 97.7</pre>
+// Percentage of weeks overstocked 97.7
+</pre>

 

@@ -364,7 +379,8 @@


<pre class="programlisting">double stock_95 = quantile(sacks, 0.95);
-cout << "Store should stock " << int(stock_95) << " sacks to meet 95% of demands." << endl;</pre>
+cout << "Store should stock " << int(stock_95) << " sacks to meet 95% of demands." << endl;
+</pre>

 

@@ -421,7 +437,8 @@
normal nominal_packs(nominal_mean, standard_deviation);
cout << "Setting the packer to " << nominal_mean << " will mean that "
 << "fraction of packs >= " << under_weight
- << " is " << cdf(complement(nominal_packs, under_weight)) << endl;</pre>
+ << " is " << cdf(complement(nominal_packs, under_weight)) << endl;
+</pre>

 

@@ -446,7 +463,8 @@

<pre class="programlisting">double p = 0.05; // wanted p th quantile.
cout << "Quantile of " << p << " = " << quantile(packs, p)
- << ", mean = " << packs.mean() << ", sd = " << packs.standard_deviation() << endl; //</pre>
+ << ", mean = " << packs.mean() << ", sd = " << packs.standard_deviation() << endl; //
+</pre>

 

@@ -470,7 +488,8 @@
cout <<"Fraction of packs >= " << under_weight << " with a mean of " << mean
 << " and standard deviation of " << pack05.standard_deviation()
 << " is " << cdf(complement(pack05, under_weight)) << endl;
-//</pre>
+//
+</pre>

 

@@ -490,7 +509,8 @@

cout <<"Fraction of packs >= " << under_weight << " with a mean of " << mean
 << " and standard deviation of " << pack06.standard_deviation()
- << " is " << cdf(complement(pack06, under_weight)) << endl;</pre>
+ << " is " << cdf(complement(pack06, under_weight)) << endl;
+</pre>

 

@@ -531,7 +551,8 @@
 << " and standard deviation of " << pack95.standard_deviation()
 << " is " << cdf(complement(pack95, under_weight)) << endl;

-// Fraction of packs >= 2.9 with a mean of 3 and standard deviation of 0.0608 is 0.95</pre>
+// Fraction of packs >= 2.9 with a mean of 3 and standard deviation of 0.0608 is 0.95
+</pre>

 

@@ -570,7 +591,6 @@
 << cdf(complement(bolts, top)) << endl;

cout << "95% of bolts are shorter than " << quantile(bolts, 0.95) << endl;
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/main_overview/history1.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/main_overview/history1.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/main_overview/history1.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -29,6 +29,14 @@
</h3></div></div></div>
<h5>
<a name="math_toolkit.main_overview.history1.h0"></a>
+ <a name="math_toolkit.main_overview.history1.boost_1_50"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_50">Boost-1.50</a>
+ </h5>
+<div class="itemizedlist"><ul class="itemizedlist" type="disc"><li class="listitem">
+ Fixed issue in <code class="computeroutput">bessel_jy</code>
+ which causes Y8.5(4π) to yield a NaN.
+ </li></ul></div>
+<h5>
+<a name="math_toolkit.main_overview.history1.h1"></a>
 <a name="math_toolkit.main_overview.history1.boost_1_49"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_49">Boost-1.49</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -72,7 +80,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h1"></a>
+<a name="math_toolkit.main_overview.history1.h2"></a>
 <a name="math_toolkit.main_overview.history1.boost_1_48"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_48">Boost-1.48</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -123,7 +131,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h2"></a>
+<a name="math_toolkit.main_overview.history1.h3"></a>
 <a name="math_toolkit.main_overview.history1.boost_1_47"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_47">Boost-1.47</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -140,7 +148,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h3"></a>
+<a name="math_toolkit.main_overview.history1.h4"></a>
 <a name="math_toolkit.main_overview.history1.boost_1_46_1"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_46_1">Boost-1.46.1</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc"><li class="listitem">
@@ -148,7 +156,7 @@
 <a href="https://svn.boost.org/trac/boost/ticket/5095" target="_top">#5113</a>.
 </li></ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h4"></a>
+<a name="math_toolkit.main_overview.history1.h5"></a>
 <a name="math_toolkit.main_overview.history1.boost_1_46_0"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_46_0">Boost-1.46.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -163,7 +171,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h5"></a>
+<a name="math_toolkit.main_overview.history1.h6"></a>
 <a name="math_toolkit.main_overview.history1.boost_1_45_0"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_45_0">Boost-1.45.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -180,7 +188,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h6"></a>
+<a name="math_toolkit.main_overview.history1.h7"></a>
 <a name="math_toolkit.main_overview.history1.boost_1_44_0"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_44_0">Boost-1.44.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -194,7 +202,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h7"></a>
+<a name="math_toolkit.main_overview.history1.h8"></a>
 <a name="math_toolkit.main_overview.history1.boost_1_41_0"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_41_0">Boost-1.41.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc"><li class="listitem">
@@ -202,7 +210,7 @@
 and its inverse.
 </li></ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h8"></a>
+<a name="math_toolkit.main_overview.history1.h9"></a>
 <a name="math_toolkit.main_overview.history1.boost_1_40_0"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_40_0">Boost-1.40.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -239,7 +247,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h9"></a>
+<a name="math_toolkit.main_overview.history1.h10"></a>
 <a name="math_toolkit.main_overview.history1.boost_1_38_0"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_38_0">Boost-1.38.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -251,14 +259,14 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h10"></a>
+<a name="math_toolkit.main_overview.history1.h11"></a>
 <a name="math_toolkit.main_overview.history1.boost_1_37_0"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_37_0">Boost-1.37.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc"><li class="listitem">
 Improved accuracy and testing of the inverse hypergeometric functions.
 </li></ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h11"></a>
+<a name="math_toolkit.main_overview.history1.h12"></a>
 <a name="math_toolkit.main_overview.history1.boost_1_36_0"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_36_0">Boost-1.36.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -291,7 +299,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h12"></a>
+<a name="math_toolkit.main_overview.history1.h13"></a>
 <a name="math_toolkit.main_overview.history1.boost_1_35_0__post_review_first_official_release"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.boost_1_35_0__post_review_first_official_release">Boost-1.35.0:
 Post Review First Official Release</a>
 </h5>
@@ -323,7 +331,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h13"></a>
+<a name="math_toolkit.main_overview.history1.h14"></a>
 <a name="math_toolkit.main_overview.history1.milestone_4__second_review_candidate__1st_march_2007_"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.milestone_4__second_review_candidate__1st_march_2007_">Milestone
 4: Second Review Candidate (1st March 2007)</a>
 </h5>
@@ -337,7 +345,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h14"></a>
+<a name="math_toolkit.main_overview.history1.h15"></a>
 <a name="math_toolkit.main_overview.history1.milestone_3__first_review_candidate__31st_dec_2006_"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.milestone_3__first_review_candidate__31st_dec_2006_">Milestone
 3: First Review Candidate (31st Dec 2006)</a>
 </h5>
@@ -365,7 +373,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h15"></a>
+<a name="math_toolkit.main_overview.history1.h16"></a>
 <a name="math_toolkit.main_overview.history1.milestone_2__released_september_10th_2006"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.milestone_2__released_september_10th_2006">Milestone
 2: Released September 10th 2006</a>
 </h5>
@@ -401,7 +409,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.main_overview.history1.h16"></a>
+<a name="math_toolkit.main_overview.history1.h17"></a>
 <a name="math_toolkit.main_overview.history1.milestone_1__released_march_31st_2006"></a><a class="link" href="history1.html#math_toolkit.main_overview.history1.milestone_1__released_march_31st_2006">Milestone
 1: Released March 31st 2006</a>
 </h5>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/discrete_quant_ref.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/discrete_quant_ref.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/discrete_quant_ref.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -78,7 +78,6 @@
double x = quantile(dist_type(20, 0.3), 0.05);
// Upper 95% quantile:
double y = quantile(complement(dist_type(20, 0.3), 0.05));
-
</pre>

 
@@ -119,7 +118,6 @@
double x = quantile(negative_binomial(20, 0.3), 0.05); // rounded up 27 from 27.3898
// Upper quantile rounded up:
double y = quantile(complement(negative_binomial(20, 0.3), 0.05)); // rounded down to 69 from 68.1584
-
</pre>

 
@@ -177,7 +175,6 @@
double x = quantile(dist_type(20, 0.3), 0.05); // 28 rounded up from 27.3898
// Upper quantile rounded down:
double y = quantile(complement(dist_type(20, 0.3), 0.05)); // 68 rounded down from 68.1584
-
</pre>

 
@@ -244,7 +241,6 @@
double x = quantile(dist_type(20, 0.3), 0.05); // 27
// Upper quantile rounded (down) to nearest:
double y = quantile(complement(dist_type(20, 0.3), 0.05)); // 68
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/error_handling_policies.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/error_handling_policies.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/error_handling_policies.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -760,7 +760,6 @@
// Construct a my_norm distribution, using default mean and standard deviation,
// and get a 0.05 or 5% quantile:
double q = quantile(my_norm(), 0.05); // = -1.64485
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/internal_promotion.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/internal_promotion.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/internal_promotion.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -111,7 +111,6 @@
// Alternatively we could use helper function make_policy,
// and concisely define everything at the call site:
double t2 = tgamma(some_value, make_policy(promote_double<false>()));
-
</pre>


@@ -136,7 +135,6 @@

// Get a quantile:
float q = quantile(my_norm(), 0.05f);
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/namespace_pol.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/namespace_pol.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/namespace_pol.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -68,7 +68,6 @@
double t = myspace::tgamma(30.0); // Will *not* throw on overflow,
// despite the large value of factorial 30 = 265252859812191058636308480000000
// unlike default policy boost::math::tgamma;
-
</pre>

 
@@ -133,7 +132,6 @@
 cout << e.what() << endl;
 }
}
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/precision_pol.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/precision_pol.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_ref/precision_pol.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -82,7 +82,6 @@
using namespace boost::math::policies;

double t = tgamma(12, policy<digits10<5> >()); // Concise make_policy.
-
</pre>

 
@@ -101,7 +100,6 @@
double q = quantile(
 normal_distribution<double, policy<digits2<bits> > >(),
 0.05); // 5% quantile.
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/ad_hoc_dist_policies.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/ad_hoc_dist_policies.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/ad_hoc_dist_policies.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -83,7 +83,6 @@
 cout << "quantile(mybinom(200, 0.25), 0.05) is: " <<
 quantile(mybinom(200, 0.25), 0.05) << endl;
}
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/ad_hoc_sf_policies.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/ad_hoc_sf_policies.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/ad_hoc_sf_policies.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -101,7 +101,6 @@
 << boost::math::tgamma(-10, c_policy()) << endl; // Negative parameter.
 cout << "errno = " << errno << endl; // error 33 Numerical argument out of domain.
} // int main()
-
</pre>

 
@@ -160,7 +159,6 @@
 // Check errno was set:
 std::cout << "errno = " << errno << std::endl;
}
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/changing_policy_defaults.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/changing_policy_defaults.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/changing_policy_defaults.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -175,7 +175,8 @@
 using std::exception;

#include <cstddef>
- // using ::errno</pre>
+ // using ::errno
+</pre>

 

@@ -217,7 +218,8 @@
 "\n""Message from thrown exception was:\n " << e.what() << std::endl;
 }
 return 0;
-} // int main()</pre>
+} // int main()
+</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/namespace_policies.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/namespace_policies.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/namespace_policies.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -53,7 +53,8 @@


<pre class="programlisting">#include <boost/math/special_functions.hpp>
-//using boost::math::tgamma; // Not needed because using C::tgamma.</pre>
+//using boost::math::tgamma; // Not needed because using C::tgamma.
+</pre>

 

@@ -82,7 +83,8 @@
 pole_error<errno_on_error>,
 overflow_error<errno_on_error>,
 evaluation_error<errno_on_error>
- > c_policy;</pre>
+ > c_policy;
+</pre>

 

@@ -93,7 +95,8 @@

<pre class="programlisting">BOOST_MATH_DECLARE_SPECIAL_FUNCTIONS(c_policy)

-} // close namespace C</pre>
+} // close namespace C
+</pre>

 

@@ -127,7 +130,8 @@
 cout << "Result of tgamma(-10) is: "
 << C::tgamma(-10) << endl;
 cout << "errno = " << errno << endl; // errno = 33, overwriting previous value of 34.
-}</pre>
+}
+</pre>

 

@@ -167,7 +171,8 @@
 evaluation_error<errno_on_error>
> c_policy;

-BOOST_MATH_DECLARE_SPECIAL_FUNCTIONS(c_policy)</pre>
+BOOST_MATH_DECLARE_SPECIAL_FUNCTIONS(c_policy)
+</pre>

 

@@ -191,7 +196,6 @@
 cout << "errno = " << errno << endl;
 // Default tgamma policy would throw an exception, and abort.
}
-
</pre>

 
@@ -230,7 +234,8 @@
 


-<pre class="programlisting">#include <boost/math/distributions.hpp></pre>
+<pre class="programlisting">#include <boost/math/distributions.hpp>
+</pre>

 

@@ -252,7 +257,8 @@
 promote_double<false>,
 // Return the closest integer result for discrete quantiles:
 discrete_quantile<integer_round_nearest>
- > my_policy;</pre>
+ > my_policy;
+</pre>

 

@@ -265,7 +271,8 @@

<pre class="programlisting">BOOST_MATH_DECLARE_DISTRIBUTIONS(double, my_policy)

-} // close namespace my_namespace</pre>
+} // close namespace my_namespace
+</pre>

 

@@ -307,7 +314,8 @@
 << quantile(binom, 0.05) << endl; // To check we get integer results.
 cout << "Result of quantile(complement(binom, 0.05)) is: "
 << quantile(complement(binom, 0.05)) << endl;
-}</pre>
+}
+</pre>

 

@@ -387,7 +395,6 @@
 cout << "Result of quantile(complement(binom, 0.05)) is: "
 << quantile(complement(binom, 0.05)) << endl;
}
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/understand_dis_quant.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/understand_dis_quant.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/understand_dis_quant.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -256,7 +256,8 @@
using std::setw;
using std::setprecision;

-#include <boost/math/distributions/binomial.hpp></pre>
+#include <boost/math/distributions/binomial.hpp>
+</pre>

 

@@ -313,7 +314,8 @@
typedef binomial_distribution<
 double,
 policy<discrete_quantile<real> > >
- binom_real_quantile;</pre>
+ binom_real_quantile;
+</pre>

 

@@ -387,7 +389,8 @@
 << setw(18) << right
 << quantile(binom_real_quantile(50, 0.5), 0.95)
 << endl;
-} // int main()</pre>
+} // int main()
+</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/user_def_err_pol.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/user_def_err_pol.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/policy/pol_tutorial/user_def_err_pol.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -84,7 +84,8 @@
 return std::numeric_limits<T>::quiet_NaN();
 }
 } // namespace policies
-}} // namespace boost{ namespace math</pre>
+}} // namespace boost{ namespace math
+</pre>

 

@@ -104,7 +105,8 @@

BOOST_MATH_DECLARE_SPECIAL_FUNCTIONS(user_error_policy)

-} // close unnamed namespace</pre>
+} // close unnamed namespace
+</pre>

 

@@ -136,7 +138,8 @@
 << erf_inv(-10) << endl;
 cout << "Result of tgamma(-10) is: "
 << tgamma(-10) << endl;
-}</pre>
+}
+</pre>

 

@@ -196,7 +199,8 @@
 


-<pre class="programlisting">#include <boost/math/special_functions.hpp></pre>
+<pre class="programlisting">#include <boost/math/special_functions.hpp>
+</pre>

 

@@ -211,7 +215,8 @@

template <class T>
T user_domain_error(const char* function, const char* message, const T& val)
-{</pre>
+{
+</pre>

 

@@ -223,7 +228,8 @@
<pre class="programlisting">if(function == 0)
 function = "Unknown function with arguments of type %1%";
if(message == 0)
- message = "Cause unknown with bad argument %1%";</pre>
+ message = "Cause unknown with bad argument %1%";
+</pre>

 

@@ -233,7 +239,8 @@


<pre class="programlisting">std::string msg("Error in function ");
-msg += (boost::format(function) % typeid(T).name()).str();</pre>
+msg += (boost::format(function) % typeid(T).name()).str();
+</pre>

 

@@ -245,7 +252,8 @@
<pre class="programlisting">msg += ": \n";
int prec = 2 + (std::numeric_limits<T>::digits * 30103UL) / 100000UL;
// int prec = std::numeric_limits<T>::max_digits10; // For C++0X Standard Library
-msg += (boost::format(message) % boost::io::group(std::setprecision(prec), val)).str();</pre>
+msg += (boost::format(message) % boost::io::group(std::setprecision(prec), val)).str();
+</pre>

 

@@ -254,7 +262,8 @@
 


-<pre class="programlisting">std::cerr << msg << std::endl;</pre>
+<pre class="programlisting">std::cerr << msg << std::endl;
+</pre>

 

@@ -264,7 +273,8 @@


<pre class="programlisting"> return std::numeric_limits<T>::quiet_NaN();
-}</pre>
+}
+</pre>

 

@@ -277,7 +287,8 @@
T user_pole_error(const char* function, const char* message, const T& val)
{
 return user_domain_error(function, message, val);
-}</pre>
+}
+</pre>

 

@@ -304,7 +315,8 @@

 // Value passed to the function is an infinity, just return it:
 return val;
-}</pre>
+}
+</pre>

 

@@ -330,7 +342,8 @@

 // Value passed to the function is zero, just return it:
 return val;
-}</pre>
+}
+</pre>

 

@@ -356,7 +369,8 @@

 // Value passed to the function is denormalised, just return it:
 return val;
-}</pre>
+}
+</pre>

 

@@ -392,7 +406,8 @@
}

} // policies
-}} // boost::math</pre>
+}} // boost::math
+</pre>

 

@@ -417,7 +432,8 @@

BOOST_MATH_DECLARE_SPECIAL_FUNCTIONS(user_error_policy)

-} // unnamed namespace</pre>
+} // unnamed namespace
+</pre>

 

@@ -460,7 +476,8 @@
 // Unfortunately we can't predicably raise a denormalised
 // result, nor can we raise an evaluation error in this example
 // since these should never really occur!
-} // int main()</pre>
+} // int main()
+</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/status/history1.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/status/history1.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/status/history1.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -28,6 +28,14 @@
</h3></div></div></div>
<h5>
<a name="math_toolkit.status.history1.h0"></a>
+ <a name="math_toolkit.status.history1.boost_1_50"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_50">Boost-1.50</a>
+ </h5>
+<div class="itemizedlist"><ul class="itemizedlist" type="disc"><li class="listitem">
+ Fixed issue in <code class="computeroutput">bessel_jy</code>
+ which causes Y8.5(4π) to yield a NaN.
+ </li></ul></div>
+<h5>
+<a name="math_toolkit.status.history1.h1"></a>
 <a name="math_toolkit.status.history1.boost_1_49"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_49">Boost-1.49</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -71,7 +79,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.status.history1.h1"></a>
+<a name="math_toolkit.status.history1.h2"></a>
 <a name="math_toolkit.status.history1.boost_1_48"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_48">Boost-1.48</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -122,7 +130,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.status.history1.h2"></a>
+<a name="math_toolkit.status.history1.h3"></a>
 <a name="math_toolkit.status.history1.boost_1_47"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_47">Boost-1.47</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -139,7 +147,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.status.history1.h3"></a>
+<a name="math_toolkit.status.history1.h4"></a>
 <a name="math_toolkit.status.history1.boost_1_46_1"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_46_1">Boost-1.46.1</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc"><li class="listitem">
@@ -147,7 +155,7 @@
 <a href="https://svn.boost.org/trac/boost/ticket/5095" target="_top">#5113</a>.
 </li></ul></div>
<h5>
-<a name="math_toolkit.status.history1.h4"></a>
+<a name="math_toolkit.status.history1.h5"></a>
 <a name="math_toolkit.status.history1.boost_1_46_0"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_46_0">Boost-1.46.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -162,7 +170,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.status.history1.h5"></a>
+<a name="math_toolkit.status.history1.h6"></a>
 <a name="math_toolkit.status.history1.boost_1_45_0"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_45_0">Boost-1.45.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -179,7 +187,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.status.history1.h6"></a>
+<a name="math_toolkit.status.history1.h7"></a>
 <a name="math_toolkit.status.history1.boost_1_44_0"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_44_0">Boost-1.44.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -193,7 +201,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.status.history1.h7"></a>
+<a name="math_toolkit.status.history1.h8"></a>
 <a name="math_toolkit.status.history1.boost_1_41_0"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_41_0">Boost-1.41.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc"><li class="listitem">
@@ -201,7 +209,7 @@
 and its inverse.
 </li></ul></div>
<h5>
-<a name="math_toolkit.status.history1.h8"></a>
+<a name="math_toolkit.status.history1.h9"></a>
 <a name="math_toolkit.status.history1.boost_1_40_0"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_40_0">Boost-1.40.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -238,7 +246,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.status.history1.h9"></a>
+<a name="math_toolkit.status.history1.h10"></a>
 <a name="math_toolkit.status.history1.boost_1_38_0"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_38_0">Boost-1.38.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -250,14 +258,14 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.status.history1.h10"></a>
+<a name="math_toolkit.status.history1.h11"></a>
 <a name="math_toolkit.status.history1.boost_1_37_0"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_37_0">Boost-1.37.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc"><li class="listitem">
 Improved accuracy and testing of the inverse hypergeometric functions.
 </li></ul></div>
<h5>
-<a name="math_toolkit.status.history1.h11"></a>
+<a name="math_toolkit.status.history1.h12"></a>
 <a name="math_toolkit.status.history1.boost_1_36_0"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_36_0">Boost-1.36.0</a>
 </h5>
<div class="itemizedlist"><ul class="itemizedlist" type="disc">
@@ -290,7 +298,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.status.history1.h12"></a>
+<a name="math_toolkit.status.history1.h13"></a>
 <a name="math_toolkit.status.history1.boost_1_35_0__post_review_first_official_release"></a><a class="link" href="history1.html#math_toolkit.status.history1.boost_1_35_0__post_review_first_official_release">Boost-1.35.0:
 Post Review First Official Release</a>
 </h5>
@@ -322,7 +330,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.status.history1.h13"></a>
+<a name="math_toolkit.status.history1.h14"></a>
 <a name="math_toolkit.status.history1.milestone_4__second_review_candidate__1st_march_2007_"></a><a class="link" href="history1.html#math_toolkit.status.history1.milestone_4__second_review_candidate__1st_march_2007_">Milestone
 4: Second Review Candidate (1st March 2007)</a>
 </h5>
@@ -336,7 +344,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.status.history1.h14"></a>
+<a name="math_toolkit.status.history1.h15"></a>
 <a name="math_toolkit.status.history1.milestone_3__first_review_candidate__31st_dec_2006_"></a><a class="link" href="history1.html#math_toolkit.status.history1.milestone_3__first_review_candidate__31st_dec_2006_">Milestone
 3: First Review Candidate (31st Dec 2006)</a>
 </h5>
@@ -364,7 +372,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.status.history1.h15"></a>
+<a name="math_toolkit.status.history1.h16"></a>
 <a name="math_toolkit.status.history1.milestone_2__released_september_10th_2006"></a><a class="link" href="history1.html#math_toolkit.status.history1.milestone_2__released_september_10th_2006">Milestone
 2: Released September 10th 2006</a>
 </h5>
@@ -400,7 +408,7 @@
 </li>
</ul></div>
<h5>
-<a name="math_toolkit.status.history1.h16"></a>
+<a name="math_toolkit.status.history1.h17"></a>
 <a name="math_toolkit.status.history1.milestone_1__released_march_31st_2006"></a><a class="link" href="history1.html#math_toolkit.status.history1.milestone_1__released_march_31st_2006">Milestone
 1: Released March 31st 2006</a>
 </h5>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/utils/fp_facets/examples.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/utils/fp_facets/examples.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/utils/fp_facets/examples.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -57,7 +57,6 @@
// with the new locale containing the nonfinite_numput facet.
// So the cout output will be "1.#INF on MS platforms
// and may be "inf" or other string representation on other platforms.
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/utils/fp_facets/intro.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/utils/fp_facets/intro.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/utils/fp_facets/intro.html 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -124,7 +124,6 @@
// with the new locale containing the nonfinite_numput facet.
// So the cout output will be "1.#INF on MS platforms
// and may be "inf" or other string representation on other platforms.
-
</pre>

Modified: trunk/libs/math/doc/sf_and_dist/roadmap.qbk
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/roadmap.qbk (original)
+++ trunk/libs/math/doc/sf_and_dist/roadmap.qbk 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -1,5 +1,9 @@
[template history[]

+[h4 Boost-1.50]
+
+* Fixed issue in `bessel_jy` which causes Y[sub 8.5](4[pi]) to yield a NaN.
+
[h4 Boost-1.49]

* Deprecated wrongly named `twothirds` math constant in favour of `two_thirds` (with underscore separator).

Modified: trunk/libs/math/test/pch.hpp
==============================================================================
--- trunk/libs/math/test/pch.hpp (original)
+++ trunk/libs/math/test/pch.hpp 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -5,14 +5,11 @@

#ifdef BOOST_BUILD_PCH_ENABLED

-#define ONE <boost/math/special_functions.hpp>
-#define TWO <boost/math/distributions.hpp>
-
#define BOOST_MATH_OVERFLOW_ERROR_POLICY ignore_error
#include <boost/test/test_exec_monitor.hpp>
#include <boost/test/floating_point_comparison.hpp>
-#include ONE
-#include TWO
+#include <boost/math/special_functions.hpp>
+#include <boost/math/distributions.hpp>
#include <boost/math/tools/test.hpp>
#include <boost/math/tools/stats.hpp>
#include <boost/type_traits.hpp>

Modified: trunk/libs/math/test/test_bessel_j.cpp
==============================================================================
--- trunk/libs/math/test/test_bessel_j.cpp (original)
+++ trunk/libs/math/test/test_bessel_j.cpp 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -271,6 +271,16 @@
       ".*", // test data group
       ".*", 40, 20); // test function
    //
+ // One set of float tests has inexact input values, so there is a slight error:
+ //
+ add_expected_result(
+ ".*", // compiler
+ ".*", // stdlib
+ ".*", // platform
+ "float", // test type(s)
+ "Bessel J: Mathworld Data", // test data group
+ ".*", 20, 20); // test function
+ //
    // Finish off by printing out the compiler/stdlib/platform names,
    // we do this to make it easier to mark up expected error rates.
    //

Modified: trunk/libs/math/test/test_bessel_j.hpp
==============================================================================
--- trunk/libs/math/test/test_bessel_j.hpp (original)
+++ trunk/libs/math/test/test_bessel_j.hpp 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -8,6 +8,7 @@
#include <boost/test/test_exec_monitor.hpp>
#include <boost/test/floating_point_comparison.hpp>
#include <boost/math/special_functions/math_fwd.hpp>
+#include <boost/math/constants/constants.hpp>
#include <boost/type_traits/is_floating_point.hpp>
#include <boost/array.hpp>
#include "functor.hpp"
@@ -215,7 +216,7 @@
 do_test_cyl_bessel_j_int<T>(j1_tricky, name, "Bessel J1: Mathworld Data (tricky cases) (Integer Version)");
 do_test_cyl_bessel_j_int<T>(jn_data, name, "Bessel JN: Mathworld Data (Integer Version)");

- static const boost::array<boost::array<T, 3>, 18> jv_data = {{
+ static const boost::array<boost::array<T, 3>, 20> jv_data = {{
 //SC_(-2.4), {{ SC_(0), std::numeric_limits<T>::infinity() }},
 {{ T(2457)/1024, T(1)/1024, SC_(3.80739920118603335646474073457326714709615200130620574875292e-9) }},
 {{ SC_(5.5), T(3217)/1024, SC_(0.0281933076257506091621579544064767140470089107926550720453038) }},
@@ -236,6 +237,8 @@
 {{ T(-10486074) / (1024*1024), SC_(1e+02), SC_(-0.0547613660316806551338637153942604550779513947674222863858713) }},
 // Bug report https://svn.boost.org/trac/boost/ticket/4812:
 {{ SC_(1.5), T(8034)/1024, SC_(0.0339477646369710610146236955872928005087352629422508823945264) }},
+ {{ SC_(8.5), boost::math::constants::pi<T>() * 4, SC_(0.0436807946352780974532519564114026730332781693877984686758680) }},
+ {{ SC_(-8.5), boost::math::constants::pi<T>() * 4, SC_(-0.257086543428224355151772807588810984369026142375675714560864) }},
 }};
 do_test_cyl_bessel_j<T>(jv_data, name, "Bessel J: Mathworld Data");
 static const boost::array<boost::array<T, 3>, 4> jv_large_data = {{
@@ -246,7 +249,7 @@
 {{ SC_(-2.5), SC_(4), SC_(-0.0145679476685218007666785535204236327832335803441449596297004) }},
 }};
 if(jv_large_data[0][1] != 0)
- do_test_cyl_bessel_j<T>(jv_data, name, "Bessel J: Mathworld Data (large values)");
+ do_test_cyl_bessel_j<T>(jv_large_data, name, "Bessel J: Mathworld Data (large values)");

#include "bessel_j_int_data.ipp"
 do_test_cyl_bessel_j<T>(bessel_j_int_data, name, "Bessel JN: Random Data");

Modified: trunk/libs/math/test/test_bessel_y.cpp
==============================================================================
--- trunk/libs/math/test/test_bessel_y.cpp (original)
+++ trunk/libs/math/test/test_bessel_y.cpp 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -229,6 +229,16 @@
       ".*", // test data group
       ".*", 80, 40); // test function
    //
+ // One set of float tests has inexact input values, so there is a slight error:
+ //
+ add_expected_result(
+ ".*", // compiler
+ ".*", // stdlib
+ ".*", // platform
+ "float", // test type(s)
+ "Yv: Mathworld Data", // test data group
+ ".*", 20, 20); // test function
+ //
    // Finish off by printing out the compiler/stdlib/platform names,
    // we do this to make it easier to mark up expected error rates.
    //

Modified: trunk/libs/math/test/test_bessel_y.hpp
==============================================================================
--- trunk/libs/math/test/test_bessel_y.hpp (original)
+++ trunk/libs/math/test/test_bessel_y.hpp 2012-03-24 07:59:21 EDT (Sat, 24 Mar 2012)
@@ -8,6 +8,7 @@
#include <boost/test/test_exec_monitor.hpp>
#include <boost/test/floating_point_comparison.hpp>
#include <boost/math/special_functions/math_fwd.hpp>
+#include <boost/math/constants/constants.hpp>
#include <boost/type_traits/is_floating_point.hpp>
#include <boost/array.hpp>
#include "functor.hpp"
@@ -176,7 +177,7 @@
 {{ SC_(-1e+03), SC_(7e+02), SC_(-1.88753109980945889960843803284345261796244752396992106755091e77) }},
 {{ SC_(-25), SC_(8), SC_(3.45113613777297661997458045843868931827873456761831907587263e8) }}
 }};
- static const boost::array<boost::array<T, 3>, 9> yv_data = {{
+ static const boost::array<boost::array<T, 3>, 11> yv_data = {{
 //SC_(2.25), {{ SC_(1) / 1024, SC_(-1.01759203636941035147948317764932151601257765988969544340275e7) }},
 {{ SC_(0.5), T(1) / (1024*1024), SC_(-817.033790261762580469303126467917092806755460418223776544122) }},
 {{ SC_(5.5), SC_(3.125), SC_(-2.61489440328417468776474188539366752698192046890955453259866) }},
@@ -187,6 +188,8 @@
 {{ SC_(141.75), SC_(1e+02), SC_(-5.38829231428696507293191118661269920130838607482708483122068e9) }},
 {{ SC_(141.75), SC_(2e+04), SC_(-0.00376577888677186194728129112270988602876597726657372330194186) }},
 {{ SC_(-141.75), SC_(1e+02), SC_(-3.81009803444766877495905954105669819951653361036342457919021e9) }},
+ {{ SC_(8.5), boost::math::constants::pi<T>() * 4, SC_(0.257086543428224355151772807588810984369026142375675714560864) }},
+ {{ SC_(-8.5), boost::math::constants::pi<T>() * 4, SC_(0.0436807946352780974532519564114026730332781693877984686758680) }},
 }};
 static const boost::array<boost::array<T, 3>, 7> yv_large_data = {{
 // Bug report https://svn.boost.org/trac/boost/ticket/5560:

Next message: jewillco_at_[hidden]: "[Boost-commit] svn:boost r77514 - in trunk: boost/property_map libs/property_map/doc"
Previous message: joel_at_[hidden]: "[Boost-commit] svn:boost r77512 - trunk/boost/spirit/include"

Date view	Thread view	Subject view	Author view

Boost-Commit list run by bdawes at acm.org, david.abrahams at rcn.com, gregod at cs.rpi.edu, cpdaniel at pacbell.net, john at johnmaddock.co.uk