https://svn.boost.org/trac10/ticket/9042 <p> In boost/math/distributions/normal.hpp: </p> <pre class="wiki">template &lt;class RealType, class Policy&gt; inline RealType cdf(const complemented2_type&lt;normal_distribution&lt;RealType, Policy&gt;, RealType&gt;&amp; c) { BOOST_MATH_STD_USING // for ADL of std functions RealType sd = c.dist.standard_deviation(); RealType mean = c.dist.mean(); RealType x = c.param; static const char* function = "boost::math::cdf(const complement(normal_distribution&lt;%1%&gt;&amp;), %1%)"; if((boost::math::isinf)(x)) { if(x &lt; 0) return 1; // cdf complement -infinity is unity. return 0; // cdf complement +infinity is zero } // These produce MSVC 4127 warnings, so the above used instead. //if(std::numeric_limits&lt;RealType&gt;::has_infinity &amp;&amp; x == std::numeric_limits&lt;RealType&gt;::infinity()) //{ // cdf complement +infinity is zero. // return 0; //} //if(std::numeric_limits&lt;RealType&gt;::has_infinity &amp;&amp; x == -std::numeric_limits&lt;RealType&gt;::infinity()) //{ // cdf complement -infinity is unity. // return 1; //} RealType result = 0; if(false == detail::check_scale(function, sd, &amp;result, Policy())) return result; if(false == detail::check_location(function, mean, &amp;result, Policy())) return result; if(false == detail::check_x(function, x, &amp;result, Policy())) return result; RealType diff = (x - mean) / (sd * constants::root_two&lt;RealType&gt;()); result = boost::math::erfc(diff, Policy()) / 2; return result; } // cdf complement </pre><p> The test if((boost::math::isinf)(x)) occurs before check_scale() or check_location() is called, returning 1 or 0 for infinite x, even if sd or mean are invalid. </p> <p> The calls to check_scale() and check_location() should be moved up just before the test if((boost::math::isinf)(x)). </p> <p> See, for example, how the tests are ordered here: </p> <pre class="wiki">template &lt;class RealType, class Policy&gt; inline RealType cdf(const normal_distribution&lt;RealType, Policy&gt;&amp; dist, const RealType&amp; x) { BOOST_MATH_STD_USING // for ADL of std functions RealType sd = dist.standard_deviation(); RealType mean = dist.mean(); static const char* function = "boost::math::cdf(const normal_distribution&lt;%1%&gt;&amp;, %1%)"; RealType result = 0; if(false == detail::check_scale(function, sd, &amp;result, Policy())) { return result; } if(false == detail::check_location(function, mean, &amp;result, Policy())) { return result; } if((boost::math::isinf)(x)) { if(x &lt; 0) return 0; // -infinity return 1; // + infinity } // These produce MSVC 4127 warnings, so the above used instead. //if(std::numeric_limits&lt;RealType&gt;::has_infinity &amp;&amp; x == std::numeric_limits&lt;RealType&gt;::infinity()) //{ // cdf +infinity is unity. // return 1; //} //if(std::numeric_limits&lt;RealType&gt;::has_infinity &amp;&amp; x == -std::numeric_limits&lt;RealType&gt;::infinity()) //{ // cdf -infinity is zero. // return 0; //} if(false == detail::check_x(function, x, &amp;result, Policy())) { return result; } RealType diff = (x - mean) / (sd * constants::root_two&lt;RealType&gt;()); result = boost::math::erfc(-diff, Policy()) / 2; return result; } // cdf </pre> en-us /htdocs/site/boost.png https://svn.boost.org/trac10/ticket/9042 Trac 1.4.3 John Maddock Sun, 08 Sep 2013 08:46:19 GMT https://svn.boost.org/trac10/ticket/9042#comment:1 https://svn.boost.org/trac10/ticket/9042#comment:1 <ul> <li><strong>status</strong> <span class="trac-field-old">new</span> → <span class="trac-field-new">closed</span> </li> <li><strong>resolution</strong> → <span class="trac-field-new">fixed</span> </li> </ul> <p> (In <a class="changeset" href="https://svn.boost.org/trac10/changeset/85601" title="Fix order of error checks. Fixes #9042.">[85601]</a>) Fix order of error checks. Fixes <a class="closed ticket" href="https://svn.boost.org/trac10/ticket/9042" title="#9042: Bugs: boost::math, cdf(complement(normal_distribution&lt;&gt;(...)) fails to catch ... (closed: fixed)">#9042</a>. </p> Ticket John Maddock Sat, 28 Sep 2013 16:19:56 GMT