Re: [Numpy-discussion] dot() performance depends on data?

2010-09-11 Thread Hagen Fürstenau
Anyway, seems it is indeed a denormal issue, as adding a small (1e-10) constant gives same speed for both timings. With adding 1e-10 or clipping to 0 at 1e-150, I still get a slowdown of about 30% compared with the random arrays. Any explanation for that? Cheers, Hagen signature.asc

Re: [Numpy-discussion] dot() performance depends on data?

2010-09-11 Thread Matthieu Brucher
Denormal numbers are a tricky beast. You may have to change the clip or the shift depending on the processor you have. It's no wonder that processors and thus compilers have options to round denormals to zero. Matthieu 2010/9/11 Hagen Fürstenau ha...@zhuliguan.net: Anyway, seems it is indeed a

[Numpy-discussion] scan array to extract min-max values (with if condition)

2010-09-11 Thread Massimo Di Stefano
Hello All, i need to extract data from an array, that are inside a rectangle area defined as : N, S, E, W = 234560.94503118, 234482.56929822, 921336.53116178, 921185.3779625 the data are in a csv (comma delimited text file, with 3 columns X,Y,Z) #X,Y,Z 3020081.5500,76.3100,0.0300

Re: [Numpy-discussion] scan array to extract min-max values (with if condition)

2010-09-11 Thread Brett Olsen
On Sat, Sep 11, 2010 at 7:45 AM, Massimo Di Stefano massimodisa...@gmail.com wrote: Hello All, i need to extract data from an array, that are inside a rectangle area defined as : N, S, E, W = 234560.94503118, 234482.56929822, 921336.53116178, 921185.3779625 the data are in a csv (comma

Re: [Numpy-discussion] scan array to extract min-max values (with if condition)

2010-09-11 Thread Massimo Di Stefano
That's awesome! masked array are defintley what i need! thanks to point my attention on it! best regards, Massimo. Il giorno 11/set/2010, alle ore 16.19, Brett Olsen ha scritto: On Sat, Sep 11, 2010 at 7:45 AM, Massimo Di Stefano massimodisa...@gmail.com wrote: Hello All, i need to

Re: [Numpy-discussion] ANNOUNCE: mahotas 0.5

2010-09-11 Thread Luis Pedro Coelho
On Friday, September 10, 2010 03:40:33 am Sebastian Haase wrote: Hi Luis, thanks for the announcement. How would you compare mahotas to scipy's ndimage ? Are you using ndimage in mahotas at all ? Hi Sebastian, In general there is little overlap (there are 1 or 2 functions which are

Re: [Numpy-discussion] scan array to extract min-max values (with if condition)

2010-09-11 Thread Massimo Di Stefano
Brett, i tried a different way to solve the problem, using : # import os fpath = '/Users/sasha/py/' input_fp = open( os.path.join(fpath, 'BE3730072600WC20050817.txt'), 'r' ) input_file = input_fp.readlines() N = 234560.94503118 S = 234482.56929822 E = 921336.53116178 W =

Re: [Numpy-discussion] scan array to extract min-max values (with if condition)

2010-09-11 Thread Pierre GM
On Sep 11, 2010, at 9:53 PM, Massimo Di Stefano wrote: have you hints on how to get working the MASK code ? as it is now it pick all the points in the mydata array. Brett's code for the mask matched the loop of your post. However, taking a second look at it, I don't see why it would work.

Re: [Numpy-discussion] scan array to extract min-max values (with if condition)

2010-09-11 Thread Massimo Di Stefano
Thanks Pierre, i tried it and all works fine and fast. my apologize :-( i used a wrong if statment to represent my needs if mydata[i,0] E or mydata[i,0] W or mydata[i,1] N or mydata[i,1] S : ^^ totally wrong for my needs^^ this if instead : if W mydata[i,0] E and S mydata[i,1] N:

Re: [Numpy-discussion] least() of datetime objects doesn't work

2010-09-11 Thread dyamins
P Sent via BlackBerry from T-Mobile -Original Message- From: Charles R Harris charlesr.har...@gmail.com Sender: numpy-discussion-boun...@scipy.org Date: Thu, 19 Aug 2010 18:03:29 To: Discussion of Numerical Pythonnumpy-discussion@scipy.org Reply-To: Discussion of Numerical Python

Re: [Numpy-discussion] scan array to extract min-max values (with if condition)

2010-09-11 Thread Brett Olsen
On Sat, Sep 11, 2010 at 4:46 PM, Massimo Di Stefano massimodisa...@gmail.com wrote: Thanks Pierre, i tried it and all works fine and fast. my apologize :-( i used a wrong if statment to represent my needs if mydata[i,0] E or mydata[i,0] W or mydata[i,1] N or mydata[i,1] S : ^^