I recently note that there is no C++ array template library shipped with Sage. I think this is an impediment for more serious numerical computations. The whole Cython thing is really nice to speed up Python code, but if you are serious about avoiding heap allocation of array entries then Cython becomes less elegant. So I would like to start some discussion about this ;-)
One possibility is blitz++ (http://www.oonumerics.org/blitz/). Others would be Eigen, boost::ublas, or boost::mtl, I made a tentative spkg here: http://www.stp.dias.ie/~vbraun/Sage/spkg/blitz-20120201.spkg To really fit into Sage one would need a Cython wrapper to access it from Python. Right now, Cython doesn't support all the necessary C++ features so its not a straightforward task. But in the longer run either Cython will become better with C++ or we'll work around it, so its not impossible. Advantages: * blitz++ has been around for a while and is essentially feature-complete * scipy.weave.blitz would benefit from it Disadvantages: * Mystical compiler error messages Comments? -- To post to this group, send an email to sage-devel@googlegroups.com To unsubscribe from this group, send an email to sage-devel+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URL: http://www.sagemath.org
