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.