Boost logo

Boost-Commit :

Subject: [Boost-commit] svn:boost r82217 - in trunk: boost/math/special_functions 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 libs/math/doc/sf_and_dist/html/math_toolkit/main_overview libs/math/doc/sf_and_dist/html/math_toolkit/special libs/math/doc/sf_and_dist/html/math_toolkit/special/expint libs/math/doc/sf_and_dist/html/math_toolkit/special/powers libs/math/doc/sf_and_dist/html/math_toolkit/status libs/math/test
From: john_at_[hidden]
Date: 2012-12-27 03:44:12


Author: johnmaddock
Date: 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
New Revision: 82217
URL: http://svn.boost.org/trac/boost/changeset/82217

Log:
More accuracy improvements based on suggestions by Rocco Romeo.
Added:
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/cos_pi.html (contents, props changed)
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/sin_pi.html (contents, props changed)
Text files modified:
   trunk/boost/math/special_functions/bessel.hpp | 4 +-
   trunk/boost/math/special_functions/detail/bessel_j0.hpp | 5 +--
   trunk/boost/math/special_functions/detail/bessel_jy.hpp | 24 ++++++++++++----
   trunk/boost/math/special_functions/detail/bessel_y0.hpp | 8 ++--
   trunk/boost/math/special_functions/detail/bessel_yn.hpp | 2
   trunk/libs/math/doc/sf_and_dist/html/index.html | 7 ++--
   trunk/libs/math/doc/sf_and_dist/html/index/s13.html | 6 ++-
   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 | 6 ++-
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/main_overview/conventions.html | 2
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/main_overview/history1.html | 4 ++
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/main_overview/navigation.html | 2
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special.html | 5 ++-
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/expint/expint_i.html | 2
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers.html | 13 +++++----
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/cbrt.html | 2
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/ct_pow.html | 2
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/expm1.html | 2
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/hypot.html | 2
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/log1p.html | 8 ++--
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/powm1.html | 2
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/sqrt1pm1.html | 2
   trunk/libs/math/doc/sf_and_dist/html/math_toolkit/status/history1.html | 4 ++
   trunk/libs/math/doc/sf_and_dist/powers.qbk | 56 +++++++++++++++++++++++++++++++++++++++
   trunk/libs/math/test/test_bessel_j.cpp | 9 -----
   trunk/libs/math/test/test_bessel_y.cpp | 7 -----
   28 files changed, 132 insertions(+), 66 deletions(-)

