[SciPy-dev] bug in scipy.stats.norm.cdf ?

John Hunter jdh2358@gmail....
Sun Mar 4 18:08:15 CST 2007


On 3/4/07, John Hunter <jdh2358@gmail.com> wrote:

> Yep, it's an isnan bug.  If I add a dummy function "erfjdh" to ndtr.c
> and _cephesmodule.c

latest update: the ndtr.c module is using cephes_isnan defined in
cephes/isnan.c.  With some experimenting, it appears I am falling into
the "size(int)==4 and  #ifdef IBMPC" branch of that function, when I
should be in the "size(int)==4 and  #ifdef MIEEE" branch as far as I
can determine.

If I manually define MIEEE in isnan.h, I avoid the unwanted isnans.
Looking at mconf.h where these macros are defined, it looks like I am
not hitting the branch

#elif defined(sel) || defined(pyr) || defined(mc68000) || defined (m68k) || \
          defined(is68k) || defined(tahoe) || defined(ibm032) || \
          defined(ibm370) || defined(MIPSEB) || defined(_MIPSEB) || \
          defined(__convex__) || defined(DGUX) || defined(hppa) || \
          defined(apollo) || defined(_CRAY) || defined(__hppa) || \
          defined(__hp9000) || defined(__hp9000s300) || \
          defined(__hp9000s700) || defined(__AIX) || defined(_AIX) ||\
          defined(__pyr__) || defined(__mc68000__) || defined(__sparc) ||\
          defined(_IBMR2) || defined (BIT_ZERO_ON_LEFT)
#define MIEEE 1     /* Motorola IEEE, high order words come first */
#define BIGENDIAN 1

Is there some additional "defined" check  that needs to be added here?
 Are other mac power pc users seeing this problem?

Thanks!
JDH


More information about the Scipy-dev mailing list