1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
|
// Copyright 2025 Matt Borland
// Use, modification and distribution are subject to 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)
//
// See: https://github.com/boostorg/math/issues/1294
#include <boost/math/distributions/logistic.hpp>
#include "math_unit_test.hpp"
int main()
{
using namespace boost::math::policies;
using boost::math::logistic_distribution;
typedef policy<
promote_float<true>,
promote_double<true>
> with_promotion;
constexpr double p = 2049.0/4096;
constexpr double ref = 9.76562577610225755e-04;
logistic_distribution<double, with_promotion> dist_promote;
const double x = quantile(dist_promote, p);
// Previously we had: 9.76562577610170027e-04
// Which is an ULP distance of 256
CHECK_ULP_CLOSE(x, ref, 1);
return boost::math::test::report_errors();
}
|
