Re: [Boost-bugs] [Boost C++ Libraries] #8621: erf function gives wrong results with pgcpp - PGI 10.4

Subject: Re: [Boost-bugs] [Boost C++ Libraries] #8621: erf function gives wrong results with pgcpp - PGI 10.4
From: Boost C++ Libraries (noreply_at_[hidden])
Date: 2013-06-13 17:41:47


#8621: erf function gives wrong results with pgcpp - PGI 10.4
-------------------------------------------------+-------------------------
  Reporter: Jacques Desfosses | Owner: johnmaddock
  <jacques.desfosses@…> | Status: new
      Type: Bugs | Component: math
 Milestone: To Be Determined | Severity: Problem
   Version: Boost 1.53.0 | Keywords: erf pgi
Resolution: |
-------------------------------------------------+-------------------------

Comment (by Jacques Desfosses <jacques.desfosses@…>):

 Hi John. Sorry for the delay. The patch works well. Did you, by any
 chance, forget to modify the following arrays?

 '''boost/math/special_functions/detail/lgamma_small.hpp'''
 {{{

 static const T Q[] = {
    static_cast<T>(0.1e1),
    static_cast<T>(BOOST_MATH_BIG_CONSTANT(T, 64,
 0.196202987197795200688e1)),
    static_cast<T>(BOOST_MATH_BIG_CONSTANT(T, 64,
 0.148019669424231326694e1)),
    static_cast<T>(BOOST_MATH_BIG_CONSTANT(T, 64,
 0.541391432071720958364e0)),
    static_cast<T>(BOOST_MATH_BIG_CONSTANT(T, 64,
 0.988504251128010129477e-1)),
    static_cast<T>(BOOST_MATH_BIG_CONSTANT(T, 64,
 0.82130967464889339326e-2)),
    static_cast<T>(BOOST_MATH_BIG_CONSTANT(T, 64,
 0.224936291922115757597e-3)),
    static_cast<T>(BOOST_MATH_BIG_CONSTANT(T, 64,
 -0.223352763208617092964e-6))
 };


 static const T Q[] = {
    static_cast<T>(0.1e1),
    static_cast<T>(BOOST_MATH_BIG_CONSTANT(T, 64,
 -0.150169356054485044494e1)),
    static_cast<T>(BOOST_MATH_BIG_CONSTANT(T, 64,
 0.846973248876495016101e0)),
    static_cast<T>(BOOST_MATH_BIG_CONSTANT(T, 64,
 -0.220095151814995745555e0)),
    static_cast<T>(BOOST_MATH_BIG_CONSTANT(T, 64,
 0.25582797155975869989e-1)),
    static_cast<T>(BOOST_MATH_BIG_CONSTANT(T, 64,
 -0.100666795539143372762e-2)),
    static_cast<T>(BOOST_MATH_BIG_CONSTANT(T, 64,
 -0.827193521891290553639e-6))
 };

 }}}

 '''boost/math/special_functions/zeta.hpp'''
 {{{
 static const T Q[8] = {
    1,
    BOOST_MATH_BIG_CONSTANT(T, 64, 0.286577739726542730421),
    BOOST_MATH_BIG_CONSTANT(T, 64, 0.0447355811517733225843),
    BOOST_MATH_BIG_CONSTANT(T, 64, 0.00430125107610252363302),
    BOOST_MATH_BIG_CONSTANT(T, 64, 0.000284956969089786662045),
    BOOST_MATH_BIG_CONSTANT(T, 64, 0.116188101609848411329e-4),
    BOOST_MATH_BIG_CONSTANT(T, 64, 0.278090318191657278204e-6),
    BOOST_MATH_BIG_CONSTANT(T, 64, -0.19683620233222028478e-8),
 };
 }}}

 I ran into another issue with the PGI compiler and the generic quantile
 finder. If it turns out to be an issue with boost, I will log another
 ticket. Thanks so much for your help,

 Jacques

-- 
Ticket URL: <https://svn.boost.org/trac/boost/ticket/8621#comment:12>
Boost C++ Libraries <http://www.boost.org/>
Boost provides free peer-reviewed portable C++ source libraries.

This archive was generated by hypermail 2.1.7 : 2017-02-16 18:50:13 UTC