boost/math
2
0
Fork 0
mirror of https://github.com/boostorg/math.git synced 2026-01-19 04:22:09 +00:00
Code Issues Projects Releases Wiki Activity

Up precision of tests.

Browse Source 
Fix some conceptual compiler errors in hyperexponential_distribution.
This commit is contained in:
jzmaddock
2014-08-24 09:29:11 +01:00
parent 4fbfc1f2ce
commit ffc483f907
2 changed files with 5 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
include/boost/math/distributions/hyperexponential.hpp
View File

@@ -157,11 +157,12 @@ bool check_rates(char const* function, std::vector<RealT> const& rates, RealT* p
template <typename RealT, typename PolicyT>
bool check_dist(char const* function, std::vector<RealT> const& probabilities, std::vector<RealT> const& rates, RealT* presult, PolicyT const& pol)
{
BOOST_MATH_STD_USING
if (probabilities.size() != rates.size())
{
*presult = policies::raise_domain_error<RealT>(function,
"The parameters \"probabilities\" and \"rates\" must have the same length, but their size differ by: %1%.",
std::abs(static_cast<RealT>(probabilities.size())-static_cast<RealT>(rates.size())),
fabs(static_cast<RealT>(probabilities.size())-static_cast<RealT>(rates.size())),
pol);
return false;
}
@@ -347,6 +348,7 @@ std::pair<RealT,RealT> support(hyperexponential_distribution<RealT,PolicyT> cons
template <typename RealT, typename PolicyT>
RealT pdf(hyperexponential_distribution<RealT, PolicyT> const& dist, RealT const& x)
{
BOOST_MATH_STD_USING
RealT result = 0;
if (!hyperexp_detail::check_x("boost::math::pdf(const boost::math::hyperexponential_distribution<%1%>&, %1%)", x, &result, PolicyT()))
@@ -363,7 +365,7 @@ RealT pdf(hyperexponential_distribution<RealT, PolicyT> const& dist, RealT const
//const exponential_distribution<RealT,PolicyT> exp(rates[i]);
//result += probs[i]*pdf(exp, x);
result += probs[i]*rates[i]*std::exp(-rates[i]*x);
result += probs[i]*rates[i]*exp(-rates[i]*x);
}
return result;

2
test/test_hyperexponential_dist.cpp
View File

@@ -56,7 +56,7 @@ RealT make_tolerance()
}
*/
const RealT tol = 1e-4;
const RealT tol = boost::math::tools::epsilon<RealT>() * 100 * 100;
//std::cout << "[" << __func__ << "] Tolerance: " << tol << "%" << std::endl;