Starting a new thread for this. On Tue, Dec 22, 2009 at 7:13 PM, Anne Archibald <peridot.face...@gmail.com> wrote:
> I think we have one major lacuna: vectorized linear algebra. If I have > to solve a whole whack of four-dimensional linear systems, right now I > need to either write a python loop and use linear algebra on them one > by one, or implement my own linear algebra. It's a frustrating lacuna, > because all the machinery is there: generalized ufuncs and LAPACK > wrappers. Somebody just needs to glue them together. I've even tried > making a start on it, but numpy's ufunc machinery and generic type > system is just too much of a pain for me to make any progress as is. Please be more specific: what (which aspects) have been "too much of a pain"? (I ask out of ignorance, not out of challenging your opinion/experience.) > I think if someone wanted to start building a low-level Again, please be more specific: what do you mean by this? (I know generally what is meant by "low level," but I'd like you to spell out a little more fully what you mean by this in this context.) > generalized > ufunc library interface to LAPACK, that would be a big improvement in > numpy/scipy's linear algebra. Pretty much everything else strikes me > as a question of notation. (Not to trivialize it: good notation makes > a tremendous difference.) Thanks, Anne. DG > > Anne > _______________________________________________ > NumPy-Discussion mailing list > NumPy-Discussion@scipy.org > http://mail.scipy.org/mailman/listinfo/numpy-discussion > _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion