RE: Dot product (nx1).*(nxm)

[Mike Miller]

On Sat, 4 Dec 1999, Heber Farnsworth wrote:

> p(:,ones(1,m)).*v
> 
> where p is nx1 and v is nxm.  What this does is to multiply v by a nxm
> natrix which has p as every column.


Apparently, Heber's method is faster in this example than the one I
suggested:

octave:1> p=rand(100,1);
octave:2> q=ones(1,50);
octave:3> p2=p*q;
octave:4> t0=cputime; for i=1:1000, p2=p*q; end , cputime-t0
ans = 2.4600
octave:5> t0=cputime; for i=1:1000, p2=p(:,q); end , cputime-t0
ans = 2.0100


Regards,

Mike

-- 
Michael B. Miller
University of Missouri--Columbia
http://taxa.psyc.missouri.edu/~mbmiller/



-----------------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.che.wisc.edu/octave/octave.html
How to fund new projects:  http://www.che.wisc.edu/octave/funding.html
Subscription information:  http://www.che.wisc.edu/octave/archive.html
-----------------------------------------------------------------------