Mailman 3 October 2005 - NumPy-Discussion

Problem with Windows version of Numeric3.0
by Colin J. Williams 08 Oct '05

08 Oct '05

Either of these import lines crashes in Pythonw, with Boa-constructor, SPE or PythonWin. # subclass.py import scipy.base as num3 import scipy.base.multiarray as M Colin W.

1 0

Re: Meaghan exhibitionism
by Fa Mantz 07 Oct '05

07 Oct '05

Hi Do you want t VE UP n your Med o SA TO 70% o dications? It's easy with us - The WebSite VCXAVL iagialanabiealiitr ra $is $xnum $a 134 (30 p.)169 (30 p.) 218 (180 p.) And much more , Good luck

1 0

Numeric3.0 - unable to retrieve data typecode
by Colin J. Williams 07 Oct '05

07 Oct '05

# temp.py import scipy.base as num3 a= num3.array([[1, 2], [9, 8]]) print a.typecode C:\>python temp.py Traceback (most recent call last): File "temp.py", line 4, in ? print a.typecode AttributeError: 'scipy.ndarray' object has no attribute 'typecode' Is there any simple way to retrieve the dataType, as with numarray's numerictypes? Colin W.

1 0

Re: [Numpy-discussion] Re: Purchasing Documentation
by Francesc Altet 07 Oct '05

07 Oct '05

A Divendres 07 Octubre 2005 18:40, Colin J. Williams va escriure: > I share these feelings and would like to express my thanks to Todd > Miller and to Perry Greenfield for their help over the years. Of course, I did not explicitely mention Todd in my previous message, but I would like to note that he has always been extremely responsive to suggestions and very efficient in solving many, many issues that appeared in numarray through the years. A big thanks to you as well, Todd :-) -- >0,0< Francesc Altet http://www.carabos.com/ V V Cárabos Coop. V. Enjoy Data "-"

2 1

wrapping question
by usr01873＠prolocation.net 07 Oct '05

07 Oct '05

Hi all, I'm using the numarray package. I would like to use the numarray C API to create a new array object and fill it with data without copying. Also, I would like numarray to take care of freeing the data once the numarray object is destructed. I've managed to find some related posts in the mailing list archive, but it is really unclear which of these posts are still accurate and what the current 'recommended' approach is. See: http://sourceforge.net/mailarchive/message.php?msg_id=11788304 http://sourceforge.net/mailarchive/message.php?msg_id=2494535 Currently, I first create a new PyArrayObject that is large enough to hold the data: array = NA_vNewArray(NULL,__numarray_type,tmp_num_dims,tmp_dims); then I pass array->data to a C function that reads the data from disk. Any comments on this approach? Also, I was wondering if numarray frees the data using this approach (because it also allocates the memory in the first place within the NA_vNewArray function.) Thanks, Joris van Zwieten

2 3

A comparative test
by Colin J. Williams 06 Oct '05

06 Oct '05

# NumericNumarrayTest.py ''' A test, provided by Francesc Altet http://sourceforge.net/mailarchive/forum.php?thread_id=6396059&forum_id=4890 ''' from timeit import Timer setup = "import Numeric; a = Numeric.arange(2000);a.shape=(1000,2)" print 'Timer-Numeric23.8:', Timer("for i in range(len(a)): row=a[i]", setup).timeit(number=100) # setup = "import numarray; a = numarray.arange(2000);a.shape=(1000,2)" print 'Timer-numarray1.3.3:', Timer("for i in range(len(a)): row=a[i]", setup).timeit(number=100) # scipy/Numeric3 added setup = "import scipy.base; a = scipy.base.arange(2000);a.shape=(1000,2)" print 'Timer-Numeric3:', Timer("for i in range(len(a)): row=a[i]", setup).timeit(number=100) RESULTS: >pythonw -u "NumericNumarrayTest.py" Timer-Numeric23.8: 0.179712784725 Timer-numarray1.3.3: 0.21674933546 Timer-Numeric3: 0.253077136899 >Exit code: 0 Colin W.

6 5

Re: [SciPy-user] Fwd: Scipy_core and numarray
by Travis Oliphant 06 Oct '05

