Mailman 3 March 2006 - SciPy-Dev

Changing return types in C API
by David M. Cooke 19 Apr '06

19 Apr '06

One thing I like about numarray's "native" C API is that for functions that return arrays, the return type is PyArrayObject *, as opposed to the Numeric-compatible API, where the return types are PyObject *. Returning PyObject * in this case basically means you're *always* doing a cast like this: PyObject *o; PyArrayObject *a; a = (PyArrayObject *)PyArray_ContiguousFromObject(o, PyArray_DOUBLE, 2, 2); Could we change the API so that PyArrayObject * is used as the return type? On a similar note: the 'data' member of PyArrayObject's is a char *, where it really should be a void *. Being a void pointer would have the advantage of not needing the cast in cases like this: double *adata = (double)(a->data); It would also mean that accidentally dereferencing the pointer would be a compiler error: currently, a->data[0] is valid, but if a->data was a void *, it wouldn't be. Both of these changes wouldn't require recompiling extensions (since we're just relabelling what the pointer types are), and would be source-compatible with old code (since the casts used would now be superfluous). Newer code will also be easier to write. I'll work on the changes in a svn branch if people think this is a good idea. -- |>|\/|< /------------------------------------------------------------------\ |David M. Cooke http://arbutus.physics.mcmaster.ca/dmc/ |cookedm(a)physics.mcmaster.ca

6 7

scipy.sparse + umfpack + system_info
by Robert Cimrman 10 Apr '06

10 Apr '06

Hi all, I have added a basic umfpack detection support to numpy.distutils.system_info() - see example site.cfg below (and replace <UMFPACKv4.4_root>...). The Umfpack wrapper still resides in the sandbox, but its setup.py now uses the detection stuff. Now I would like to move the wrapper somewhere under scipy.sparse and I hit a problem: umpfack.py depends on scipy.sparse, since it uses the CSR/CSC matrix formats, but also scipy.sparse uses umfpack.py for solving linear equations.if present. How to get out of this circular dependency problem? Do you have any suggestions on how to organize the files? r. -- site.cfg: [umfpack] library_dirs = <UMFPACKv4.4_root>/UMFPACK/Lib:<UMFPACKv4.4_root>/AMD/Lib include_dirs = <UMFPACKv4.4_root>/UMFPACK/Include umfpack_libs = umfpack, amd

4 32

[PATCH] d1mach problem
by Neal Becker 04 Apr '06

04 Apr '06

I think these patches will fix the problem with d1mach miscompiling with gcc4.1:

2 2

Deleting a recipe in the cookbook wiki
by Francesc Altet 01 Apr '06

01 Apr '06

Hi, I've ended changing the name of a recipe of mine in the SciPy cookbook wiki, i.e. from "A Pyrex Agnostic Class" to "A Numerical Agnostic Pyrex Class", which I find better. Anyone with privileges can delete the former one? Thanks, -- >0,0< Francesc Altet http://www.carabos.com/ V V Cárabos Coop. V. Enjoy Data "-"

2 1

problem importing numpy.lib.scimath from gplt
by Brett Olivier 31 Mar '06

31 Mar '06

