[Numpy-discussion] efficient 3d histogram creation

Dag Sverre Seljebotn dagss@student.matnat.uio...
Thu May 7 06:35:13 CDT 2009

Dag Sverre Seljebotn wrote:
> Stéfan van der Walt wrote:
>> 2009/5/7 Chris Colbert <sccolbert@gmail.com>:
>>> This was really my first attempt at doing anything constructive with Cython.
>>> It was actually unbelievably easy to work with. I think i spent less time
>>> working on this, than I did trying to find an optimized solution using pure
>>> numpy and python.
>> One aspect we often overlook is how easy it is to write a for-loop in
>> comparison to vectorisation.  Besides, for-loops are sometimes easier
>> to read as well!
>> I think the Cython guys are planning some sort of templating, but I'll
>> CC Dag so that he can tell us more.
> We were discussing how it would/should look like, but noone's committed 
> to implementing it so it's pretty much up in the blue I think -- someone 
> might jump in and do it next week, or it might go another year, I can't 
> tell.

BTW the consensus pretty much ended on:

cdef class MyClass[T](Ancestor):
    cdef T evaluate(T x): ...

And then instantiate with

cdef MyClass[int] obj = MyClass[int]()

Only class templates would be targeted at first.

Dag Sverre

More information about the Numpy-discussion mailing list