Mailman 3 June 2006 - NumPy-Discussion

How do I use numpy to do this?
by Christopher Barker 03 Jun '06

03 Jun '06

I'm trying to get the (x,y) coords for all the points in a grid, bound by xmin, xmax, ymin, ymax. This list comprehension does it fine: Points = [(x,y) for x in xrange(minx, maxx) for y in xrange(miny, maxy)] But I can't think at the moment how to do it with numpy. Any ideas? Thanks, -Chris -- Christopher Barker, Ph.D. Oceanographer NOAA/OR&R/HAZMAT (206) 526-6959 voice 7600 Sand Point Way NE (206) 526-6329 fax Seattle, WA 98115 (206) 526-6317 main reception Chris.Barker(a)noaa.gov

5 13

numpy vs numeric benchmarks
by Eric Jonas 02 Jun '06

02 Jun '06

Hello! I've been using numeric for a while, and the recent list traffic prompted me to finally migrate all my old code. On a whim, we were benchmarking numpy vs numeric and have been lead to the conclusion that numpy is at least 50x slower; a 1000x1000 matmul takes 16 sec in numpy but 300 ms in numeric. Now, of course, I don't believe this, but I can't figure out what we're doing wrong; I'm not the only person who has looked at this code, so can anyone tell me what we're doing wrong? We run both benchmarks twice to try and mitigate any start-up and cache effects. This is with debian-amd64's packaged numeric 24.2-2 and a locally built numpy-0.9.8. /usr/bin/python import time import numpy import random import Numeric def numpytest(): N = 1000 x = numpy.zeros((N,N),'f') y = numpy.zeros((N,N),'f') for i in range(N): for j in range(N): x[i, j] = random.random() y[i, j] = random.random() t1 = time.clock() z = numpy.matrixmultiply(x, y) t2 = time.clock() print (((t2 - t1)*1000)) def numerictest(): N = 1000 x = Numeric.zeros((N,N),'f') y = Numeric.zeros((N,N),'f') for i in range(N): for j in range(N): x[i, j] = random.random() y[i, j] = random.random() t1 = time.clock() z = Numeric.matrixmultiply(x, y) t2 = time.clock() print (((t2 - t1)*1000)) numerictest() numpytest() numpytest() numerictest() on our hardware a call to numerictest() takes 340 ms and a numpytest takes around 13 sec (!). Any advice on what we're doing wrong would be very helpful. ...Eric

8 13

Re: [Numpy-discussion] Free SciPy 2006 porting service
by r.demaria＠tiscali.it 02 Jun '06

02 Jun '06

Hi, maybe is not what you meant, but presently I'm looking for a sparse eigenvalue solver. As far as I've understood the ARPACK bindings are still missing. This library is one of the most used, so I think it would be very useful to have integrated in numpy. Riccardo La gara più entusiasmante dell'anno! Gioca e corri alla velocità della luce sui 18 circuiti di Intel Speed Contest 2006! I più bravi vincono Notebook Sony VAIO, iPod da 60 GB e altro ancora... Sfida gli amici! http://intelspeedcontest2006.tiscali.it/

2 1

Re: [Numpy-discussion] numpy vs numeric benchmarks
by RayS 02 Jun '06

02 Jun '06

favorable numpy creates arrays much faster, fft seems a tad faster a useful metric, I think, for O-scope and ADC apps I get 0.0039054614015815738 0.0019759541205486885 0.023268623246481726 0.0023570392204637913 from the below on a PIII 600... from time import * n=4096 r = range(n) #numpy import numpy arr = numpy.array # array creation t0 = clock() for i in r: a = arr(r) (clock()-t0)/float(n) #fft of n fftn = numpy.fft t0 = clock() for i in r: f = fftn(a) (clock()-t0)/float(n) #Numeric import Numeric arr = Numeric.array # array creation t0 = clock() for i in r: a = arr(r) (clock()-t0)/float(n) #fft of n from FFT import * t0 = clock() for i in r: f = fft(a) (clock()-t0)/float(n)

1 0

lexsort
by Tom Denniston 02 Jun '06

02 Jun '06

This function is really useful but it seems to only take tuples not ndarrays. This seems kinda strange. Does one have to convert the ndarray into a tuple to use it? This seems extremely inefficient. Is there an efficient way to argsort a 2d array based upon multiple columns if lexsort is not the correct way to do this? The only way I have found to do this is to construct a list of tuples and sort them using python's list sort. This is inefficient and convoluted so I was hoping lexsort would provide a simple solution. --Tom

3 3

Free SciPy 2006 porting service
by Travis Oliphant 01 Jun '06

01 Jun '06

