![](https://secure.gravatar.com/avatar/3659900a66bb9cfc6b575b1d106db079.jpg?s=120&d=mm&r=g)
On Sun, 18 Nov 2007 00:36:06 +0100 Gael Varoquaux <gael.varoquaux@normalesup.org> wrote:
On Sat, Nov 17, 2007 at 05:30:18PM -0600, Alan Jackson wrote:
I'm not sure where the blitz errors came from - the same code compiled clean with inline, so that is a mystery (I did modify the arrays for blitz before compiling).
Just checking: you did replace the "M[x, y]" by "M(x, y)" in all your code?
Yep...
I also created all the arrays I wanted returned in numpy first - which makes for a really long function call, but it works and it is not too awkward.
I always do that. It removes all memory management from the C level. Tht's one less hurdle. I usual wrap my function in a another one, to hide this ugly function call, though.
I got a speed up of x37 on my real test.
Nice !
Gaƫl _______________________________________________ SciPy-user mailing list SciPy-user@scipy.org http://projects.scipy.org/mailman/listinfo/scipy-user
-- ----------------------------------------------------------------------- | Alan K. Jackson | To see a World in a Grain of Sand | | alan@ajackson.org | And a Heaven in a Wild Flower, | | www.ajackson.org | Hold Infinity in the palm of your hand | | Houston, Texas | And Eternity in an hour. - Blake | -----------------------------------------------------------------------