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 can

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 UPDATEIFCOPY : False

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.

[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 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 would

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

2009-03-16 Thread Charles R Harris
2009/3/16 Ryan May rma...@gmail.com 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

[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] 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)

Re: [Numpy-discussion] Overlapping ranges

2009-03-16 Thread Robert Kern
2009/3/16 Peter Saffrey p...@dcs.gla.ac.uk: 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

Re: [Numpy-discussion] Overlapping ranges

2009-03-16 Thread josef . pktd
On Mon, Mar 16, 2009 at 5:29 PM, Robert Kern robert.k...@gmail.com wrote: 2009/3/16 Peter Saffrey p...@dcs.gla.ac.uk: 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