Re: [Numpy-discussion] Improvement of performance

2010-05-05 Thread Guilherme P. de Freitas
On Tue, May 4, 2010 at 9:23 PM, josef.p...@gmail.com wrote: In [2] I didn't see anything about higher derivatives, so to get the Hessian I still had to do a finite difference (Jacobian) on the complex_step_grad. Even then the results look pretty good. Yes, the traditional complex step does

Re: [Numpy-discussion] Improvement of performance

2010-05-05 Thread Guilherme P. de Freitas
Just to make this thread more useful for someone interested in these topics, this seems to be the book on automatic differentiation (it's one of the references in the autodiff website) Evaluating Derivatives: Principles and Techniques of Algorithmic Differentiation (2nd ed.), by Andreas Griewank

[Numpy-discussion] 3d plane to point cloud fitting using SVD

2010-05-05 Thread Peter Schmidtke
Dear Numpy Users, I want to fit a 3d plane into a 3d point cloud and I saw that one could use svd for this purpose. So as I am very fond of numpy I saw that svd was implementented in the linalg module. Currently I have a numpy array called xyz with n lines (number of points) and 3 columns

[Numpy-discussion] BUG: NumPy exposes the wrong include directory

2010-05-05 Thread Stéfan van der Walt
Hi all, Under Ubuntu, the NumPy headers are dumped under /usr/include/python2.6. This is not the location that should be exposed when another copy of NumPy is installed. The attached patch lets from numpy.distutils.system_info import get_info print get_info['numpy'] return the correct path.

Re: [Numpy-discussion] 3d plane to point cloud fitting using SVD

2010-05-05 Thread Stéfan van der Walt
Hi Peter On 5 May 2010 10:02, Peter Schmidtke pschmid...@mmb.pcb.ub.es wrote: u,s,vh=numpy.linalg.linalg.svd(M) Then in the matlab analog they use the last column of vh to get the a,b,c coefficients for the equation a,b,c=vh[:, -1] in numpy Note that vh is the conjugate transpose of v. You

[Numpy-discussion] embedding/extending numpy:import _array() problem

2010-05-05 Thread tino
Hi guys, I've a C++ code that calls a python script. The python script, in turn, should use C modules to retrieve data (arrays). I'm having some problem with import_array. I've placed the import_array() just after the initialization of the interpreter and after PyInit_module

[Numpy-discussion] Fwd: [SciPy-Dev] numpy 1.4.1 f2py on Snow Leopard.

2010-05-05 Thread Charles سمير Doutriaux
I'm forwarding this to the numpy list as well since it seems to be f2py related therefore belonging here really.C.Begin forwarded message:From: Charles سمير Doutriaux doutria...@llnl.govDate: April 30, 2010 9:10:08 AM PDTTo: SciPy Developers List scipy-...@scipy.orgSubject: Re: [SciPy-Dev] numpy

Re: [Numpy-discussion] PY_ARRAY_UNIQUE_SYMBOL is too far reaching?

2010-05-05 Thread David
On 05/04/2010 04:38 PM, Austin Bingham wrote: I admit I'm having trouble formulating questions to address my problems, so please bear with me. Say I've got a shared library of utilities for working with numpy arrays. It's intended to be used in multiple extension modules and in some places

Re: [Numpy-discussion] Brooken Toolchain (Mac 10.6)

2010-05-05 Thread Patrick Marsh
I apologize for not following up on this thread right away. I was sick for a couple of days and then there were several major severe weather outbreaks that required me to spend more time at work. With this said, I've resumed trying to build a DMG with Python2.5. The export suggestion below