On Tue, Jun 30, 2009 at 12:31 PM, Russell E. Owen
<rowen@u.washington.edu> wrote:
In article
<e06186140906291710s34865590p38032012f12d0f60@mail.gmail.com>,
> > > How big is the extension and what does it do?
> >
> > It basically contains 2 functions:
> > 1: radProfile: given a masked image (2d array), a radius and a desired
> > center: compute a new 1d array whose value at index r is the sum of all
> > unmasked pixels at radius r.
> >
> > 2: radAsymm: given the same inputs as radProfile, return a (scalar)
> > measure of radial asymmetry by computing the variance of unmasked pixels
> > at each radius and combining the results.
> >
> > The original source file is about 1000 lines long, of which 1/3 to 1/2
> > is the basic C code and the rest is Python wrapper.
>
> It sounds small enough that you should be able to update it to the numpy
> interface. What functions do you need? You should also be able to attach a
> copy (zipped) if it is small enough, which might help us help you.
It is the PyGuide package
<http://www.astro.washington.edu/rowen/PyGuide/files/PyGuide.zip>
a 525k zip file. The extension code is in the src directory.
I would certainly be grateful for any pointers to how the old numarray C
API functions map to the new numpy ones. I would prefer to use the new
numpy API if I can figure out what to do.
This doesn't look too bad, I only count 5 functions/macros.
NA_InputArray
NA_OutputArray
NA_ShapeEqual
NA_NewArray
NA_OFFSETDATA
The quick and dirty solution would be to just copy those functions in at the top of your code. You might want to fix up the NumarrayType enum instead of including it and a few other such.
The code looks like it would go over into cython fairly nicely since it is split between interface code, which would look good in python, and a couple of pure c functions. If you have the time that might be a good way to go.
Chuck