boost/math
2
0
Fork 0
mirror of https://github.com/boostorg/math.git synced 2026-02-13 12:32:15 +00:00
Code Issues Projects Releases Wiki Activity

More review from codereview.stackexchange.com

Browse Source
This commit is contained in:
Nick Thompson
2018-01-12 22:59:56 -06:00
parent a90aac3656
commit 7341a149de
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
include/boost/math/quadrature/naive_monte_carlo.hpp
View File

@@ -130,11 +130,11 @@ public:
{
for (size_t j = 0; j < m_lbs.size(); ++j)
{
x[j] = (gen()+1)*inv_denom;
while (x[j] < numeric_limits<Real>::epsilon() || x[j] > 1 - numeric_limits<Real>::epsilon())
do
{
x[j] = (gen()+1)*inv_denom;
}
while (x[j] < numeric_limits<Real>::epsilon() || x[j] > 1 - numeric_limits<Real>::epsilon());
}
Real y = m_integrand(x);
m_thread_averages.emplace(i, y);
@@ -308,11 +308,11 @@ private:
{
for (size_t i = 0; i < m_lbs.size(); ++i)
{
x[i] = (gen()+1)*inv_denom;
while (x[i] < numeric_limits<Real>::epsilon() || x[i] > 1 - numeric_limits<Real>::epsilon())
do
{
x[i] = (gen()+1)*inv_denom;
}
while (x[i] < numeric_limits<Real>::epsilon() || x[i] > 1 - numeric_limits<Real>::epsilon());
}
Real f = m_integrand(x);
++k;