Hi I picked this up using the the gplt module in recent SVN numpy/scipy versions. When trying to import stuff from numpy.lib.scimath (although it isn't specific to gplt): +++++++++++ In [2]: from numpy.lib.scimath import * --------------------------------------------------------------------------- exceptions.AttributeError Traceback (most recent call last) AttributeError: 'module' object has no attribute 'ERR_CALL' +++++++++++ I'm not sure exactly where to look for the problem. Thanks in advance Brett -- Brett G. Olivier Triple-J Group for Molecular Cell Physiology Stellenbosch University

2 1

FFT with DJBFFT failing
by John Byrnes 30 Mar '06

30 Mar '06

Hello all, I've found that DJBFFT causes scipy to fail several tests. These failures only occur when DJBFFT is found by distutils when building. I'm running the current SVN copies of both scipy and numpy. Regards, John Report is below: ====================================================================== FAIL: check_normal (scipy.stats.tests.test_morestats.test_anderson) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/stats/tests/test_morestats.py", line 51, in check_normal assert_array_less(A, crit[-2:]) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 255, in assert_array_less assert cond,\ AssertionError: Arrays are not less-ordered (mismatch 50.0%): Array 1: 0.90160265659743288 Array 2: [ 0.858 1.0209999999999999] ====================================================================== FAIL: check_definition (scipy.fftpack.tests.test_pseudo_diffs.test_diff) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py", line 85, in check_definition assert_array_almost_equal(diff(sin(x)),direct_diff(sin(x))) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 87.5%): Array 1: [ 1.0000000e+00 9.2387953e-01 7.0710678e-01 3.8268343e-01 8.8009815e-17 -3.8268343e-01 -7.0710678e-01 -9.23... Array 2: [ -0.0000000e+00 5.0000000e-01 7.7662794e-17 1.2967736e-16 1.2490009e-16 2.1612893e-16 1.4972165e-16 -9.71... ====================================================================== FAIL: check_definition (scipy.fftpack.tests.test_pseudo_diffs.test_hilbert) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py", line 301, in check_definition assert_array_almost_equal (y,y1) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 87.5%): Array 1: [ 1.0000000e+00 9.2387953e-01 7.0710678e-01 3.8268343e-01 8.2304423e-17 -3.8268343e-01 -7.0710678e-01 -9.23... Array 2: [ -0.0000000e+00 +0.0000000e+00j 5.0000000e-01 -7.3105720e-17j 3.8831397e-17 -2.4238713e-17j 4.3225785e-17 -6.813... ====================================================================== FAIL: check_random_even (scipy.fftpack.tests.test_pseudo_diffs.test_hilbert) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py", line 332, in check_random_even assert_array_almost_equal(direct_hilbert(direct_ihilbert(f)),f) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 100.0%): Array 1: [ 0. -0.j -0.0003392+0.0003636j -0.0010602-0.0002132j -0.0005238+0.0002887j 0.0003091-0.0020875j 0.00076... Array 2: [-0.1106071 -0.497691 0.4783423 -0.337373 0.2458101 0.0210909 -0.4631084 0.192853 0.4792201 -0.3906347 0.22212... ====================================================================== FAIL: check_tilbert_relation (scipy.fftpack.tests.test_pseudo_diffs.test_hilbert) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py", line 311, in check_tilbert_relation assert_array_almost_equal (y,y1) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 62.5%): Array 1: [ 1.0000000e+00 6.5328148e-01 2.1280927e-16 -2.7059805e-01 -1.0626020e-16 2.7059805e-01 1.5393080e-16 -6.53... Array 2: [ -0.0000000e+00 +0.0000000e+00j 2.5000000e-01 -3.5284004e-17j -1.1905260e-18 +4.4746813e-17j 2.5000000e-01 -9.941... ====================================================================== FAIL: check_definition (scipy.fftpack.tests.test_pseudo_diffs.test_ihilbert) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py", line 370, in check_definition assert_array_almost_equal (y,y1) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 87.5%): Array 1: [ -1.0000000e+00 -9.2387953e-01 -7.0710678e-01 -3.8268343e-01 -8.2304423e-17 3.8268343e-01 7.0710678e-01 9.23... Array 2: [ 0.0000000e+00 -0.0000000e+00j -5.0000000e-01 +7.3105720e-17j -3.8831397e-17 +2.4238713e-17j -4.3225785e-17 +6.813... ====================================================================== FAIL: check_itilbert_relation (scipy.fftpack.tests.test_pseudo_diffs.test_ihilbert) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py", line 380, in check_itilbert_relation assert_array_almost_equal (y,y1) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 62.5%): Array 1: [ -1.0000000e+00 -6.5328148e-01 -2.1280927e-16 2.7059805e-01 1.0626020e-16 -2.7059805e-01 -1.5393080e-16 6.53... Array 2: [ 0.0000000e+00 -0.0000000e+00j -2.5000000e-01 +3.5284004e-17j 1.1905260e-18 -4.4746813e-17j -2.5000000e-01 +9.941... ====================================================================== FAIL: check_definition (scipy.fftpack.tests.test_pseudo_diffs.test_itilbert) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py", line 288, in check_definition assert_array_almost_equal (y,y1) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 87.5%): Array 1: [ -9.9667995e-02 -9.2081220e-02 -7.0475915e-02 -3.8141290e-02 -7.1261911e-18 3.8141290e-02 7.0475915e-02 9.20... Array 2: [ -0.0000000e+00 +0.0000000e+00j -4.9833997e-02 +7.2863005e-18j -7.6643594e-18 +4.7841238e-18j -1.2592216e-17 +1.984... ====================================================================== FAIL: check_definition (scipy.fftpack.tests.test_pseudo_diffs.test_shift) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py", line 390, in check_definition assert_array_almost_equal(shift(sin(x),a),direct_shift(sin(x),a)) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 100.0%): Array 1: [ 0.0998334 0.1968802 0.2920308 0.3843691 0.4730057 0.5570869 0.6358031 0.7083962 0.7741671 0.8324823 0.88278... Array 2: [ 1.2321788e-17 4.9916708e-02 2.4648710e-19 8.8465810e-18 5.9412517e-18 2.0994071e-17 1.4776750e-17 1.22... ====================================================================== FAIL: check_definition (scipy.fftpack.tests.test_pseudo_diffs.test_tilbert) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py", line 227, in check_definition assert_array_almost_equal (y,y1) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 87.5%): Array 1: [ 1.0033311e+01 9.2695708e+00 7.0946223e+00 3.8395819e+00 1.1013865e-15 -3.8395819e+00 -7.0946223e+00 -9.26... Array 2: [ -0.0000000e+00 +0.0000000e+00j 5.0166556e+00 -7.3349243e-16j 1.9673887e-16 -1.2280519e-16j 1.4838281e-16 -2.338... ====================================================================== FAIL: check_random_even (scipy.fftpack.tests.test_pseudo_diffs.test_tilbert) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py", line 240, in check_random_even assert_array_almost_equal(direct_tilbert(direct_itilbert(f,h),h),f) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 100.0%): Array 1: [ -0.0000000e+00 +0.0000000e+00j 1.4398428e-03 +1.8373263e-03j -8.7680356e-04 +9.9658995e-04j -7.3017115e-04 +1.287... Array 2: [ 0.5075109 -0.2380892 -0.2038641 -0.3047494 -0.2114436 0.4541402 -0.1193784 -0.0768124 -0.1561373 0.1053 0.03414... ====================================================================== FAIL: check_definition (scipy.fftpack.tests.test_basic.test_fft) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py", line 98, in check_definition assert_array_almost_equal(y,y1) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 100.0%): Array 1: [ 1.+0.j 2.+0.j 3.+0.j 4.+1.j 1.+0.j 2.+0.j 3.+0.j 4.+2.j] Array 2: [ 20. +3.j -0.7071068+0.7071068j -7. +4.j -0.7071068-0.7071068j -4. -3.j 0.707106... ====================================================================== FAIL: check_definition (scipy.fftpack.tests.test_basic.test_fftn) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py", line 431, in check_definition assert_array_almost_equal(fftn(x),direct_dftn(x)) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 100.0%): Array 1: [[[[ 5.9716569e+02 +0.j -4.7701973e+00 -6.248201j 7.2274468e+00 +2.123387j ..., -1.8242164e+01 +4.20633... Array 2: [[[[ 1.4378537e+02+0.j 3.8003180e-01-0.9697448j 6.8443765e+00+3.6280536j ..., -2.9207977e+00+2.7763892j... ====================================================================== FAIL: check_definition (scipy.fftpack.tests.test_basic.test_ifft) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py", line 183, in check_definition assert_array_almost_equal(y,y1) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 100.0%): Array 1: [ 0.125+0.j 0.25 +0.j 0.375+0.j 0.5 +0.125j 0.125+0.j 0.25 +0.j 0.375+0.j 0.5 +0.25j ] Array 2: [ 2.5 +0.375j 0.0883883+0.0883883j -0.125 -0.5j 0.0883883-0.0883883j -0.5 -0.375j -0.0883883-0.0... ====================================================================== FAIL: check_random_complex (scipy.fftpack.tests.test_basic.test_ifft) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py", line 211, in check_random_complex assert_array_almost_equal (ifft(fft(x)),x) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 100.0%): Array 1: [ 0.000438 +1.2422571e-03j 0.0081279 +1.3153027e-03j 0.0058322 +2.5472850e-03j 0.0098841 +1.5452180e-02j 0.007895... Array 2: [ 0.0280313+0.0795045j 0.5201845+0.0841794j 0.3732589+0.1630262j 0.6325855+0.9889395j 0.5052878+0.7382801j 0.12457... ====================================================================== FAIL: check_random_real (scipy.fftpack.tests.test_basic.test_ifft) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py", line 217, in check_random_real assert_array_almost_equal (ifft(fft(x)),x) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 100.0%): Array 1: [ 0.4665868+0.j 0.0158205-0.0104689j -0.0024869-0.0038622j 0.0255389-0.0128942j 0.0246149-0.0121173j 0.00113... Array 2: [ 0.5673542 0.8727196 0.7445887 0.5776637 0.910135 0.2396252 0.4235496 0.0873539 0.9231279 0.0897232 0.15871... ====================================================================== FAIL: check_definition (scipy.fftpack.tests.test_basic.test_ifftn) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py", line 601, in check_definition assert_array_almost_equal(ifftn(x),direct_idftn(x)) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 100.0%): Array 1: [[[[ 4.9883390e-01 +0.0000000e+00j 7.9039717e-03 +1.7999071e-03j -2.5307390e-03 -2.3994048e-03j ..., 5.8921717... Array 2: [[[[ 1.2726888e-01 +0.0000000e+00j 8.6321110e-04 +2.0907331e-03j 2.2428627e-03 +4.2203592e-04j ..., -2.7428294... ====================================================================== FAIL: check_definition (scipy.fftpack.tests.test_basic.test_irfft) ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py", line 341, in check_definition assert_array_almost_equal(y,ifft(x1)) File "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py", line 231, in assert_array_almost_equal assert cond,\ AssertionError: Arrays are not almost equal (mismatch 100.0%): Array 1: [ 2.625 -1.6856602 -0.375 -1.1856602 0.625 0.4356602 -0.375 0.9356602] Array 2: [ 0.125+0.j 0.25 +0.375j 0.5 +0.125j 0.25 +0.375j 0.5 +0.j 0.25 -0.375j 0.5 -0.125j 0.25 -0.375j] ---------------------------------------------------------------------- Ran 1508 tests in 2.893s FAILED (failures=18)

1 1

1/0
by "Johannes Löhnert" 30 Mar '06

30 Mar '06

Hi, In [48]: int32(1)/int32(0) Out[48]: 0 Why is the result not inf? (I looked up IEEE 854 but do not understand a word...) Johannes -- Echte DSL-Flatrate dauerhaft f�r 0,- Euro*! "Feel free" mit GMX DSL! http://www.gmx.net/de/go/dsl

3 2

masked_array drops fill_value
by Tim Leslie 25 Mar '06

25 Mar '06

Hi all, In ma.py the masked_array function the fill_value of the array passed in is not used. The code currently looks like this: def masked_array (a, mask=nomask, fill_value=None): """masked_array(a, mask=nomask) = array(a, mask=mask, copy=0, fill_value=fill_value) """ return array(a, mask=mask, copy=0, fill_value=fill_value) It seems to me that using the fill_value from a (if it is a MaskedArray) would be the sane thing to do? Something like if fill_value == None and isinstance(a, MaskedArray): fill_value = a.fill_value() return array(a, mask=mask, copy=0, fill_value=fill_value) As it stands, all the ma functions such as transpose, reshape, etc lose the fill_value which seems wrong to me. Tim

1 0

fitpack memory corruption
by Paul Janzen 25 Mar '06

25 Mar '06

The following fragment crashes the Python interpreter on both Windows and Linux due to memory corruption: from scipy.interpolate import splrep splrep(arange(10),arange(10),k=3,task=-1,t=[0,0,0,0,1,1,1,1]) Here is a simple patch against 1712 that at least fixes this case: --- Lib\interpolate\__fitpack.h~ 2006-03-17 04:43:36.000000000 -0800 +++ Lib\interpolate\__fitpack.h 2006-03-25 00:16:45.810340800 -0800 @@ -374,11 +374,11 @@ ap_wrk = (PyArrayObject *)PyArray_FromDims(1,&n,PyArray_DOUBLE); ap_iwrk = (PyArrayObject *)PyArray_FromDims(1,&n,PyArray_INT); if (ap_wrk == NULL || ap_iwrk == NULL) goto fail; + memcpy(ap_wrk->data,wrk,n*sizeof(double)); + memcpy(ap_iwrk->data,iwrk,n*sizeof(int)); } memcpy(ap_t->data,t,n*sizeof(double)); memcpy(ap_c->data,c,lc*sizeof(double)); - memcpy(ap_wrk->data,wrk,n*sizeof(double)); - memcpy(ap_iwrk->data,iwrk,n*sizeof(int)); if (wa) free(wa); Py_DECREF(ap_x); Py_DECREF(ap_y); Similar scenario at line 215:fitpack_surfit. I still don't understand why the condition for copying into ap_{i,}wrk->data (line 373) is iopt==0. Doesn't the curfit documentation imply that you want to persist wrk/iwrk iff iopt==1? Also, it looks like the two assignments at 355 and 370 risk leaking references? ap_t=(PyArrayObject*)PyArray_ContiguousFromObject(t_py,PyArray_DOUBLE, 0, 1); ... ap_t = (PyArrayObject *)PyArray_FromDims(1,&n,PyArray_DOUBLE); As do the second assignments to ap_wrk/ap_iwrk. -- Paul

1 0

PyFITS 1.1 "alpha" release -- NUMPY now supported
by Christopher Hanley 25 Mar '06

25 Mar '06

------------------ | PyFITS Release | ------------------ Space Telescope Science Institute is pleased to announce the "alpha" release of PyFITS 1.1. This release includes support for both the NUMPY and NUMARRAY array packages. This software can be downloaded at: http://www.stsci.edu/resources/software_hardware/pyfits/Download The NUMPY support in PyFITS is not nearly as well tested as the NUMARRAY support. We expect that you will encounter bugs. Please send bug reports to "help(a)stsci.edu". We intend to support NUMARRAY and NUMPY simultaneously for a transition period of no less than 1 year. Eventually, however, support for NUMARRAY will disappear. During this period, it is likely that new features will appear only for NUMPY. The support for NUMARRAY will primarily be to fix serious bugs and handle platform updates. ----------- | Version | ----------- Version 1.1a1; March 24, 2006 ------------------------------- | Major Changes since v1.0.1 | ------------------------------- * Added support for the NUMPY array package * Use of the NUMPY or NUMARRAY array package is controlled through the use of an environmental variable (NUMERIX) * Access to private methods and attributes is now through the pyfits.core module * Now installed as a packaged instead of a singe file ------------------------- | Software Requirements | ------------------------- PyFITS Version 1.0 REQUIRES: * Python 2.3 or later * NUMPY or NUMARRAY --------------------- | Installing PyFITS | --------------------- PyFITS 1.0 is distributed as a Python distutils module. Installation simply involves unpacking the package and executing % python setup.py install to install it in Python's site-packages directory. Alternatively the command %python setup.py install --local="/destination/directory/" will install PyFITS in an arbitrary directory which should be placed on PYTHONPATH. Once numarray or numpy has been installed, then PyFITS should be available for use under Python. ----------------- | Download Site | ----------------- http://www.stsci.edu/resources/software_hardware/pyfits/Download ---------- | Usage | ---------- Users will issue an "import pyfits" command as in the past. However, the use of the NUMPY or NUMARRAY version of PyFITS will be controlled by an environment variable called NUMERIX. Set NUMERIX to 'numarray' for the NUMARRAY version of PyFITS. Set NUMERIX to 'numpy' for the NUMPY version of pyfits. If only one array package is installed, that package's version of PyFITS will be imported. If both packages are installed the NUMERIX value is used to decide which version to import. If no NUMERIX value is set then the NUMARRAY version of PyFITS will be imported. Anything else will raise an exception upon import. --------------- | Bug Reports | --------------- Please send all PyFITS bug reports to help(a)stsci.edu ------------------ | Advanced Users | ------------------ Users who would like to follow the "bleeding" edge of PyFITS development can retrieve the software from our SUBVERSION repository hosted at: http://astropy.scipy.org/svn/pyfits/trunk We also provide a Trac site at: http://projects.scipy.org/astropy/pyfits/wiki -- Christopher Hanley Systems Software Engineer Space Telescope Science Institute 3700 San Martin Drive Baltimore MD, 21218 (410) 338-4338

1 0