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

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 t

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 m

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

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