Mailman 3 August 2009 - NumPy-Discussion

reduce function of vectorize doesn't respect dtype?
by T J 12 Aug '09

12 Aug '09

The reduce function of ufunc of a vectorized function doesn't seem to respect the dtype. >>> def a(x,y): return x+y >>> b = vectorize(a) >>> c = array([1,2]) >>> b(c, c) # use once to populate b.ufunc >>> d = b.ufunc.reduce(c) >>> c.dtype, type(d) dtype('int32'), <type 'int'> >>> c = array([[1,2,3],[4,5,6]]) >>> b.ufunc.reduce(c) array([5, 7, 9], dtype=object) My goal is to use the output of vectorize() as if it is actually a ufunc. So I'd really like to just type: b.reduce, b.accumulate, etc. And I don't want to have to write (after using it once to populate b.ufunc): b.ufunc.reduce(c).astype(numpy.int32)

2 2

Faulty behavior of numpy.histogram?
by Danny Handoko 12 Aug '09

12 Aug '09

Dear all, We try to use numpy.histogram with combination of matplotlib. We are using numpy 1.3.0, but a somewhat older matplotlib version of 0.91.2. Matplotlib's axes.hist() function calls the numpy.histogram, passing through the 'normed' parameter. However, this version of matplotlib uses '0' as the default value of 'normed' (I see it fixed in higher version). What I found strange is that if the 'normed' parameter of numpy.histogram is set with other object than 'True' or 'False', the output becomes None, but no exceptions are raised. As a result, the matplotlib code that does something like this: >>> n, bins = numpy.histogram([1,2,3], 10, range = None, normed = 0) Traceback (most recent call last): File "<stdin>", line 1, in <module> TypeError: 'NoneType' object is not iterable results in the above exception. Secondly, this matplotlib version also expects both outputs to be of the same length, which is no longer true with the new histogram semantics. This can be easily reverted using the parameter 'new = False' in numpy.histogram, but this parameter is not available for the caller of axes.hist() function in matplotlib. Is there any way to tell numpy to use the old semantics? Upgrading to the newer matplotlib is a rather longer term solution, and we hope to be able to find some workaround/short-term solution Thank you, -- Danny Handoko System Architecture and Generics Room 7G2.003 -- ph: x2968 email: danny.handoko(a)asml.com -- The information contained in this communication and any attachments is confidential and may be privileged, and is for the sole use of the intended recipient(s). Any unauthorized review, use, disclosure or distribution is prohibited. Unless explicitly stated otherwise in the body of this communication or the attachment thereto (if any), the information is provided on an AS-IS basis without any express or implied warranties or liabilities. To the extent you are relying on this information, you are doing so at your own risk. If you are not the intended recipient, please notify the sender immediately by replying to this message and destroy all copies of this message and any attachments. ASML is neither liable for the proper and complete transmission of the information contained in this communication, nor for any delay in its receipt.

2 1

speeding up getting a subset of a data array
by Brennan Williams 10 Aug '09

10 Aug '09

Hi No doubt asked many times before so apologies.... I'm pulling a subset array out of a data array where I have a list of the indices I want (could be an array rather than a list actually - I have it in both). Potentially the number of points and the number of times I do this can get very large so any saving in time is good. So, paraphrasing what I've currently got.... say I have... subsetpointerlist=[0,1,2,5,8,15,25...] subsetsize=len(subsetpointerlist) subsetarray=zeros(subsetsize,dtype=float) for index,pos in enumerate(subsetpointerlist): subsetarray[index]=dataarray[pos] How do I speed this up in numpy, i.e. by removing the for loop? Do I set up some sort of a subsetpointerarray as a mask and then somehow apply that to dataarray to get the values into subsetarray? Thanks Brennan

2 3

vectorize problem with f2py and gfortran 4.3
by David Huard 10 Aug '09

10 Aug '09

Hi all, A user on the pymc user list has reported a problem with f2py wrapped fortran functions compiled with gfortran 4.3, which is the standard Ubuntu Jaunty fortran compiler. I noticed the same bug in some of my own routines. The problem, as far as I can understand, is that vectorize tries to find the number of arguments by calling the function with no arguments and parsing the error message. With numpy 1.3, python 2.6 and gfortran 4.3, the error message is not what numpy expects, and does not contain the expected number of arguments. So I am wondering if there is a reliable way to introspect compiled extensions to provide the number of arguments needed by vectorize ? Thanks, David

1 0

SWIG, numpy.i and errno: comments?
by Egor Zindy 10 Aug '09

10 Aug '09

