Re: [Numpy-discussion] Help - numpy / scipy binary compatibility

2014-08-08 Thread David Cournapeau
On Sat, Aug 9, 2014 at 10:41 AM, Matthew Brett wrote: > Hi, > > I would be very happy of some help trying to work out a numpy package > binary incompatibility. > > I'm trying to work out what's happening for this ticket: > > https://github.com/scipy/scipy/issues/3863 > > which I summarized at the

[Numpy-discussion] Help - numpy / scipy binary compatibility

2014-08-08 Thread Matthew Brett
Hi, I would be very happy of some help trying to work out a numpy package binary incompatibility. I'm trying to work out what's happening for this ticket: https://github.com/scipy/scipy/issues/3863 which I summarized at the end: https://github.com/scipy/scipy/issues/3863#issuecomment-51669861

Re: [Numpy-discussion] OpenBLAS and dotblas

2014-08-08 Thread Sturla Molden
Charles R Harris wrote: > It looks like numpy dot only uses BLAS if ATLAS is present, see > numpy/core/setup.py. Has anyone done the mods needed to use OpenBLAS? What > is the current status of using OpenBLAS with numpy? I thought it also uses BLAS if MKL or Accerate Framework is present, but I

Re: [Numpy-discussion] Calculation of a hessian

2014-08-08 Thread Eelco Hoogendoorn
Do it in pure numpy? How about copying the source of numdifftools? What exactly is the obstacle to using numdifftools? There seem to be no licensing issues. In my experience, its a crafty piece of work; and calculating a hessian correctly, accounting for all kinds of nasty floating point issues, i

Re: [Numpy-discussion] Calculation of a hessian

2014-08-08 Thread Jose Gomez-Dans
Your function looks fairly simple to differentiate by hand, but if you have access to the gradient (or you estimate it numerically using scipy...), this function might do the job: def hessian ( x, the_func, epsilon=1e-8): """Numerical approximation to the Hessian Parameters ---

[Numpy-discussion] Calculation of a hessian

2014-08-08 Thread Kiko
Hi all, I am trying to calculate a Hessian. I am using numdifftools for this ( https://pypi.python.org/pypi/Numdifftools). My question is, is it possible to make it using pure numpy?. The actual code is like this: *import numdifftools as nd* *import numpy as np* *def log_likelihood(params):*