Re: [Numpy-discussion] Different results from repeated calculation, part 2

2008-08-14 Thread Holger Rapp
Hi, Am 14.08.2008 um 19:48 schrieb Alok Singhal: > On 14/08/08: 10:20, Keith Goodman wrote: >> A unit test is attached. It contains three tests: >> >> In test1, I construct matrices x and y and then repeatedly >> calculate z >> = calc(x,y). The result z is the same every time. So this test >>

Re: [Numpy-discussion] SSEPlus + Framewave

2008-08-13 Thread Holger Rapp
> >> Or do it get you completly wrong? Is your intention to make a plugin >> architecture in the sense of: copy some directory with libs and >> config >> in your site-packages and then your multiplications are much >> faster? I >> would consider such a framework a bit overengineered, since spee

Re: [Numpy-discussion] SSEPlus + Framewave

2008-08-13 Thread Holger Rapp
Hello David, > The problem is not so much the build part, but the clear separation I > was talking about. My experience with ATLAS convinced me the only way > to make sse work reliably is to detect the CPU arch at runtime; > compiling binaries incompatible on different arch is just not scalable >

[Numpy-discussion] SSEPlus + Framewave

2008-08-11 Thread Holger Rapp
Hi everyone, I have a question concerning performance of numpy. I'm using it for heavy image processing tasks and often need more speed than the stock numpy delivers. Especially in numeric calculations with big arrays (in my current case shape = (8,4,1200,1600), 'float'). So I often rely on