I will be available during the SciPy 2006 conference to help port open-source applications to NumPy for no charge. (I'm always available for porting commercial code for a reasonable fee). Others who want to assist will be welcome. Conference attendees will get first priority, but others who want to email their request can do so. Offer will be on a first come, first serve basis but I will reserve the liberty to rearrange the order to serve as many projects as possible. I'll place a note on the Wiki Coding Sprint page to this effect. -Travis O.

1 0

Re: Any Numeric or numarray users on this list?
by Travis Oliphant 01 Jun '06

01 Jun '06

Berthold Höllmann wrote: > Travis Oliphant <oliphant.travis(a)ieee.org> writes: > > >> 2) Will you transition within the next 6 months? (if you answered No to #1) >> > > Unlikely > > >> 3) Please, explain your reason(s) for not making the switch. (if you >> answered No to #2) >> > > Lack of resources (Numeric is used in hand coded extensions; are > arrays of type PyObject supported in NumPy, they were not in numarray) > Yes, NumPy is actually quite similar to Numeric. Most C-extensions are easily ported simply by replacing #include Numeric/arrayobject.h with #include numpy/arrayobject.h (and making sure you get the right location for the headers). -Travis

2 1

What am I missing about concatenate?
by Christopher Barker 01 Jun '06

01 Jun '06

I want to take two (2,) arrays and put them together into one (2,2) array. I thought one of these would work: >>> N.concatenate(((1,2),(3,4)),0) array([1, 2, 3, 4]) >>> N.concatenate(((1,2),(3,4)),1) array([1, 2, 3, 4]) Is this the best I can do? >>> N.concatenate(((1,2),(3,4))).reshape(2,2) array([[1, 2], [3, 4]]) Is it because the arrays I'm putting together are rank-1? >>> N.__version__ '0.9.6' -Chris -- Christopher Barker, Ph.D. Oceanographer NOAA/OR&R/HAZMAT (206) 526-6959 voice 7600 Sand Point Way NE (206) 526-6329 fax Seattle, WA 98115 (206) 526-6317 main reception Chris.Barker(a)noaa.gov

5 6

{Reply} Nothing like it INFX
by zpxknsxt wwacjenon 01 Jun '06

01 Jun '06

INFX**INFX**INFX**INFX**INFX**INFX**INFX**INFX** Infinex Ventures Inc. (INFX) Current Price: $0.52 The Rally has begun Watch this one like a hawk, this report is sent because the potential is incredible This is AS sure as it gets H U G E N E W S read below COMPANY OVERVIEW Aggressive and energetic, Infinex boasts a dynamic and diversified portfolio of operations across North America, with an eye on international expansion. Grounded in natural resource exploration, Inifinex also offers investors access to exciting new developments in the high-tech sector and the booming international real estate market. Our market based experience, tenacious research techniques, and razor sharp analytical skills allow us to leverage opportunities in emerging markets and developing technologies. Identifying these opportunities in the earliest stages allows us to accelerate business development and fully realize the companys true potential. Maximizing overall profitability and in turn enhancing shareholder value. Current Press Release Infinex Announces Extension to Its Agreement in Chile LAS VEGAS, NV, May 9 /PRNewswire-FirstCall/ - Infinex Ventures Inc. (INFX:OB - News; "the Company") and its Board of Directors are pleased to announce that the Company has received an extension (90 days) to its Agreement for the due diligence period, in an effort to fully verify the offered title and all additional documentation, including but not limited to, Trial C-1912- 2001 at the 14th Civil Court of Santiago and Criminal Trial 1160-2002 at the 19th Court of Crime of Santiago of Chile, Ministry of Mines of Chile over its sole and exclusive right to acquire a 50% interest in the Tesoro 1-12 Mining Claims. Infinex Announces Joint Venture and Option Agreement Extension LAS VEGAS, May 5 /PRNewswire-FirstCall/ - Infinex Ventures Inc. (INFX:OB - "the Company") and its Board of Directors are please to announce that the Company has been granted an extension of 120 days to fulfill its contractual obligations under the Joint Venture and Option Agreement dated June 14, 2004 on the Texada Island "Yew Gr0up" Mining Claims: Shake like a leaf. We'll hand you out to dry. We'll hand you out to dry. Stand your ground. Ugly as a mud fence. Scraping the bottom of the barrel. Sly as a fox. A snail's pace. Your ass is grass. There is always next year. Rise and shine. Sly as a fox. Putting the cart before the horse. Walking on thin ice. Stand your ground. Root it out. The stronger the breeze the stronger the trees. Putting it in a nutshell. This is for the birds. Wrinkled as a prune. Up one side and down the other. Rise and shine. Sour as a green apple. You say potayto, I say potahto. She's a mother hen. Say it with flowers. A thorn in my side. Weed out.

1 0

NumPy-Discussion June 2006