[SciPy-user] scipy.f2py issues

Travis Oliphant oliphant at ee.byu.edu
Fri Nov 4 12:39:34 CST 2005


Chris Fonnesbeck wrote:

>Again, not sure if this belongs on this list or the f2py list, but here goes.
>
>Delving deeper into the f2py array handling issue from my previous
>post, I get a very strange result using a Numeric array with an f2py
>module. The function _poisson() is simply an f2py function that
>returns the log-likelihood of a poisson distribution for particular
>passed values. Here is what happens:
>
>(Pdb) _poisson([3.],[4.])
>Out[2]: -1.6328763961791992
>(Pdb) _poisson(array([3.]),[4.])
>Out[2]: -19775746048.0
>
>Why should a list and an array return different values? Ten orders of
>magnitude different!
>  
>
Good question.  They should both be converted to the same thing 
internally (an ndarray).  Apparently they are not.   We'll look into this. 

-Travis



More information about the SciPy-user mailing list