Re: [Numpy-discussion] Further matrix oddities: no inner product

2007-03-29 Thread Sven Schreiber
Charles R Harris schrieb: > > Nope, vdot doesn't work for row and column vectors. So there is *no* > builtin inner product that works for matrices. I wonder if we should > have one, and if so, what it should be called. I think that vdot should > probably be modified to do the job. There is also

[Numpy-discussion] Further matrix oddities: no inner product

2007-03-27 Thread Charles R Harris
In mathematics, and inner product is a sesquilinear form on pairs of vectors, so at the least it should return a scalar. In numpy inner is a sum over the last indices. OK, so we have In [10]: inner(ones(2),ones(2)) Out[10]: 2.0 This doesn't work as an inner product for column vectors, which woul