Re: [Numpy-discussion] Overlapping ranges

2009-03-16 Thread josef . pktd
On Mon, Mar 16, 2009 at 5:29 PM, Robert Kern wrote: > 2009/3/16 Peter Saffrey : > >> At the moment, I'm using a fairly naive approach that finds roughly in the >> genome (which gene) each point might be and then checking it against the >> bins in that gene. If I split the problem into chromosomes,

Re: [Numpy-discussion] Overlapping ranges

2009-03-16 Thread Robert Kern
2009/3/16 Peter Saffrey : > At the moment, I'm using a fairly naive approach that finds roughly in the > genome (which gene) each point might be and then checking it against the > bins in that gene. If I split the problem into chromosomes, I feel sure > there must be some super-fast matrix approac

[Numpy-discussion] Overlapping ranges

2009-03-16 Thread Peter Saffrey
I'm trying to file a set of data points, defined by genome coordinates, into bins, also based on genome coordinates. Each data point is (chromosome, start, end, point) and each bin is (chromosome, start, end). I have about 140 million points to file into around 100,000 bins. Both are (roughly)

[Numpy-discussion] 1.3.x branch created - trunk now opened for 1.4

2009-03-16 Thread David Cournapeau
Hi, I have just started the 1.3.x branch - as such, any change done to the trunk will not end up in the 1.3 release. I will announce the 1.3 beta release within the day, hopefully, cheers, David ___ Numpy-discussion mailing list Numpy-discussion@scipy.

Re: [Numpy-discussion] svn and tickets email status

2009-03-16 Thread Charles R Harris
2009/3/16 Ryan May > Hi, > > What's the status on SVN and ticket email notifications? The only messages > I'm seeing since the switch is the occasional spam. Should I try > re-subscribing? > I get the ticket notifications but I think the svn notifications are still broken. I needed to update m

Re: [Numpy-discussion] Superfluous array transpose (cf. ticket #1054)

2009-03-16 Thread Pearu Peterson
On Mon, March 16, 2009 4:05 pm, Sturla Molden wrote: > On 3/16/2009 9:27 AM, Pearu Peterson wrote: > >> If a operation produces new array then the new array should have the >> storage properties of the lhs operand. > > That would not be enough, as 1+a would behave differently from a+1. The > former

[Numpy-discussion] svn and tickets email status

2009-03-16 Thread Ryan May
Hi, What's the status on SVN and ticket email notifications? The only messages I'm seeing since the switch is the occasional spam. Should I try re-subscribing? Ryan -- Ryan May Graduate Research Assistant School of Meteorology University of Oklahoma Sent from: Norman Oklahoma United States. _

Re: [Numpy-discussion] Superfluous array transpose (cf. ticket #1054)

2009-03-16 Thread Sturla Molden
On 3/16/2009 9:27 AM, Pearu Peterson wrote: > If a operation produces new array then the new array should have the > storage properties of the lhs operand. That would not be enough, as 1+a would behave differently from a+1. The former would change storage order and the latter would not. Broadca

Re: [Numpy-discussion] Superfluous array transpose (cf. ticket #1054)

2009-03-16 Thread Pearu Peterson
On Sun, March 15, 2009 8:57 pm, Sturla Molden wrote: > > Regarding ticket #1054. What is the reason for this strange behaviour? > a = np.zeros((10,10),order='F') a.flags > C_CONTIGUOUS : False > F_CONTIGUOUS : True > OWNDATA : True > WRITEABLE : True > ALIGNED : True > UPDATEI

Re: [Numpy-discussion] inplace dot products

2009-03-16 Thread David Warde-Farley
On 20-Feb-09, at 6:41 AM, Olivier Grisel wrote: > Alright, thanks for the reply. > > Is there a canonical way /sample code to gain low level access to > blas / lapack > atlas routines using ctypes from numpy / scipy code? > > I don't mind fixing the dimensions and the ndtype of my array if it