Mailman 3 March 2009 - NumPy-Discussion

Intel MKL on Core2 system
by Ryan May March 12, 2009

March 12, 2009

Hi, I noticed the following in numpy/distutils/system_info.py while trying to get numpy to build against MKL: if cpu.is_Itanium(): plt = '64' #l = 'mkl_ipf' elif cpu.is_Xeon(): plt = 'em64t' #l = 'mkl_em64t' else: plt = '32' #l = 'mkl_ia32' So in the autodetection for MKL, the only way to get plt (platform) set to 'em64t' is to test true for a Xeon. This function returns false on my Core2 Duo system, even though the platform is very much 'em64t'. I think that check should instead read: elif cpu.is_Xeon() or cpu.is_Core2(): Thoughts? Ryan -- Ryan May Graduate Research Assistant School of Meteorology University of Oklahoma Sent from: Norman Oklahoma United States.

4 12

1.3 release: getting rid of sourceforge ?
by David Cournapeau March 12, 2009

March 12, 2009

Hi, I was wondering if there was any reason for still using sourceforge ? AFAIK, we only use it to put the files there, and dealing with sourceforge to upload files is less than optimal to say the least. Is there any drawback to directly put the files to scipy.org ? cheers, David

1 0

image processing using numpy-scipy?
by Prashant Saxena March 12, 2009

March 12, 2009

Hi, This a little wiered problem. I am having a black and white image. (black background) Entire image is filled with noisy white patterns of different size and shape. I need to fill the white patches if there area is more then given one. Logically this could possible to use a quickfill algorithm for every pixel and if the filled area generated by that pixel is more then given area(in pixel) then paint that patch (filled area) with black. I have read the docs of PIL but there is no function for this. Can I use numpy-scipy for the matter? The image size is 1K. Any help would be greatly apperciated. Regards Prashant Connect with friends all over the world. Get Yahoo! India Messenger at http://in.messenger.yahoo.com/?wm=n/

4 5

What is the logical value of nan?
by Charles R Harris March 11, 2009

March 11, 2009

It isn't 0 so it should be True. Any disagreement?... Chuck

12 20

Don't call e Euler's constant.
by Charles R Harris March 11, 2009

March 11, 2009

Traditionally, Euler's constant is 0.57721 56649 01532 86060 65120 90082 40243 10421 59335 93992... see wikipedia<http://en.wikipedia.org/wiki/Euler%E2%80%93Mascheroni_constant>. The constant e is sometimes called Euler's number -- shouldn't that be Napier or Bernoulli in a pc world -- but I think e is more universally understood and the distinction between "constant" and "number" is rather obscure. Chuck

3 2

array 2 string
by Mark Bakker March 11, 2009

March 11, 2009

Hello, I want to convert an array to a string. I like array2string, but it puts these annoying square brackets around the array, like [[1 2 3], [3 4 5]] Anyway we can suppress the square brackets and get (this is what is written with savetxt, but I cannot get it to store in a variable) 1 2 3 4 5 6 Thanks, Mark

6 5

Numpy documentation: status and distribution for 1.3.0
by David Cournapeau March 11, 2009

March 11, 2009

Hi, For the upcoming 1.3.0 release, I would like to distribute the (built) documentation in some way. But first, I need to be able to build it :) What are the exact requirements to build the documentation ? Is sphinx 0.5 enough ? I can't manage to build it on either mac os x or linux: ... dumping search index... Exception occurred: File "/Users/david/local/stow/sphinx.dev/lib/python2.5/site-packages/sphinx/search.py", line 151, in get_descrefs pdict[name] = (fn2index[doc], i) KeyError: 'reference/c-api.types-and-structures' The full traceback has been saved in /var/folders/b-/b-BC2bPYFouYhoybrvprFE+++TI/-Tmp-/sphinx-err-PKglvL.log, if you want to report the issue to the author. Please also report this if it was a user error, so that a better error message can be provided next time. There are also some errors on mac os x about too many opened files (which can be alleviated by running the make html again, but obviously, that's not great). I don't know if there are easy solutions to that problem, cheers, David

4 6

Automatic differentiation (was Re: second-order gradient)
by Rob Clewley March 11, 2009

March 11, 2009

> Maybe we should focus on writing a decent 'deriv' function then. I > know Konrad Hinsen's Scientific had a derivatives package > (Scientific.Functions.Derivatives) that implemented automatic > differentiation: > > http://en.wikipedia.org/wiki/Automatic_differentiation That would be great, but wouldn't that be best suited as a utility requiring Sympy? You'll want to take advantage of all sorts of symbolic classes, especially for any source code transformation approach. IMO Hinsen's implementation isn't a very efficient or attractive solution to AD given the great existing C/C++ codes out there. Maybe we should be looking to provide a python interface to an existing open source package such as ADOL-C, but I'm all in favour of a new pure python approach too. What would be perfect is to have a single interface to a python AD package that would support a faster implementation if the user wished to install a C/C++ package, otherwise would default to a pure python equivalent. -Rob

4 4

new numpy error in 1.3.0.dev6618
by Christopher Hanley March 11, 2009

March 11, 2009

====================================================================== ERROR: test_float_repr (test_scalarmath.TestRepr) ---------------------------------------------------------------------- Traceback (most recent call last): File "/Users/chanley/dev/site-packages/lib/python/numpy/core/tests/test_scalarmath.py", line 101, in test_float_repr val2 = t(eval(val_repr)) File "<string>", line 1, in <module> NameError: name 'nan' is not defined ---------------------------------------------------------------------- Ran 2018 tests in 10.311s FAILED (KNOWNFAIL=1, SKIP=1, errors=1) <nose.result.TextTestResult run=2018 errors=1 failures=0> >>> numpy.__version__ '1.3.0.dev6618' >>> This was run on a Intel Mac running OS X 10.5.6. Chris -- Christopher Hanley Senior Systems Software Engineer Space Telescope Science Institute 3700 San Martin Drive Baltimore MD, 21218 (410) 338-4338

3 3

Buildbot issues
by Pauli Virtanen March 10, 2009

March 10, 2009

Hi, There seem to be some problems with the buildbot: - It's not building on new commits automatically. IIRC this could be fixed by removing some (all?) of buildmaster's cache files, and/or switching to PersistentSVNPoller. - The buildmaster apparently has the old 'scipy.org' SVN url that doesn't work any more in its config. - The FreeBSD_64 slave (+ maybe others) has a wrong Numpy test stanza in its Makefile, so tests won't run. -- Pauli Virtanen

3 2