Modified: trunk/boost/math/special_functions/bessel.hpp
==============================================================================
--- trunk/boost/math/special_functions/bessel.hpp (original)
+++ trunk/boost/math/special_functions/bessel.hpp 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -301,7 +301,7 @@
 
    if(floor(v) == v)
    {
- if((fabs(x) > asymptotic_bessel_y_limit<T>(tag_type())) && (fabs(x) > 5 * abs(v)))
+ if((fabs(x) > asymptotic_bessel_y_limit<T>(tag_type())) && (4 * fabs(x) * sqrt(sqrt(sqrt(14 * tools::epsilon<T>() * (fabs(x) + fabs(v) / 2) / (5120 * fabs(x))))) > fabs(v)))
       {
          T r = asymptotic_bessel_y_large_x_2(static_cast<T>(abs(v)), x);
          if((v < 0) && (itrunc(v, pol) & 1))
@@ -330,7 +330,7 @@
    BOOST_MATH_INSTRUMENT_VARIABLE(v);
    BOOST_MATH_INSTRUMENT_VARIABLE(x);
 
- if((fabs(x) > asymptotic_bessel_y_limit<T>(tag_type())) && (fabs(x) > 5 * abs(v)))
+ if((fabs(x) > asymptotic_bessel_y_limit<T>(tag_type())) && (4 * fabs(x) * sqrt(sqrt(sqrt(14 * tools::epsilon<T>() * (fabs(x) + abs(v) / 2) / (5120 * x)))) > abs(v)))
    {
       T r = asymptotic_bessel_y_large_x_2(static_cast<T>(abs(v)), x);
       if((v < 0) && (v & 1))

Modified: trunk/boost/math/special_functions/detail/bessel_j0.hpp
==============================================================================
--- trunk/boost/math/special_functions/detail/bessel_j0.hpp (original)
+++ trunk/boost/math/special_functions/detail/bessel_j0.hpp 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -169,7 +169,7 @@
         BOOST_ASSERT(sizeof(PS) == sizeof(QS));
         rc = evaluate_rational(PC, QC, y2);
         rs = evaluate_rational(PS, QS, y2);
- factor = sqrt(2 / (x * pi<T>()));
+ factor = constants::one_div_root_pi<T>() / sqrt(x);
         //
         // What follows is really just:
         //
@@ -181,8 +181,7 @@
         //
         T sx = sin(x);
         T cx = cos(x);
- value = factor * (rc * (cx * constants::one_div_root_two<T>() + sx * constants::half_root_two<T>())
- - y * rs * (sx * constants::one_div_root_two<T>() - cx * constants::half_root_two<T>()));
+ value = factor * (rc * (cx + sx) - y * rs * (sx - cx));
     }
 
     return value;

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-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -307,7 +307,7 @@
     // x is positive until reflection
     W = T(2) / (x * pi<T>()); // Wronskian
     T Yv_scale = 1;
- if((x > 8) && (x < 1000) && hankel_PQ(v, x, &p, &q, pol))
+ if((x > 8) && hankel_PQ(v, x, &p, &q, pol))
     {
        //
        // Hankel approximation: note that this method works best when x
@@ -315,10 +315,22 @@
        // of large values, with horrendous resulting accuracy. It is fast though
        // when it works....
        //
- T chi = x - fmod(T(v / 2 + 0.25f), T(2)) * boost::math::constants::pi<T>();
- T sc = sin(chi);
- T cc = cos(chi);
- chi = sqrt(2 / (boost::math::constants::pi<T>() * x));
+ // Normally we calculate sin/cos(chi) where:
+ //
+ // chi = x - fmod(T(v / 2 + 0.25f), T(2)) * boost::math::constants::pi<T>();
+ //
+ // But this introduces large errors, so use sin/cos addition formulae to
+ // improve accuracy:
+ //
+ T mod_v = fmod(T(v / 2 + 0.25f), T(2));
+ T sx = sin(x);
+ T cx = cos(x);
+ T sv = sin_pi(mod_v);
+ T cv = cos_pi(mod_v);
+
+ T sc = sx * cv - sv * cx; // == sin(chi);
+ T cc = cx * cv + sx * sv; // == cos(chi);
+ T chi = boost::math::constants::root_two<T>() / (boost::math::constants::root_pi<T>() * sqrt(x)); //sqrt(2 / (boost::math::constants::pi<T>() * x));
        Yv = chi * (p * sc + q * cc);
        Jv = chi * (p * cc - q * sc);
     }
@@ -459,7 +471,7 @@
               for (k = n; k > 0; k--) // backward recurrence for J
               {
                   T t = 2 * (u + k) / x;
- if(tools::max_value<T>() / t < current)
+ if((t > 1) && (tools::max_value<T>() / t < current))
                   {
                      over = true;
                      break;

Modified: trunk/boost/math/special_functions/detail/bessel_y0.hpp
==============================================================================
--- trunk/boost/math/special_functions/detail/bessel_y0.hpp (original)
+++ trunk/boost/math/special_functions/detail/bessel_y0.hpp 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -199,7 +199,7 @@
         T y2 = y * y;
         rc = evaluate_rational(PC, QC, y2);
         rs = evaluate_rational(PS, QS, y2);
- factor = sqrt(2 / (x * pi<T>()));
+ factor = constants::one_div_root_pi<T>() / sqrt(x);
         //
         // The following code is really just:
         //
@@ -207,12 +207,12 @@
         // value = factor * (rc * sin(z) + y * rs * cos(z));
         //
         // But using the sin/cos addition formulae and constant values for
- // sin/cos of PI/4:
+ // sin/cos of PI/4 which then cancel part of the "factor" term as they're all
+ // 1 / sqrt(2):
         //
         T sx = sin(x);
         T cx = cos(x);
- value = factor * (rc * (sx * constants::one_div_root_two<T>() - cx * constants::half_root_two<T>())
- + y * rs * (cx * constants::one_div_root_two<T>() + sx * constants::half_root_two<T>()));
+ value = factor * (rc * (sx - cx) + y * rs * (cx + sx));
     }
 
     return value;

Modified: trunk/boost/math/special_functions/detail/bessel_yn.hpp
==============================================================================
--- trunk/boost/math/special_functions/detail/bessel_yn.hpp (original)
+++ trunk/boost/math/special_functions/detail/bessel_yn.hpp 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -78,7 +78,7 @@
        do
        {
            T fact = 2 * k / x;
- if((tools::max_value<T>() - fabs(prev)) / fact < fabs(current))
+ if((fact > 1) && ((tools::max_value<T>() - fabs(prev)) / fact < fabs(current)))
            {
               prev /= current;
               factor /= current;

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-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -355,9 +355,10 @@
 <dt><span class="section"><a href="math_toolkit/special/expint/expint_i.html">Exponential Integral
         Ei</a></span></dt>
 </dl></dd>
-<dt><span class="section"><a href="math_toolkit/special/powers.html">Logs, Powers, Roots and
- Exponentials</a></span></dt>
+<dt><span class="section">Basic Functions</span></dt>
 <dd><dl>
+<dt><span class="section">sin_pi</span></dt>
+<dt><span class="section">cos_pi</span></dt>
 <dt><span class="section">log1p</span></dt>
 <dt><span class="section">expm1</span></dt>
 <dt><span class="section">cbrt</span></dt>
@@ -580,7 +581,7 @@
   </p>
 </div>
 <table xmlns:rev="http://www.cs.rpi.edu/~gregod/boost/tools/doc/revision" width="100%"><tr>
-<td align="left"><p><small>Last revised: December 12, 2012 at 13:18:44 GMT</small></p></td>
+<td align="left"><p><small>Last revised: December 24, 2012 at 12:22:13 GMT</small></p></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-12-27 03:44:08 EST (Thu, 27 Dec 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 id1334003">
+<div class="section id1323319">
 <div class="titlepage"><div><div><h2 class="title" style="clear: both">
-<a name="id1334003"></a>Function Index</h2></div></div></div>
+<a name="id1323319"></a>Function Index</h2></div></div></div>
 <p><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">J</a> <a class="link" href="s13.html#idx_id_10">K</a> <a class="link" href="s13.html#idx_id_11">L</a> <a class="link" href="s13.html#idx_id_12">M</a> <a class="link" href="s13.html#idx_id_13">N</a> <a class="link" href="s13.html#idx_id_14">O</a> <a class="link" href="s13.html#idx_id_15">P</a> <a class="link" href="s13.html#idx_id_16">Q</a> <a class="link" href="s13.html#idx_id_17">R</a> <a class="link" href="s13.html#idx_id_18">S</a> <a class="link" href="s13.html#idx_id_19">T</a> <a class="link" href="s13.html#idx_id_20">U</a> <a class="link" href=
"s13.html#idx_id_21">V</a> <a class="link" href="s13.html#idx_id_22">W</a> <a class="link" href="s13.html#idx_id_23">Z</a></p>
 <div class="variablelist"><dl class="variablelist">
 <dt>
@@ -408,6 +408,7 @@
 <li class="listitem" style="list-style-type: none"><p><a class="link" href="../math_toolkit/extern_c/c99.html" title="C99 C Functions"><span class="index-entry-level-1">C99 C Functions</span></a></p></li>
 </ul></div>
 </li>
+<li class="listitem" style="list-style-type: none"><p><a class="link" href="../math_toolkit/special/powers/cos_pi.html" title="cos_pi"><span class="index-entry-level-0">cos_pi</span></a></p></li>
 <li class="listitem" style="list-style-type: none">
 <p><span class="index-entry-level-0">cyl_bessel_i</span></p>
 <div class="index"><ul class="index" style="list-style-type: none; ">
@@ -2021,6 +2022,7 @@
 </li>
 <li class="listitem" style="list-style-type: none"><p><a class="link" href="../math_toolkit/special/sinc/sinc_pi.html" title="sinc_pi"><span class="index-entry-level-0">sinc_pi</span></a></p></li>
 <li class="listitem" style="list-style-type: none"><p><a class="link" href="../math_toolkit/special/sinc/sinhc_pi.html" title="sinhc_pi"><span class="index-entry-level-0">sinhc_pi</span></a></p></li>
+<li class="listitem" style="list-style-type: none"><p><a class="link" href="../math_toolkit/special/powers/sin_pi.html" title="sin_pi"><span class="index-entry-level-0">sin_pi</span></a></p></li>
 <li class="listitem" style="list-style-type: none">
 <p><span class="index-entry-level-0">skewness</span></p>
 <div class="index"><ul class="index" style="list-style-type: none; ">

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-12-27 03:44:08 EST (Thu, 27 Dec 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 id1351545">
+<div class="section id1341560">
 <div class="titlepage"><div><div><h2 class="title" style="clear: both">
-<a name="id1351545"></a>Class Index</h2></div></div></div>
+<a name="id1341560"></a>Class Index</h2></div></div></div>
 <p><a class="link" href="s14.html#idx_id_25">B</a> <a class="link" href="s14.html#idx_id_26">C</a> <a class="link" href="s14.html#idx_id_27">D</a> <a class="link" href="s14.html#idx_id_28">E</a> <a class="link" href="s14.html#idx_id_29">F</a> <a class="link" href="s14.html#idx_id_30">G</a> <a class="link" href="s14.html#idx_id_31">H</a> <a class="link" href="s14.html#idx_id_32">I</a> <a class="link" href="s14.html#idx_id_35">L</a> <a class="link" href="s14.html#idx_id_36">M</a> <a class="link" href="s14.html#idx_id_37">N</a> <a class="link" href="s14.html#idx_id_39">P</a> <a class="link" href="s14.html#idx_id_41">R</a> <a class="link" href="s14.html#idx_id_42">S</a> <a class="link" href="s14.html#idx_id_43">T</a> <a class="link" href="s14.html#idx_id_44">U</a> <a class="link" href="s14.html#idx_id_46">W</a></p>
 <div class="variablelist"><dl class="variablelist">
 <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-12-27 03:44:08 EST (Thu, 27 Dec 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 id1352486">
+<div class="section id1343048">
 <div class="titlepage"><div><div><h2 class="title" style="clear: both">
-<a name="id1352486"></a>Typedef Index</h2></div></div></div>
+<a name="id1343048"></a>Typedef Index</h2></div></div></div>
 <p><a class="link" href="s15.html#idx_id_48">A</a> <a class="link" href="s15.html#idx_id_49">B</a> <a class="link" href="s15.html#idx_id_50">C</a> <a class="link" href="s15.html#idx_id_51">D</a> <a class="link" href="s15.html#idx_id_52">E</a> <a class="link" href="s15.html#idx_id_53">F</a> <a class="link" href="s15.html#idx_id_54">G</a> <a class="link" href="s15.html#idx_id_55">H</a> <a class="link" href="s15.html#idx_id_56">I</a> <a class="link" href="s15.html#idx_id_59">L</a> <a class="link" href="s15.html#idx_id_61">N</a> <a class="link" href="s15.html#idx_id_62">O</a> <a class="link" href="s15.html#idx_id_63">P</a> <a class="link" href="s15.html#idx_id_65">R</a> <a class="link" href="s15.html#idx_id_66">S</a> <a class="link" href="s15.html#idx_id_67">T</a> <a class="link" href="s15.html#idx_id_68">U</a> <a class="link" href="s15.html#idx_id_69">V</a> <a class="link" href="s15.html#idx_id_70">W</a></p>
 <div class="variablelist"><dl class="variablelist">
 <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-12-27 03:44:08 EST (Thu, 27 Dec 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 id1354986">
+<div class="section id1344864">
 <div class="titlepage"><div><div><h2 class="title" style="clear: both">
-<a name="id1354986"></a>Macro Index</h2></div></div></div>
+<a name="id1344864"></a>Macro Index</h2></div></div></div>
 <p><a class="link" href="s16.html#idx_id_73">B</a> <a class="link" href="s16.html#idx_id_77">F</a></p>
 <div class="variablelist"><dl class="variablelist">
 <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-12-27 03:44:08 EST (Thu, 27 Dec 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 id864598">
+<div class="section id854253">
 <div class="titlepage"><div><div><h2 class="title" style="clear: both">
-<a name="id864598"></a>Index</h2></div></div></div>
+<a name="id854253"></a>Index</h2></div></div></div>
 <p><a class="link" href="s17.html#idx_id_96">A</a> <a class="link" href="s17.html#idx_id_97">B</a> <a class="link" href="s17.html#idx_id_98">C</a> <a class="link" href="s17.html#idx_id_99">D</a> <a class="link" href="s17.html#idx_id_100">E</a> <a class="link" href="s17.html#idx_id_101">F</a> <a class="link" href="s17.html#idx_id_102">G</a> <a class="link" href="s17.html#idx_id_103">H</a> <a class="link" href="s17.html#idx_id_104">I</a> <a class="link" href="s17.html#idx_id_105">J</a> <a class="link" href="s17.html#idx_id_106">K</a> <a class="link" href="s17.html#idx_id_107">L</a> <a class="link" href="s17.html#idx_id_108">M</a> <a class="link" href="s17.html#idx_id_109">N</a> <a class="link" href="s17.html#idx_id_110">O</a> <a class="link" href="s17.html#idx_id_111">P</a> <a class="link" href="s17.html#idx_id_112">Q</a> <a class="link" href="s17.html#idx_id_113">R</a> <a class="link" href="s17.html#idx_id_114">S</a> <a class="link" href="s17.html#idx_id_115">T</a> <a class="link" href="s17.html#idx_id_116">
U</a> <a class="link" href="s17.html#idx_id_117">V</a> <a class="link" href="s17.html#idx_id_118">W</a> <a class="link" href="s17.html#idx_id_119">Z</a></p>
 <div class="variablelist"><dl class="variablelist">
 <dt>
@@ -1374,6 +1374,7 @@
 <li class="listitem" style="list-style-type: none"><p><a class="link" href="../math_toolkit/extern_c/c99.html" title="C99 C Functions"><span class="index-entry-level-1">C99 C Functions</span></a></p></li>
 </ul></div>
 </li>
+<li class="listitem" style="list-style-type: none"><p><a class="link" href="../math_toolkit/special/powers/cos_pi.html" title="cos_pi"><span class="index-entry-level-0">cos_pi</span></a></p></li>
 <li class="listitem" style="list-style-type: none">
 <p><span class="index-entry-level-0">Credits and Acknowledgements</span></p>
 <div class="index"><ul class="index" style="list-style-type: none; ">
@@ -4300,6 +4301,7 @@
 </li>
 <li class="listitem" style="list-style-type: none"><p><a class="link" href="../math_toolkit/special/sinc/sinc_pi.html" title="sinc_pi"><span class="index-entry-level-0">sinc_pi</span></a></p></li>
 <li class="listitem" style="list-style-type: none"><p><a class="link" href="../math_toolkit/special/sinc/sinhc_pi.html" title="sinhc_pi"><span class="index-entry-level-0">sinhc_pi</span></a></p></li>
+<li class="listitem" style="list-style-type: none"><p><a class="link" href="../math_toolkit/special/powers/sin_pi.html" title="sin_pi"><span class="index-entry-level-0">sin_pi</span></a></p></li>
 <li class="listitem" style="list-style-type: none">
 <p><span class="index-entry-level-0">Skew Normal Distribution</span></p>
 <div class="index"><ul class="index" style="list-style-type: none; ">

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/main_overview/conventions.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/main_overview/conventions.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/main_overview/conventions.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -27,7 +27,7 @@
 <a name="math_toolkit.main_overview.conventions"></a><a class="link" href="conventions.html" title="Document Conventions">Document Conventions</a>
 </h3></div></div></div>
 <p>
- <a class="indexterm" name="id857525"></a>
+ <a class="indexterm" name="id846602"></a>
       </p>
 <p>
         This documentation aims to use of the following naming and formatting conventions.

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-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -61,6 +61,10 @@
             and <a class="link" href="../special/sf_gamma/igamma.html" title="Incomplete Gamma Functions">gamma_p</a>/<a class="link" href="../special/sf_gamma/igamma.html" title="Incomplete Gamma Functions">gamma_q</a> thanks
             to suggestions from Rocco Romeo.
           </li>
+<li class="listitem">
+ Improve accuracy of Bessel J and Y for integer orders thanks to suggestions
+ from Rocco Romeo.
+ </li>
 </ul></div>
 <h5>
 <a name="math_toolkit.main_overview.history1.h1"></a>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/main_overview/navigation.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/main_overview/navigation.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/main_overview/navigation.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -27,7 +27,7 @@
 <a name="math_toolkit.main_overview.navigation"></a><a class="link" href="navigation.html" title="Navigation">Navigation</a>
 </h3></div></div></div>
 <p>
- <a class="indexterm" name="id857395"></a>
+ <a class="indexterm" name="id846473"></a>
       </p>
 <p>
         Boost.Math documentation is provided in both HTML and PDF formats.

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -157,9 +157,10 @@
 <dt><span class="section"><a href="special/expint/expint_i.html">Exponential Integral
         Ei</a></span></dt>
 </dl></dd>
-<dt><span class="section"><a href="special/powers.html">Logs, Powers, Roots and
- Exponentials</a></span></dt>
+<dt><span class="section">Basic Functions</span></dt>
 <dd><dl>
+<dt><span class="section">sin_pi</span></dt>
+<dt><span class="section">cos_pi</span></dt>
 <dt><span class="section">log1p</span></dt>
 <dt><span class="section">expm1</span></dt>
 <dt><span class="section">cbrt</span></dt>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/expint/expint_i.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/expint/expint_i.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/expint/expint_i.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -7,7 +7,7 @@
 <link rel="home" href="../../../index.html" title="Math Toolkit">
 <link rel="up" href="../expint.html" title="Exponential Integrals">
 <link rel="prev" href="expint_n.html" title="Exponential Integral En">
-<link rel="next" href="../powers.html" title="Logs, Powers, Roots and Exponentials">
+<link rel="next" href="../powers.html" title="Basic Functions">
 </head>
 <body bgcolor="white" text="black" link="#0000FF" vlink="#840084" alink="#0000FF">
 <table cellpadding="2" width="100%"><tr>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -1,13 +1,13 @@
 <html>
 <head>
 <meta http-equiv="Content-Type" content="text/html; charset=US-ASCII">
-<title>Logs, Powers, Roots and Exponentials</title>
+<title>Basic Functions</title>
 <link rel="stylesheet" href="../../../../../../../doc/src/boostbook.css" type="text/css">
 <meta name="generator" content="DocBook XSL Stylesheets V1.77.1">
 <link rel="home" href="../../index.html" title="Math Toolkit">
 <link rel="up" href="../special.html" title="Special Functions">
 <link rel="prev" href="expint/expint_i.html" title="Exponential Integral Ei">
-<link rel="next" href="powers/log1p.html" title="log1p">
+<link rel="next" href="powers/sin_pi.html" title="sin_pi">
 </head>
 <body bgcolor="white" text="black" link="#0000FF" vlink="#840084" alink="#0000FF">
 <table cellpadding="2" width="100%"><tr>
@@ -20,14 +20,15 @@
 </tr></table>
 <hr>
 <div class="spirit-nav">
-<a accesskey="p" href="expint/expint_i.html"><img src="../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../special.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="powers/log1p.html"><img src="../../../../../../../doc/src/images/next.png" alt="Next"></a>
+<a accesskey="p" href="expint/expint_i.html"><img src="../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../special.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="powers/sin_pi.html"><img src="../../../../../../../doc/src/images/next.png" alt="Next"></a>
 </div>
 <div class="section math_toolkit_special_powers">
 <div class="titlepage"><div><div><h3 class="title">
-<a name="math_toolkit.special.powers"></a><a class="link" href="powers.html" title="Logs, Powers, Roots and Exponentials">Logs, Powers, Roots and
- Exponentials</a>
+<a name="math_toolkit.special.powers"></a><a class="link" href="powers.html" title="Basic Functions">Basic Functions</a>
 </h3></div></div></div>
 <div class="toc"><dl>
+<dt><span class="section">sin_pi</span></dt>
+<dt><span class="section">cos_pi</span></dt>
 <dt><span class="section">log1p</span></dt>
 <dt><span class="section">expm1</span></dt>
 <dt><span class="section">cbrt</span></dt>
@@ -50,7 +51,7 @@
 </tr></table>
 <hr>
 <div class="spirit-nav">
-<a accesskey="p" href="expint/expint_i.html"><img src="../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../special.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="powers/log1p.html"><img src="../../../../../../../doc/src/images/next.png" alt="Next"></a>
+<a accesskey="p" href="expint/expint_i.html"><img src="../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../special.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="powers/sin_pi.html"><img src="../../../../../../../doc/src/images/next.png" alt="Next"></a>
 </div>
 </body>
 </html>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/cbrt.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/cbrt.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/cbrt.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -5,7 +5,7 @@
 <link rel="stylesheet" href="../../../../../../../../doc/src/boostbook.css" type="text/css">
 <meta name="generator" content="DocBook XSL Stylesheets V1.77.1">
 <link rel="home" href="../../../index.html" title="Math Toolkit">
-<link rel="up" href="../powers.html" title="Logs, Powers, Roots and Exponentials">
+<link rel="up" href="../powers.html" title="Basic Functions">
 <link rel="prev" href="expm1.html" title="expm1">
 <link rel="next" href="sqrt1pm1.html" title="sqrt1pm1">
 </head>

Added: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/cos_pi.html
==============================================================================
--- (empty file)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/cos_pi.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -0,0 +1,78 @@
+<html>
+<head>
+<meta http-equiv="Content-Type" content="text/html; charset=US-ASCII">
+<title>cos_pi</title>
+<link rel="stylesheet" href="../../../../../../../../doc/src/boostbook.css" type="text/css">
+<meta name="generator" content="DocBook XSL Stylesheets V1.77.1">
+<link rel="home" href="../../../index.html" title="Math Toolkit">
+<link rel="up" href="../powers.html" title="Basic Functions">
+<link rel="prev" href="sin_pi.html" title="sin_pi">
+<link rel="next" href="log1p.html" title="log1p">
+</head>
+<body bgcolor="white" text="black" link="#0000FF" vlink="#840084" alink="#0000FF">
+<table cellpadding="2" width="100%"><tr>
+<td valign="top"><img alt="Boost C++ Libraries" width="277" height="86" src="../../../../../../../../boost.png"></td>
+<td align="center">Home</td>
+<td align="center">Libraries</td>
+<td align="center">People</td>
+<td align="center">FAQ</td>
+<td align="center">More</td>
+</tr></table>
+<hr>
+<div class="spirit-nav">
+<a accesskey="p" href="sin_pi.html"><img src="../../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../powers.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="log1p.html"><img src="../../../../../../../../doc/src/images/next.png" alt="Next"></a>
+</div>
+<div class="section math_toolkit_special_powers_cos_pi">
+<div class="titlepage"><div><div><h4 class="title">
+<a name="math_toolkit.special.powers.cos_pi"></a><a class="link" href="cos_pi.html" title="cos_pi">cos_pi</a>
+</h4></div></div></div>
+<p>
+</p>
+<pre class="programlisting"><span class="preprocessor">#include</span> <span class="special">&lt;</span><span class="identifier">boost</span><span class="special">/</span><span class="identifier">math</span><span class="special">/</span><span class="identifier">special_functions</span><span class="special">/</span><span class="identifier">cos_pi</span><span class="special">.</span><span class="identifier">hpp</span><span class="special">&gt;</span>
+</pre>
+<p>
+ </p>
+<pre class="programlisting"><span class="keyword">namespace</span> <span class="identifier">boost</span><span class="special">{</span> <span class="keyword">namespace</span> <span class="identifier">math</span><span class="special">{</span>
+
+<span class="keyword">template</span> <span class="special">&lt;</span><span class="keyword">class</span> <span class="identifier">T</span><span class="special">&gt;</span>
+<a class="link" href="../../main_overview/result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">cos_pi</span><span class="special">(</span><span class="identifier">T</span> <span class="identifier">x</span><span class="special">);</span>
+
+<span class="keyword">template</span> <span class="special">&lt;</span><span class="keyword">class</span> <span class="identifier">T</span><span class="special">,</span> <span class="keyword">class</span> <a class="link" href="../../policy.html" title="Policies">Policy</a><span class="special">&gt;</span>
+<a class="link" href="../../main_overview/result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">cos_pi</span><span class="special">(</span><span class="identifier">T</span> <span class="identifier">x</span><span class="special">,</span> <span class="keyword">const</span> <a class="link" href="../../policy.html" title="Policies">Policy</a><span class="special">&amp;);</span>
+
+<span class="special">}}</span> <span class="comment">// namespaces</span>
+</pre>
+<p>
+ Returns the cosine of <span class="emphasis"><em>&#960;x</em></span>.
+ </p>
+<p>
+ The return type of this function is computed using the <a class="link" href="../../main_overview/result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>result
+ type calculation rules</em></span></a>: the return is <code class="computeroutput"><span class="keyword">double</span></code>
+ when <span class="emphasis"><em>x</em></span> is an integer type and T otherwise.
+ </p>
+<p>
+ The final <a class="link" href="../../policy.html" title="Policies">Policy</a> argument is
+ optional and can be used to control the behaviour of the function: how
+ it handles errors, what level of precision to use etc. Refer to the <a class="link" href="../../policy.html" title="Policies">policy documentation for more details</a>.
+ </p>
+<p>
+ This function performs exact all-integer arithmetic argument reduction
+ before computing the cosine of <span class="emphasis"><em>&#960;x</em></span>.
+ </p>
+</div>
+<table xmlns:rev="http://www.cs.rpi.edu/~gregod/boost/tools/doc/revision" width="100%"><tr>
+<td align="left"></td>
+<td align="right"><div class="copyright-footer">Copyright &#169; 2006-2010, 2012 John Maddock, Paul A. Bristow, Hubert Holin, Xiaogang
+ Zhang, Bruno Lalande, Johan R&#229;de, Gautam Sewani, Thijs van den Berg and Benjamin
+ Sobotta<p>
+ Distributed under the Boost Software License, Version 1.0. (See accompanying
+ file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
+ </p>
+</div></td>
+</tr></table>
+<hr>
+<div class="spirit-nav">
+<a accesskey="p" href="sin_pi.html"><img src="../../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../powers.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="log1p.html"><img src="../../../../../../../../doc/src/images/next.png" alt="Next"></a>
+</div>
+</body>
+</html>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/ct_pow.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/ct_pow.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/ct_pow.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -5,7 +5,7 @@
 <link rel="stylesheet" href="../../../../../../../../doc/src/boostbook.css" type="text/css">
 <meta name="generator" content="DocBook XSL Stylesheets V1.77.1">
 <link rel="home" href="../../../index.html" title="Math Toolkit">
-<link rel="up" href="../powers.html" title="Logs, Powers, Roots and Exponentials">
+<link rel="up" href="../powers.html" title="Basic Functions">
 <link rel="prev" href="hypot.html" title="hypot">
 <link rel="next" href="../sinc.html" title="Sinus Cardinal and Hyperbolic Sinus Cardinal Functions">
 </head>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/expm1.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/expm1.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/expm1.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -5,7 +5,7 @@
 <link rel="stylesheet" href="../../../../../../../../doc/src/boostbook.css" type="text/css">
 <meta name="generator" content="DocBook XSL Stylesheets V1.77.1">
 <link rel="home" href="../../../index.html" title="Math Toolkit">
-<link rel="up" href="../powers.html" title="Logs, Powers, Roots and Exponentials">
+<link rel="up" href="../powers.html" title="Basic Functions">
 <link rel="prev" href="log1p.html" title="log1p">
 <link rel="next" href="cbrt.html" title="cbrt">
 </head>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/hypot.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/hypot.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/hypot.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -5,7 +5,7 @@
 <link rel="stylesheet" href="../../../../../../../../doc/src/boostbook.css" type="text/css">
 <meta name="generator" content="DocBook XSL Stylesheets V1.77.1">
 <link rel="home" href="../../../index.html" title="Math Toolkit">
-<link rel="up" href="../powers.html" title="Logs, Powers, Roots and Exponentials">
+<link rel="up" href="../powers.html" title="Basic Functions">
 <link rel="prev" href="powm1.html" title="powm1">
 <link rel="next" href="ct_pow.html" title="Compile Time Power of a Runtime Base">
 </head>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/log1p.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/log1p.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/log1p.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -5,8 +5,8 @@
 <link rel="stylesheet" href="../../../../../../../../doc/src/boostbook.css" type="text/css">
 <meta name="generator" content="DocBook XSL Stylesheets V1.77.1">
 <link rel="home" href="../../../index.html" title="Math Toolkit">
-<link rel="up" href="../powers.html" title="Logs, Powers, Roots and Exponentials">
-<link rel="prev" href="../powers.html" title="Logs, Powers, Roots and Exponentials">
+<link rel="up" href="../powers.html" title="Basic Functions">
+<link rel="prev" href="cos_pi.html" title="cos_pi">
 <link rel="next" href="expm1.html" title="expm1">
 </head>
 <body bgcolor="white" text="black" link="#0000FF" vlink="#840084" alink="#0000FF">
@@ -20,7 +20,7 @@
 </tr></table>
 <hr>
 <div class="spirit-nav">
-<a accesskey="p" href="../powers.html"><img src="../../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../powers.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="expm1.html"><img src="../../../../../../../../doc/src/images/next.png" alt="Next"></a>
+<a accesskey="p" href="cos_pi.html"><img src="../../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../powers.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="expm1.html"><img src="../../../../../../../../doc/src/images/next.png" alt="Next"></a>
 </div>
 <div class="section math_toolkit_special_powers_log1p">
 <div class="titlepage"><div><div><h4 class="title">
@@ -115,7 +115,7 @@
 </tr></table>
 <hr>
 <div class="spirit-nav">
-<a accesskey="p" href="../powers.html"><img src="../../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../powers.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="expm1.html"><img src="../../../../../../../../doc/src/images/next.png" alt="Next"></a>
+<a accesskey="p" href="cos_pi.html"><img src="../../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../powers.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="expm1.html"><img src="../../../../../../../../doc/src/images/next.png" alt="Next"></a>
 </div>
 </body>
 </html>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/powm1.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/powm1.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/powm1.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -5,7 +5,7 @@
 <link rel="stylesheet" href="../../../../../../../../doc/src/boostbook.css" type="text/css">
 <meta name="generator" content="DocBook XSL Stylesheets V1.77.1">
 <link rel="home" href="../../../index.html" title="Math Toolkit">
-<link rel="up" href="../powers.html" title="Logs, Powers, Roots and Exponentials">
+<link rel="up" href="../powers.html" title="Basic Functions">
 <link rel="prev" href="sqrt1pm1.html" title="sqrt1pm1">
 <link rel="next" href="hypot.html" title="hypot">
 </head>

Added: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/sin_pi.html
==============================================================================
--- (empty file)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/sin_pi.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -0,0 +1,78 @@
+<html>
+<head>
+<meta http-equiv="Content-Type" content="text/html; charset=US-ASCII">
+<title>sin_pi</title>
+<link rel="stylesheet" href="../../../../../../../../doc/src/boostbook.css" type="text/css">
+<meta name="generator" content="DocBook XSL Stylesheets V1.77.1">
+<link rel="home" href="../../../index.html" title="Math Toolkit">
+<link rel="up" href="../powers.html" title="Basic Functions">
+<link rel="prev" href="../powers.html" title="Basic Functions">
+<link rel="next" href="cos_pi.html" title="cos_pi">
+</head>
+<body bgcolor="white" text="black" link="#0000FF" vlink="#840084" alink="#0000FF">
+<table cellpadding="2" width="100%"><tr>
+<td valign="top"><img alt="Boost C++ Libraries" width="277" height="86" src="../../../../../../../../boost.png"></td>
+<td align="center">Home</td>
+<td align="center">Libraries</td>
+<td align="center">People</td>
+<td align="center">FAQ</td>
+<td align="center">More</td>
+</tr></table>
+<hr>
+<div class="spirit-nav">
+<a accesskey="p" href="../powers.html"><img src="../../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../powers.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="cos_pi.html"><img src="../../../../../../../../doc/src/images/next.png" alt="Next"></a>
+</div>
+<div class="section math_toolkit_special_powers_sin_pi">
+<div class="titlepage"><div><div><h4 class="title">
+<a name="math_toolkit.special.powers.sin_pi"></a><a class="link" href="sin_pi.html" title="sin_pi">sin_pi</a>
+</h4></div></div></div>
+<p>
+</p>
+<pre class="programlisting"><span class="preprocessor">#include</span> <span class="special">&lt;</span><span class="identifier">boost</span><span class="special">/</span><span class="identifier">math</span><span class="special">/</span><span class="identifier">special_functions</span><span class="special">/</span><span class="identifier">sin_pi</span><span class="special">.</span><span class="identifier">hpp</span><span class="special">&gt;</span>
+</pre>
+<p>
+ </p>
+<pre class="programlisting"><span class="keyword">namespace</span> <span class="identifier">boost</span><span class="special">{</span> <span class="keyword">namespace</span> <span class="identifier">math</span><span class="special">{</span>
+
+<span class="keyword">template</span> <span class="special">&lt;</span><span class="keyword">class</span> <span class="identifier">T</span><span class="special">&gt;</span>
+<a class="link" href="../../main_overview/result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">sin_pi</span><span class="special">(</span><span class="identifier">T</span> <span class="identifier">x</span><span class="special">);</span>
+
+<span class="keyword">template</span> <span class="special">&lt;</span><span class="keyword">class</span> <span class="identifier">T</span><span class="special">,</span> <span class="keyword">class</span> <a class="link" href="../../policy.html" title="Policies">Policy</a><span class="special">&gt;</span>
+<a class="link" href="../../main_overview/result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">sin_pi</span><span class="special">(</span><span class="identifier">T</span> <span class="identifier">x</span><span class="special">,</span> <span class="keyword">const</span> <a class="link" href="../../policy.html" title="Policies">Policy</a><span class="special">&amp;);</span>
+
+<span class="special">}}</span> <span class="comment">// namespaces</span>
+</pre>
+<p>
+ Returns the sine of <span class="emphasis"><em>&#960;x</em></span>.
+ </p>
+<p>
+ The return type of this function is computed using the <a class="link" href="../../main_overview/result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>result
+ type calculation rules</em></span></a>: the return is <code class="computeroutput"><span class="keyword">double</span></code>
+ when <span class="emphasis"><em>x</em></span> is an integer type and T otherwise.
+ </p>
+<p>
+ The final <a class="link" href="../../policy.html" title="Policies">Policy</a> argument is
+ optional and can be used to control the behaviour of the function: how
+ it handles errors, what level of precision to use etc. Refer to the <a class="link" href="../../policy.html" title="Policies">policy documentation for more details</a>.
+ </p>
+<p>
+ This function performs exact all-integer arithmetic argument reduction
+ before computing the sine of <span class="emphasis"><em>&#960;x</em></span>.
+ </p>
+</div>
+<table xmlns:rev="http://www.cs.rpi.edu/~gregod/boost/tools/doc/revision" width="100%"><tr>
+<td align="left"></td>
+<td align="right"><div class="copyright-footer">Copyright &#169; 2006-2010, 2012 John Maddock, Paul A. Bristow, Hubert Holin, Xiaogang
+ Zhang, Bruno Lalande, Johan R&#229;de, Gautam Sewani, Thijs van den Berg and Benjamin
+ Sobotta<p>
+ Distributed under the Boost Software License, Version 1.0. (See accompanying
+ file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
+ </p>
+</div></td>
+</tr></table>
+<hr>
+<div class="spirit-nav">
+<a accesskey="p" href="../powers.html"><img src="../../../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../powers.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="cos_pi.html"><img src="../../../../../../../../doc/src/images/next.png" alt="Next"></a>
+</div>
+</body>
+</html>

Modified: trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/sqrt1pm1.html
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/sqrt1pm1.html (original)
+++ trunk/libs/math/doc/sf_and_dist/html/math_toolkit/special/powers/sqrt1pm1.html 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -5,7 +5,7 @@
 <link rel="stylesheet" href="../../../../../../../../doc/src/boostbook.css" type="text/css">
 <meta name="generator" content="DocBook XSL Stylesheets V1.77.1">
 <link rel="home" href="../../../index.html" title="Math Toolkit">
-<link rel="up" href="../powers.html" title="Logs, Powers, Roots and Exponentials">
+<link rel="up" href="../powers.html" title="Basic Functions">
 <link rel="prev" href="cbrt.html" title="cbrt">
 <link rel="next" href="powm1.html" title="powm1">
 </head>

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-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -60,6 +60,10 @@
             and <a class="link" href="../special/sf_gamma/igamma.html" title="Incomplete Gamma Functions">gamma_p</a>/<a class="link" href="../special/sf_gamma/igamma.html" title="Incomplete Gamma Functions">gamma_q</a> thanks
             to suggestions from Rocco Romeo.
           </li>
+<li class="listitem">
+ Improve accuracy of Bessel J and Y for integer orders thanks to suggestions
+ from Rocco Romeo.
+ </li>
 </ul></div>
 <h5>
 <a name="math_toolkit.status.history1.h1"></a>

Modified: trunk/libs/math/doc/sf_and_dist/powers.qbk
==============================================================================
--- trunk/libs/math/doc/sf_and_dist/powers.qbk (original)
+++ trunk/libs/math/doc/sf_and_dist/powers.qbk 2012-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -1,4 +1,58 @@
-[section:powers Logs, Powers, Roots and Exponentials]
+[section:powers Basic Functions]
+
+[section:sin_pi sin_pi]
+
+``
+#include <boost/math/special_functions/sin_pi.hpp>
+``
+
+ namespace boost{ namespace math{
+
+ template <class T>
+ ``__sf_result`` sin_pi(T x);
+
+ template <class T, class ``__Policy``>
+ ``__sf_result`` sin_pi(T x, const ``__Policy``&);
+
+ }} // namespaces
+
+Returns the sine of ['[pi]x].
+
+The return type of this function is computed using the __arg_pomotion_rules:
+the return is `double` when /x/ is an integer type and T otherwise.
+
+[optional_policy]
+
+This function performs exact all-integer arithmetic argument reduction before computing the sine of ['[pi]x].
+
+[endsect]
+
+[section:cos_pi cos_pi]
+
+``
+#include <boost/math/special_functions/cos_pi.hpp>
+``
+
+ namespace boost{ namespace math{
+
+ template <class T>
+ ``__sf_result`` cos_pi(T x);
+
+ template <class T, class ``__Policy``>
+ ``__sf_result`` cos_pi(T x, const ``__Policy``&);
+
+ }} // namespaces
+
+Returns the cosine of ['[pi]x].
+
+The return type of this function is computed using the __arg_pomotion_rules:
+the return is `double` when /x/ is an integer type and T otherwise.
+
+[optional_policy]
+
+This function performs exact all-integer arithmetic argument reduction before computing the cosine of ['[pi]x].
+
+[endsect]
 
 [section:log1p log1p]
 

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-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -254,14 +254,7 @@
       ".*", // stdlib
       ".*", // platform
       largest_type, // test type(s)
- ".*JN.*Integer.*", // test data group
- ".*", 30000, 10000); // test function
- add_expected_result(
- ".*", // compiler
- ".*", // stdlib
- ".*", // platform
- largest_type, // test type(s)
- ".*(JN|j).*|.*Tricky.*", // test data group
+ "Bessel j:.*|Bessel JN: Mathworld.*|.*Tricky.*", // test data group
       ".*", 1500, 700); // test function
    add_expected_result(
       ".*", // compiler

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-12-27 03:44:08 EST (Thu, 27 Dec 2012)
@@ -191,13 +191,6 @@
       ".*", // compiler
       ".*", // stdlib
       ".*", // platform
- largest_type, // test type(s)
- ".*Y0.*Random.*", // test data group
- ".*", 600, 400); // test function
- add_expected_result(
- ".*", // compiler
- ".*", // stdlib
- ".*", // platform
       "real_concept", // test type(s)
       ".*(Y[nv]|y).*Random.*", // test data group
       ".*", 2000, 2000); // test function


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