Hello list, this is my attempt at generating python exceptions in SWIG/C using the errno mechanism: http://www.scipy.org/Cookbook/SWIG_NumPy_examples#head-10f49a0f5ea6b313127d… Used together with numpy.i, this has been useful for notifying (in a pythonic way) memory allocation errors or array index problems. A change in the errno global variable is detected in the %exception part of the SWIG interface file, and Python exceptions are generated after $action depending on the errno error code value. %exception { errno = 0; $action if (errno != 0) { switch(errno) { case EPERM: PyErr_Format(PyExc_IndexError, "Index out of range"); break; case ENOMEM: PyErr_Format(PyExc_MemoryError, "Failed malloc()"); break; default: PyErr_Format(PyExc_Exception, "Unknown exception"); } return NULL; } } If there's a better way of doing this, I'll update the cookbook recipe. Regards, Egor

2 3

problem during installation
by Giuseppe Aprea 10 Aug '09

10 Aug '09

Sorry if I posted this twice. I wonder if anyone can suggest me how to fix this error which I find during numpy build: ......... creating build/temp.linux-i686-2.6/numpy/linalg compile options: '-DNO_ATLAS_INFO=1 -Inumpy/core/include -Ibuild/src.linux-i686-2.6/numpy/core/include/numpy -Inumpy/core/src -Inumpy/core/include -I/home/gaprea/usr/local/include/python2.6 -c' gcc: numpy/linalg/lapack_litemodule.c gcc: numpy/linalg/python_xerbla.c /usr/bin/gfortran -Wall -L/home/gaprea/usr/local/lib build/temp.linux-i686-2.6/numpy/linalg/lapack_litemodule.o build/temp.linux-i686-2.6/numpy/linalg/python_xerbla.o -L/home/gaprea/usr/local/lib -L/usr/lib -L/home/gaprea/usr/local/lib -Lbuild/temp.linux-i686-2.6 -llapack -lblas -lpython2.6 -lgfortran -o build/lib.linux-i686-2.6/numpy/linalg/lapack_lite.so /usr/lib/gcc/i486-linux-gnu/4.2.4/libgfortranbegin.a(fmain.o): In function `main': (.text+0x23): undefined reference to `MAIN__' collect2: ld returned 1 exit status /usr/lib/gcc/i486-linux-gnu/4.2.4/libgfortranbegin.a(fmain.o): In function `main': (.text+0x23): undefined reference to `MAIN__' collect2: ld returned 1 exit status error: Command "/usr/bin/gfortran -Wall -L/home/gaprea/usr/local/lib build/temp.linux-i686-2.6/numpy/linalg/lapack_litemodule.o build/temp.linux-i686-2.6/numpy/linalg/python_xerbla.o -L/home/gaprea/usr/local/lib -L/usr/lib -L/home/gaprea/usr/local/lib -Lbuild/temp.linux-i686-2.6 -llapack -lblas -lpython2.6 -lgfortran -o build/lib.linux-i686-2.6/numpy/linalg/lapack_lite.so" failed with exit status 1 It seems that the installation program is using gfortran to link while it should have used gcc. I am using Kubuntu 8.04 and I have the following version of gfortran and gcc: $ gfortran -v Using built-in specs. Target: i486-linux-gnu Configured with: ../src/configure -v --enable-languages=c,c++,fortran,objc,obj-c++,treelang --prefix=/usr --enable-shared --with-system-zlib --libexecdir=/usr/lib --without-included-gettext --enable-threads=posix --enable-nls --with-gxx-include-dir=/usr/include/c++/4.2 --program-suffix=-4.2 --enable-clocale=gnu --enable-libstdcxx-debug --enable-objc-gc --enable-mpfr --enable-targets=all --enable-checking=release --build=i486-linux-gnu --host=i486-linux-gnu --target=i486-linux-gnu Thread model: posix gcc version 4.2.4 (Ubuntu 4.2.4-1ubuntu4) $ gcc -v Using built-in specs. Target: i486-linux-gnu Configured with: ../src/configure -v --enable-languages=c,c++,fortran,objc,obj-c++,treelang --prefix=/usr --enable-shared --with-system-zlib --libexecdir=/usr/lib --without-included-gettext --enable-threads=posix --enable-nls --with-gxx-include-dir=/usr/include/c++/4.2 --program-suffix=-4.2 --enable-clocale=gnu --enable-libstdcxx-debug --enable-objc-gc --enable-mpfr --enable-targets=all --enable-checking=release --build=i486-linux-gnu --host=i486-linux-gnu --target=i486-linux-gnu Thread model: posix gcc version 4.2.4 (Ubuntu 4.2.4-1ubuntu4) does anybody have a suggestion? thanks giuseppe

1 0

Sanity checklist for those attending the SciPy'09 tutorials
by Fernando Perez 09 Aug '09

09 Aug '09

Hi all, [ sorry for spamming the list, but even though I sent this to all the email addresses I have on file for tutorial attendees, I know I am missing a few, so I hope they see this message. ] In order to make your experience at the scipy tutorials as smooth as possible, we strongly recommend that you take a little time to install the necessary tools in advance. For both introductory and advanced tutorials: http://conference.scipy.org/intro_tutorial http://conference.scipy.org/advanced_tutorials you will find instructions on what to install and where to download it from. In addition (this is also mentioned on those pages), we encourage you to run according to your tutorial of choice, a little checklist script: https://cirl.berkeley.edu/fperez/tmp/intro_tut_checklist.py https://cirl.berkeley.edu/fperez/tmp/adv_tut_checklist.py This will try to spot any problems early, and we'll do our best to help you with them before you arrive to the conference. Best regards, Dave Peterson and Fernando Perez. ps - for those of you who may find fixes for the checklist scripts, the sources are hosted on github: http://github.com/fperez/scipytut/

2 2

Indexing with a list...
by T J 09 Aug '09

09 Aug '09

>>> z = array([1,2,3,4]) >>> z[[1]] array([1]) >>> z[(1,)] 1 I'm just curious: What is the motivation for this differing behavior? Is it a necessary consequence of, for example, the following: >>> z[z<3] array([1,2])

3 4

Test failures r7300
by lukshuntim＠gmail.com 09 Aug '09

09 Aug '09

Hi, I got 16 test failures after building r7300 from svn on debian/sid/i386. Seems all related to complex linear algebra modules. Here's the error messages: <quote> Running unit tests for numpy NumPy version 1.4.0.dev7300 NumPy is installed in /var/opt/py/lib/python2.5/site-packages/numpy Python version 2.5.4 (r254:67916, Feb 17 2009, 20:16:45) [GCC 4.3.3] nose version 0.11.1 ... FAIL: test_cdouble (test_linalg.TestCond2) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 44, in test_cdouble self.do(a, b) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 114, in do old_assert_almost_equal(s[0]/s[-1], linalg.cond(a,2), decimal=5) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 421, in assert_almost_equal raise AssertionError(msg) AssertionError: Arrays are not almost equal ACTUAL: 9.4348091510413177 DESIRED: 22.757141876814547 ====================================================================== FAIL: test_csingle (test_linalg.TestCond2) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 39, in test_csingle self.do(a, b) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 114, in do old_assert_almost_equal(s[0]/s[-1], linalg.cond(a,2), decimal=5) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 421, in assert_almost_equal raise AssertionError(msg) AssertionError: Arrays are not almost equal ACTUAL: 9.4348097 DESIRED: 22.757143 ====================================================================== FAIL: test_cdouble (test_linalg.TestDet) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 44, in test_cdouble self.do(a, b) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 135, in do assert_almost_equal(d, multiply.reduce(ev)) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: (8.881784197e-16-4j) DESIRED: (5.28-11.04j) ====================================================================== FAIL: test_csingle (test_linalg.TestDet) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 39, in test_csingle self.do(a, b) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 135, in do assert_almost_equal(d, multiply.reduce(ev)) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: (8.881784197e-16-4j) DESIRED: (5.28-11.04j) ====================================================================== FAIL: test_cdouble (test_linalg.TestEig) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 44, in test_cdouble self.do(a, b) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 94, in do assert_almost_equal(dot(a, evectors), multiply(evectors, evalues)) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: [[ 2.72530404+2.67511327j 1.91551601+1.28537403j] [ 5.95809316+4.79684551j 3.39598770+1.39546789j]] DESIRED: [[ 2.01388405+1.03693361j -1.39855180+1.88751398j] [ 1.78601662+0.01838201j -0.10378837-3.53101635j]] ====================================================================== FAIL: test_csingle (test_linalg.TestEig) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 39, in test_csingle self.do(a, b) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 94, in do assert_almost_equal(dot(a, evectors), multiply(evectors, evalues)) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: [[ 2.72530413+2.6751132j 1.91551602+1.28537405j] [ 5.95809317+4.79684544j 3.39598775+1.39546788j]] DESIRED: [[ 2.01388407+1.03693354j -1.39855182+1.887514j ] [ 1.78601658+0.01838197j -0.10378837-3.53101635j]] ====================================================================== FAIL: test_cdouble (test_linalg.TestEigh) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 221, in test_cdouble self.do(a) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 259, in do assert_almost_equal(ev, evalues) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: [-2.60555128 4.60555128] DESIRED: [-1.71080202-1.00413682j 3.01849433+1.46567528j] ====================================================================== FAIL: test_csingle (test_linalg.TestEigh) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 217, in test_csingle self.do(a) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 259, in do assert_almost_equal(ev, evalues) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: [-2.60555124 4.60555124] DESIRED: [-1.71080208-1.0041368j 3.01849437+1.46567523j] ====================================================================== FAIL: test_cdouble (test_linalg.TestEigvalsh) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 221, in test_cdouble self.do(a) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 249, in do assert_almost_equal(ev, evalues) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: [-2.60555128+0.j 4.60555128+0.j] DESIRED: [-1.71080202-1.00413682j 3.01849433+1.46567528j] ====================================================================== FAIL: test_csingle (test_linalg.TestEigvalsh) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 217, in test_csingle self.do(a) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 249, in do assert_almost_equal(ev, evalues) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: [-2.60555124+0.j 4.60555124+0.j] DESIRED: [-1.71080208-1.0041368j 3.01849437+1.46567523j] ====================================================================== FAIL: test_cdouble (test_linalg.TestLstsq) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 44, in test_cdouble self.do(a, b) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 141, in do assert_almost_equal(b, dot(a, x)) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: [ 2.+1.j 1.+2.j] DESIRED: [ 0.95920929+0.98311952j 1.23494444+0.67346351j] ====================================================================== FAIL: test_csingle (test_linalg.TestLstsq) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 39, in test_csingle self.do(a, b) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 141, in do assert_almost_equal(b, dot(a, x)) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: [ 2.+1.j 1.+2.j] DESIRED: [ 0.95920926+0.98311943j 1.23494434+0.67346334j] ====================================================================== FAIL: test_cdouble (test_linalg.TestPinv) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 44, in test_cdouble self.do(a, b) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 124, in do assert_almost_equal(dot(a, a_ginv), identity(asarray(a).shape[0])) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: [[ 0.29169056-0.07799046j 0.17767375-0.01332484j] [ 0.04125021-0.38255608j 0.73402869+0.62377356j]] DESIRED: [[ 1. 0.] [ 0. 1.]] ====================================================================== FAIL: test_csingle (test_linalg.TestPinv) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 39, in test_csingle self.do(a, b) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 124, in do assert_almost_equal(dot(a, a_ginv), identity(asarray(a).shape[0])) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: [[ 0.29169053-0.07799049j 0.17767370-0.0133248j ] [ 0.04125014-0.38255614j 0.73402858+0.62377363j]] DESIRED: [[ 1. 0.] [ 0. 1.]] ====================================================================== FAIL: test_cdouble (test_linalg.TestSVD) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 44, in test_cdouble self.do(a, b) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 100, in do assert_almost_equal(a, dot(multiply(u, s), vt)) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: [[ 1.+2.j 2.+3.j] [ 3.+4.j 4.+5.j]] DESIRED: [[ 1.00000000+2.j 2.36670415+2.98574489j] [ 3.00000000+4.j 2.80882652+6.25521741j]] ====================================================================== FAIL: test_csingle (test_linalg.TestSVD) ---------------------------------------------------------------------- Traceback (most recent call last): File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 39, in test_csingle self.do(a, b) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 100, in do assert_almost_equal(a, dot(multiply(u, s), vt)) File "/var/opt/py/lib/python2.5/site-packages/numpy/linalg/tests/test_linalg.py", line 23, in assert_almost_equal old_assert_almost_equal(a, b, decimal=decimal, **kw) File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py", line 400, in assert_almost_equal "DESIRED: %s\n" % (str(actual), str(desired))) AssertionError: Items are not equal: ACTUAL: [[ 1.+2.j 2.+3.j] [ 3.+4.j 4.+5.j]] DESIRED: [[ 0.99999994+2.j 2.36670423+2.98574495j] [ 3.00000000+4.j 2.80882668+6.25521755j]] ---------------------------------------------------------------------- Ran 2186 tests in 10.594s FAILED (KNOWNFAIL=1, SKIP=11, failures=16) </quote> Regards, ST --

2 4

Merging datetime, Yes or No
by Charles R Harris 08 Aug '09

08 Aug '09

I ask again, Datetime is getting really stale and hasn't been touched recently. Do the datetime folks want it merged or not, because it's getting to be a bit of work. Chuck

3 3