[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

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

From: Mike Miller
Subject: RE: Dot product (nx1).*(nxm)
Date: Sat, 4 Dec 1999 22:27:45 -0600 (CST)

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

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



Michael B. Miller
University of Missouri--Columbia

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

Octave's home on the web:
How to fund new projects:
Subscription information:

reply via email to

[Prev in Thread] Current Thread [Next in Thread]