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-09 03:37:52

#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@…>):

 Thanks again for your help! The compiler is okay with long initialization.
 The compiler seems to dislike the mixing of 1L and BOOST_MATH_BIG_CONSTANT
 initializers. I was able to obtain correct results by defining Q using
 these 2 methods:

 1) Without BOOST_MATH_BIG_CONSTANT

 {{{
 static const T Q[] = {
             1L,
             1.84759070983002217845,
             1.42628004845511324508,
             0.578052804889902404909,
             0.12385097467900864233,
             0.0113385233577001411017,
             0.337511472483094676155e-5,
          };
 }}}


 2) Only with BOOST_MATH_BIG_CONSTANT

 {{{
 static const T Q[] = {
             BOOST_MATH_BIG_CONSTANT(T, 53, 1),
             BOOST_MATH_BIG_CONSTANT(T, 53, 1.84759070983002217845),
             BOOST_MATH_BIG_CONSTANT(T, 53, 1.42628004845511324508),
             BOOST_MATH_BIG_CONSTANT(T, 53, 0.578052804889902404909),
             BOOST_MATH_BIG_CONSTANT(T, 53, 0.12385097467900864233),
             BOOST_MATH_BIG_CONSTANT(T, 53, 0.0113385233577001411017),
             BOOST_MATH_BIG_CONSTANT(T, 53, 0.337511472483094676155e-5),
          };
 }}}

 I also had the opportunity to try out an older PGI compiler (version
 7.1-3) and I also ran into the same problem. So the issue is not specific
 to PGI version 10.4.

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