06 Oct '05

Andrew Jaffe wrote: >Hi All, > >I've tried to post via the newsgroup interface, but my messages have >never appeared, so here's another try!... > >I'm currently a numarray user from the astrophysics community and I >suppose I want to understand a bit more about the relationship between >the two projects. Do the two projects exist side-by-side? Will numarray >be superseded by scipy_core? I like the speed advantage of scipy_core, >but I there may be features of numarray I don't want to give up (or >legacy code I need to interfact with). > > I'm interested in the "features" you don't want to give up. SciPy core should have all features of numarray. Also, scipy core, numarray, and Numeric can all live together. If you are using sufficient versions of each, you can exchange data via the array interface. The ONLY reason I started scipy core was to unify Numeric and numarray users. My selling documentation, should not be construed as anything else but trying to help people who can't find time to contribute in other ways to contribute to that goal. If I was "just trying to make money," I would be doing something far different than this... > >Along the latter lines, I have a more practical question: one >numarray capability that I used a lot was the ability to create >uninitialized arrays, such as > > arr = numarray.array(shape=(3,5), type=Float64) > > It's actually there in Numeric (since 23.7 I think). Numeric.empty((3,5),'d') and also in (new) scipy scipy.empty((3,5),scipy.float64) In (new) scipy you can also create new arrays from buffers very easily. -Travis

2 1

A possible problem with PythonWin and Numeric3
by Colin J. Williams 06 Oct '05

06 Oct '05

