Re: [Numpy-discussion] Meta: too many numerical libraries doing the same thing?

"Paul F. Dubois" wrote:
Personally, I think this could be a big gain. At the moment, if you don't get the performance you need with NumPy, you have to write some of your code in C, and using the Numeric and Python C API is a whole lot of work, particularly if you want your function to work on non-contiguous arrays and/or arrays of any type. I don't know much C++, and I have no idea if Blitz++ fits this bill, but it seemed to me that using an object oriented framework that could take care of reference counting, and allow you to work with generic arrays, and index them naturally, etc, would be a great improvement, even if the performance was the same as the current C API. Perhaps NumPy2 has accomplished that, it sounds like it is a step in the right direction, at least. In a sentence: the most important reason for using a C++ object oriented multi-dimensional array package would be easy of use, not speed. It's nice to hear Blitz++ was considered, it was proably rejected for good reason, but it just looked very promising to me. -Chris -- Christopher Barker, Ph.D. ChrisHBarker@home.net --- --- --- http://members.home.net/barkerlohmann ---@@ -----@@ -----@@ ------@@@ ------@@@ ------@@@ Oil Spill Modeling ------ @ ------ @ ------ @ Water Resources Engineering ------- --------- -------- Coastal and Fluvial Hydrodynamics -------------------------------------- ------------------------------------------------------------------------

Travis Oliphant wrote:
Yes, it does. That's where I heard about it. That also brings up a good point. Paul mentioned that using something like Blitz++ would only help performance if you could pass it an entire expression, like: x = a+b+c+d. That is exactly what Eric's compiler module does, and it would sure be easier if NumPy already used Blitz++! In Fact, I suppose Eric's compiler is a start towards a tool that could comp9le en entire NumPy function or module. I'd love to be able to just do that (with some tweeking perhaps) rather than having to code it all by hand. My fantasies continue... -Chris -- Christopher Barker, Ph.D. ChrisHBarker@home.net --- --- --- http://members.home.net/barkerlohmann ---@@ -----@@ -----@@ ------@@@ ------@@@ ------@@@ Oil Spill Modeling ------ @ ------ @ ------ @ Water Resources Engineering ------- --------- -------- Coastal and Fluvial Hydrodynamics -------------------------------------- ------------------------------------------------------------------------

Travis Oliphant wrote:
Yes, it does. That's where I heard about it. That also brings up a good point. Paul mentioned that using something like Blitz++ would only help performance if you could pass it an entire expression, like: x = a+b+c+d. That is exactly what Eric's compiler module does, and it would sure be easier if NumPy already used Blitz++! In Fact, I suppose Eric's compiler is a start towards a tool that could comp9le en entire NumPy function or module. I'd love to be able to just do that (with some tweeking perhaps) rather than having to code it all by hand. My fantasies continue... -Chris -- Christopher Barker, Ph.D. ChrisHBarker@home.net --- --- --- http://members.home.net/barkerlohmann ---@@ -----@@ -----@@ ------@@@ ------@@@ ------@@@ Oil Spill Modeling ------ @ ------ @ ------ @ Water Resources Engineering ------- --------- -------- Coastal and Fluvial Hydrodynamics -------------------------------------- ------------------------------------------------------------------------
participants (2)
-
Chris Barker
-
Travis Oliphant