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!
Thanks,
--
Chris Fonnesbeck
Atlanta, GA