[Boost-commit] svn:boost r67212 - trunk/boost/math/distributions

Subject: [Boost-commit] svn:boost r67212 - trunk/boost/math/distributions
From: pbristow_at_[hidden]
Date: 2010-12-13 11:36:02

Author: pbristow
Date: 2010-12-13 11:35:59 EST (Mon, 13 Dec 2010)
New Revision: 67212
URL: http://svn.boost.org/trac/boost/changeset/67212

Log:
use raise_overflow error in place of return NaN
Text files modified: 
   trunk/boost/math/distributions/laplace.hpp |    30 +++++++++++++++++-------------          
   1 files changed, 17 insertions(+), 13 deletions(-)

Modified: trunk/boost/math/distributions/laplace.hpp
==============================================================================
--- trunk/boost/math/distributions/laplace.hpp	(original)
+++ trunk/boost/math/distributions/laplace.hpp	2010-12-13 11:35:59 EST (Mon, 13 Dec 2010)
@@ -66,23 +66,17 @@
          return true;
    }
 
-
 private:
    RealType m_location;
    RealType m_scale;
-
 }; // class laplace_distribution
 
-
-
-//
-// Convenient type synonym
 //
+// Convenient type synonym for double
 typedef laplace_distribution<double> laplace;
 
 //
 // Non member functions
-//
 template <class RealType, class Policy>
 inline const std::pair<RealType, RealType> range(const laplace_distribution<RealType, Policy>&)
 {
@@ -137,7 +131,7 @@
    if (false == dist.check_parameters(function, &result)) return result;
    if (false == detail::check_x(function, x, &result, Policy())) return result;
 
-   // Special cdf values
+   // Special cdf values:
    if((boost::math::isinf)(x))
    {
      if(x < 0) return 0; // -infinity
@@ -160,7 +154,7 @@
 template <class RealType, class Policy>
 inline RealType quantile(const laplace_distribution<RealType, Policy>& dist, const RealType& p)
 {
-   BOOST_MATH_STD_USING // for ADL of std functions
+   BOOST_MATH_STD_USING // for ADL of std functions.
 
    // Checking function argument
    RealType result;
@@ -168,10 +162,20 @@
    if (false == dist.check_parameters(function, &result)) return result;
    if(false == detail::check_probability(function, p, &result, Policy())) return result;
 
-   // extreme values
-   if(p == 0) return -std::numeric_limits<RealType>::infinity();
-   if(p == 1) return std::numeric_limits<RealType>::infinity();
-
+   // Extreme values of p:
+   if(p == 0)
+   {
+      result = policies::raise_overflow_error<RealType>(function,
+        "probability parameter is 0, but must be > 0!", Policy());
+      return -result; // -std::numeric_limits<RealType>::infinity();
+   }
+  
+   if(p == 1)
+   {
+      result = policies::raise_overflow_error<RealType>(function,
+        "probability parameter is 1, but must be < 1!", Policy());
+      return result; // std::numeric_limits<RealType>::infinity();
+   }
    // Calculate Quantile
    RealType scale( dist.scale() );
    RealType location( dist.location() );