[Numpy-discussion] LA improvements (was: dot function or dot notation, matrices, arrays?)

[David Goldsmith]

Starting a new thread for this.

On Tue, Dec 22, 2009 at 7:13 PM, Anne Archibald
<peridot.faceted at 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 at scipy.org
> http://mail.scipy.org/mailman/listinfo/numpy-discussion
>