Re: [Numpy-discussion] doing zillions of 3x3 determinants fast

2008-08-25 Thread Dan Lenski
On Sun, 24 Aug 2008 23:49:45 -0600, Charles R Harris wrote: On Sun, Aug 24, 2008 at 9:48 PM, Daniel Lenski [EMAIL PROTECTED] wrote: Hi all, I need to take the determinants of a large number of 3x3 matrices, in order to determine for each of N points, in which of M tetrahedral cells they

[Numpy-discussion] doing zillions of 3x3 determinants fast

2008-08-24 Thread Daniel Lenski
Hi all, I need to take the determinants of a large number of 3x3 matrices, in order to determine for each of N points, in which of M tetrahedral cells they lie. I arrange the matrices in an ndarray of shape (N,M,5,3,3). As far as I can tell, Numpy doesn't have a function to do determinants

Re: [Numpy-discussion] doing zillions of 3x3 determinants fast

2008-08-24 Thread Daniel Lenski
On Mon, 25 Aug 2008 03:48:54 +, Daniel Lenski wrote: * it's fast enough for 100,000 determinants, but it bogs due to all the temporary arrays when I try to do 1,000,000 determinants (=72 MB array) I've managed to reduce the memory usage significantly by getting the number of

Re: [Numpy-discussion] doing zillions of 3x3 determinants fast

2008-08-24 Thread Anne Archibald
2008/8/25 Daniel Lenski [EMAIL PROTECTED]: On Mon, 25 Aug 2008 03:48:54 +, Daniel Lenski wrote: * it's fast enough for 100,000 determinants, but it bogs due to all the temporary arrays when I try to do 1,000,000 determinants (=72 MB array) I've managed to reduce the memory

Re: [Numpy-discussion] doing zillions of 3x3 determinants fast

2008-08-24 Thread Charles R Harris
On Sun, Aug 24, 2008 at 9:48 PM, Daniel Lenski [EMAIL PROTECTED] wrote: Hi all, I need to take the determinants of a large number of 3x3 matrices, in order to determine for each of N points, in which of M tetrahedral cells they lie. I arrange the matrices in an ndarray of shape (N,M,5,3,3).