Re: [Boost-bugs] [Boost C++ Libraries] #5095: incorrect results from hypergeometric pdf

Subject: Re: [Boost-bugs] [Boost C++ Libraries] #5095: incorrect results from hypergeometric pdf
From: Boost C++ Libraries (noreply_at_[hidden])
Date: 2011-01-20 17:07:59

#5095: incorrect results from hypergeometric pdf
------------------------------------------+---------------------------------
  Reporter:  David Koes <dkoes@…>         |       Owner:  johnmaddock
      Type:  Bugs                         |      Status:  new        
 Milestone:  To Be Determined             |   Component:  math       
   Version:  Boost 1.45.0                 |    Severity:  Problem    
Resolution:                               |    Keywords:             
------------------------------------------+---------------------------------

Comment (by David Koes <dkoes@…>):

 Attached.  Here are some interesting results that may partially explain
 the reproducibility problem:

 {{{
 dkoes_at_quasar:~/tmp$ g++  hyper.cpp; ./a.out
 1.11519e-09
 dkoes_at_quasar:~/tmp$ g++ -m32 hyper.cpp; ./a.out
 0.00639305
 dkoes_at_quasar:~/tmp$ g++ -mfpmath=387 hyper.cpp; ./a.out
 0.00639305
 dkoes_at_quasar:~/tmp$ g++ -mfpmath=sse hyper.cpp; ./a.out
 1.11519e-09
 dkoes_at_quasar:~/tmp$ g++ -mno-sse hyper.cpp; ./a.out
 In file included from /usr/local/include/boost/config/no_tr1/cmath.hpp:21,
                  from
 /usr/local/include/boost/math/policies/error_handling.hpp:15,
                  from
 /usr/local/include/boost/math/distributions/detail/common_error_handling.hpp:12,
                  from
 /usr/local/include/boost/math/distributions/hypergeometric.hpp:12,
                  from hyper.cpp:1:
 /usr/include/c++/4.4/cmath: In function ‘double std::abs(double)’:
 /usr/include/c++/4.4/cmath:94: error: SSE register return with SSE
 disabled
 1.11519e-09
 dkoes_at_quasar:~/tmp$ g++ -mno-sse2 hyper.cpp; ./a.out
 1.11519e-09
 dkoes_at_quasar:~/tmp$ g++ -mno-sse3 hyper.cpp; ./a.out
 1.11519e-09
 dkoes_at_quasar:~/tmp$ g++ -mno-sse4 hyper.cpp; ./a.out
 1.11519e-09
 dkoes_at_quasar:~/tmp$ g++ -DBOOST_MATH_SPECIAL_FUNCTIONS_LANCZOS_SSE2
 hyper.cpp; ./a.out
 1.11519e-09

 }}}

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