PythonWin 2.4.1 (#65, Mar 30 2005, 09:13:57) [MSC v.1310 32 bit (Intel)] on win32. Portions Copyright 1994-2004 Mark Hammond (mhammond(a)skippinet.com.au) - see 'Help/About PythonWin' for further copyright information. >>> import scipy <== Executing this causes a crash. >>> import scipy.base <== This has the same effect Colin W.

1 0

pysci base formerly Numeric3
by Colin J. Williams 05 Oct '05

05 Oct '05

I have just downloaded and installed the Windows version. The install appeared to go smoothly but test_basic.py gave some errors, is this what one should expect? --------------------------------------------------------------------------------------------------------------------- >pythonw -u "test_basic.py" !! No test file 'test_scimath.py' found for <module 'scipy.base.scimath' from '...ackages\\scipy\\base\\scimath.pyc'> !! No test file 'test_umath.py' found for <module 'scipy.base.umath' from '...-packages\\scipy\\base\\umath.pyd'> Found 23 tests for scipy.base.function_base !! No test file 'test_machar.py' found for <module 'scipy.base.machar' from '...packages\\scipy\\base\\machar.pyc'> !! No test file 'test_ma.py' found for <module 'scipy.base.ma' from '...ite-packages\\scipy\\base\\ma.pyc'> !! No test file 'test_numerictypes.py' found for <module 'scipy.base.numerictypes' from '...es\\scipy\\base\\numerictypes.pyc'> !! No test file 'test_getlimits.py' found for <module 'scipy.base.getlimits' from '...kages\\scipy\\base\\getlimits.pyc'> Found 9 tests for scipy.base.twodim_base !! No test file 'test__compiled_base.py' found for <module 'scipy.base._compiled_base' from '...\\scipy\\base\\_compiled_base.pyd'> !! No test file 'test_info_scipy_base.py' found for <module 'scipy.base.info_scipy_base' from '...scipy\\base\\info_scipy_base.pyc'> !! No test file 'test_base.py' found for <module 'scipy.base' from '...ckages\\scipy\\base\\__init__.pyc'> !! No test file 'test_convertcode.py' found for <module 'scipy.base.convertcode' from '...ges\\scipy\\base\\convertcode.pyc'> !! No test file 'test_multiarray.py' found for <module 'scipy.base.multiarray' from '...ages\\scipy\\base\\multiarray.pyd'> !! No test file 'test_matrix.py' found for <module 'scipy.base.matrix' from '...packages\\scipy\\base\\matrix.pyc'> !! No test file 'test_oldnumeric.py' found for <module 'scipy.base.oldnumeric' from '...ages\\scipy\\base\\oldnumeric.pyc'> Found 44 tests for scipy.base.shape_base !! No test file 'test_numeric.py' found for <module 'scipy.base.numeric' from '...ackages\\scipy\\base\\numeric.pyc'> !! No test file 'test_polynomial.py' found for <module 'scipy.base.polynomial' from '...ages\\scipy\\base\\polynomial.pyc'> Found 42 tests for scipy.base.type_check !! No test file 'test_arrayprint.py' found for <module 'scipy.base.arrayprint' from '...ages\\scipy\\base\\arrayprint.pyc'> Found 4 tests for scipy.base.index_tricks Found 0 tests for __main__ ................................E.................................................E........E.E.E...........EEEEEE.E..E.... ====================================================================== ERROR: check_simple (scipy.base.shape_base.test_shape_base.test_apply_along_axis) ---------------------------------------------------------------------- Traceback (most recent call last): File "C:\Python24\Lib\site-packages\scipy\base\tests\test_shape_base.py", line 14, in check_simple assert_array_equal(apply_along_axis(len,0,a),len(a)*ones(shape(a)[1])) File "C:\Python24\Lib\site-packages\scipy\base\shape_base.py", line 31, in apply_along_axis res = func1d(arr[tuple(i)],*args) IndexError: each subindex must be either a slice, an integer, Ellipsis, or newaxis ====================================================================== ERROR: check_complex1 (scipy.base.type_check.test_type_check.test_isfinite) ---------------------------------------------------------------------- Traceback (most recent call last): File "C:\Python24\Lib\site-packages\scipy\base\tests\test_type_check.py", line 162, in check_complex1 assert_all(isfinite(array(1+1j)/0.) == 0) ZeroDivisionError: complex division ====================================================================== ERROR: check_neginf_scalar (scipy.base.type_check.test_type_check.test_isinf) ---------------------------------------------------------------------- Traceback (most recent call last): File "C:\Python24\Lib\site-packages\scipy\base\tests\test_type_check.py", line 176, in check_neginf_scalar assert_all(isinf(array(-1.)/0.) == 1) ZeroDivisionError: float division ====================================================================== ERROR: check_posinf_scalar (scipy.base.type_check.test_type_check.test_isinf) ---------------------------------------------------------------------- Traceback (most recent call last): File "C:\Python24\Lib\site-packages\scipy\base\tests\test_type_check.py", line 172, in check_posinf_scalar assert_all(isinf(array(1.,)/0.) == 1) ZeroDivisionError: float division ====================================================================== ERROR: check_complex1 (scipy.base.type_check.test_type_check.test_isnan) ---------------------------------------------------------------------- Traceback (most recent call last): File "C:\Python24\Lib\site-packages\scipy\base\tests\test_type_check.py", line 142, in check_complex1 assert_all(isnan(array(0+0j)/0.) == 1) ZeroDivisionError: complex division ====================================================================== ERROR: check_default_1 (scipy.base.type_check.test_type_check.test_mintypecode) ---------------------------------------------------------------------- Traceback (most recent call last): File "C:\Python24\Lib\site-packages\scipy\base\tests\test_type_check.py", line 18, in check_default_1 assert_equal(mintypecode(itype),'d') File "C:\Python24\Lib\site-packages\scipy\base\type_check.py", line 30, in mintypecode typecodes = [(type(t) is type('') and t) or asarray(t).dtypechar\ UnboundLocalError: local variable 'typecodes' referenced before assignment ====================================================================== ERROR: check_default_2 (scipy.base.type_check.test_type_check.test_mintypecode) ---------------------------------------------------------------------- Traceback (most recent call last): File "C:\Python24\Lib\site-packages\scipy\base\tests\test_type_check.py", line 26, in check_default_2 assert_equal(mintypecode(itype+'f'),'f') File "C:\Python24\Lib\site-packages\scipy\base\type_check.py", line 30, in mintypecode typecodes = [(type(t) is type('') and t) or asarray(t).dtypechar\ UnboundLocalError: local variable 'typecodes' referenced before assignment ====================================================================== ERROR: check_default_3 (scipy.base.type_check.test_type_check.test_mintypecode) ---------------------------------------------------------------------- Traceback (most recent call last): File "C:\Python24\Lib\site-packages\scipy\base\tests\test_type_check.py", line 50, in check_default_3 assert_equal(mintypecode('fdF'),'D') File "C:\Python24\Lib\site-packages\scipy\base\type_check.py", line 30, in mintypecode typecodes = [(type(t) is type('') and t) or asarray(t).dtypechar\ UnboundLocalError: local variable 'typecodes' referenced before assignment ====================================================================== ERROR: check_complex_bad (scipy.base.type_check.test_type_check.test_nan_to_num) ---------------------------------------------------------------------- Traceback (most recent call last): File "C:\Python24\Lib\site-packages\scipy\base\tests\test_type_check.py", line 211, in check_complex_bad v += array(0+1.j)/0. ZeroDivisionError: complex division ====================================================================== ERROR: check_complex_bad2 (scipy.base.type_check.test_type_check.test_nan_to_num) ---------------------------------------------------------------------- Traceback (most recent call last): File "C:\Python24\Lib\site-packages\scipy\base\tests\test_type_check.py", line 217, in check_complex_bad2 v += array(-1+1.j)/0. ZeroDivisionError: complex division ====================================================================== ERROR: check_complex_good (scipy.base.type_check.test_type_check.test_nan_to_num) ---------------------------------------------------------------------- Traceback (most recent call last): File "C:\Python24\Lib\site-packages\scipy\base\tests\test_type_check.py", line 207, in check_complex_good vals = nan_to_num(1+1j) File "C:\Python24\Lib\site-packages\scipy\base\type_check.py", line 99, in nan_to_num y = nan_to_num(x.real) + 1j * nan_to_num(x.imag) File "C:\Python24\Lib\site-packages\scipy\base\type_check.py", line 106, in nan_to_num y[are_nan] = 0 TypeError: object does not support item assignment ====================================================================== ERROR: check_integer (scipy.base.type_check.test_type_check.test_nan_to_num) ---------------------------------------------------------------------- Traceback (most recent call last): File "C:\Python24\Lib\site-packages\scipy\base\tests\test_type_check.py", line 204, in check_integer vals = nan_to_num(1) File "C:\Python24\Lib\site-packages\scipy\base\type_check.py", line 105, in nan_to_num maxf, minf = _getmaxmin(y.dtype) File "C:\Python24\Lib\site-packages\scipy\base\type_check.py", line 86, in _getmaxmin f = getlimits.finfo(t) File "C:\Python24\Lib\site-packages\scipy\base\getlimits.py", line 29, in __init__ raise ValueError, "data type not inexact" ValueError: data type not inexact ====================================================================== ERROR: check_basic (scipy.base.type_check.test_type_check.test_real_if_close) ---------------------------------------------------------------------- Traceback (most recent call last): File "C:\Python24\Lib\site-packages\scipy\base\tests\test_type_check.py", line 230, in check_basic b = real_if_close(a+1e-15j) File "C:\Python24\Lib\site-packages\scipy\base\type_check.py", line 120, in real_if_close tol = f.epsilon * tol AttributeError: 'finfo' object has no attribute 'epsilon' ---------------------------------------------------------------------- Ran 122 tests in 0.521s FAILED (errors=13) >Exit code: -1073741819 ------------------------------------------------------------------------------------------------------------------- Colin W.

1 0

Re: [SciPy-user] Release of SciPy Core 0.4 (Beta)
by Travis Oliphant 04 Oct '05

04 Oct '05

Rob Managan wrote: > Does the new scipy_core support the Numeric function > PyArray_FromDimsAndData? > > I use that to create a front end to a stand alone program that > generates a lot of arrays that I want to be able to query. Not the > best approach maybe but it works! > Yes, All of the old Numeric C-API is available (I believe...). Direct access to descr->one and descr->zero does not work anymore, though (it's replaced with a function call). If you look at the source, you will see that these older calls are all special cases of the call to PyArray_New() -Travis

1 0