Thread Subject: How to do vectorize this simple code?

I want to vectorize this simple snippet:

for i=1:num_elements
    y(i)=my_matrix(row(i),col(i));
end

I tried y=my_matrix(row,col)
and y=my_matrix([row col])

but nothing seems to work. Any pointers? Thanks.

"Henry " <henry@fullspectrumengineering.com> wrote in message <gfkstu$1id$1@fred.mathworks.com>...
> I want to vectorize this simple snippet:
>
> for i=1:num_elements
> y(i)=my_matrix(row(i),col(i));
> end
>
> I tried y=my_matrix(row,col)
> and y=my_matrix([row col])
>
> but nothing seems to work. Any pointers? Thanks.

Check out the SUB2IND function.

hth,
Ken

If I understand you correctly, all you need is the diagonal. I believe this is the simplest way (at least that I'm aware of).

y(i)=diag(my_matrix(row(i),col(i)));

"Scott " <crazyivan84__remove__@hotmail.com> wrote in message <gfktep$9al$1@fred.mathworks.com>...
> If I understand you correctly, all you need is the diagonal. I believe this is the simplest way (at least that I'm aware of).
>
> y(i)=diag(my_matrix(row(i),col(i)));

This is an extremely inefficient way to do it of course.

John

"Henry " <henry@fullspectrumengineering.com> wrote in message <gfkstu$1id$1@fred.mathworks.com>...
> I want to vectorize this simple snippet:
>
> for i=1:num_elements
> y(i)=my_matrix(row(i),col(i));
> end
>
> I tried y=my_matrix(row,col)
> and y=my_matrix([row col])
>
> but nothing seems to work. Any pointers? Thanks.

Hi Henry

I guess this works:

U = mymatrix(row,:);
V = U(:,col);
y = diag(V);

-Mah

Hmm all those suggestions don't seem to be very efficient compared to my original code of:
for i=1:num_elements
    y(i)=my_matrix(row(i),col(i));
end

I just thought I could write it as a single line like y=my_matrix(row,col) where row and col are vectors and it would match them up as pairs and return a vector.

Seems like a common procedure and I'm surprised matlab doesn't have it built in...

"Mahdieh" <mahdieh.emrani@capitalhealth.ca> wrote in message <gfl0tf$srs$1@fred.mathworks.com>...
> "Henry " <henry@fullspectrumengineering.com> wrote in message <gfkstu$1id$1@fred.mathworks.com>...
> > I want to vectorize this simple snippet:
> >
> > for i=1:num_elements
> > y(i)=my_matrix(row(i),col(i));
> > end
> >
> > I tried y=my_matrix(row,col)
> > and y=my_matrix([row col])
> >
> > but nothing seems to work. Any pointers? Thanks.
>
> Hi Henry
>
> I guess this works:
>
> U = mymatrix(row,:);
> V = U(:,col);
> y = diag(V);
>
> -Mah

"John D'Errico" <woodchips@rochester.rr.com> wrote in message <gfkuen$m7v$1@fred.mathworks.com>...
> "Scott " <crazyivan84__remove__@hotmail.com> wrote in message <gfktep$9al$1@fred.mathworks.com>...
> > If I understand you correctly, all you need is the diagonal. I believe this is the simplest way (at least that I'm aware of).
> >
> > y(i)=diag(my_matrix(row(i),col(i)));
>
> This is an extremely inefficient way to do it of course.
>
> John

Actually, I've always found sub2ind to be quite slow.
Consider:

rows=randint(5,100000,[1,10]);
cols=randint(5,100000,[1,10]);
testMat=magic(10);

%Sub2ind
tic
resultsInd=zeros(100000,5);
for i=1:100000
    index=sub2ind(size(testMat),rows(:,i),cols(:,i));
    resultsInd(i,:)=testMat(index);
end
toc

%Diag
tic
resultsDiag=zeros(100000,5);
for i=1:100000
    resultsDiag(i,:)=diag(testMat(rows(:,i),cols(:,i)));
end
toc

isequal(resultsInd,resultsDiag)


With the result:
Elapsed time is 6.405682 seconds. (sub2ind)
Elapsed time is 0.293780 seconds. (diag)

ans =

     1

Your comparison is invalid.

"Henry " <henry@fullspectrumengineering.com> wrote in message <gfl212$m0t$1@fred.mathworks.com>...
> Hmm all those suggestions don't seem to be very efficient compared to my original code of:
> for i=1:num_elements
> y(i)=my_matrix(row(i),col(i));
> end
>
> I just thought I could write it as a single line like y=my_matrix(row,col) where row and col are vectors and it would match them up as pairs and return a vector.
>
> Seems like a common procedure and I'm surprised matlab doesn't have it built in...

I don't know if SUB2IND is vectorized. If not, try this:

N = size(my_matrix,1);
ind = N*(col-1)+row;
y = my_matrix(size(my_matrix,1)*(col-1)+row);

hth,
Ken

"Kenneth Eaton" <Kenneth.dot.Eaton@cchmc.dot.org> wrote in message <gflhc9$h1d$1@fred.mathworks.com>...
>
> I don't know if SUB2IND is vectorized. If not, try this:
>
> N = size(my_matrix,1);
> ind = N*(col-1)+row;
> y = my_matrix(size(my_matrix,1)*(col-1)+row);

Oops, mistyped that a little; you can remove the first 2 lines, using the third as a one-liner, or make the third line: y = my_matrix(ind);

Ken

"Henry " <henry@fullspectrumengineering.com> wrote in message <gfkstu$1id$1@fred.mathworks.com>...
> I want to vectorize this simple snippet:
>
> for i=1:num_elements
> y(i)=my_matrix(row(i),col(i));
> end
>
> I tried y=my_matrix(row,col)
> and y=my_matrix([row col])
>
> but nothing seems to work. Any pointers? Thanks.

This works, but its not so obvious and simple as it should be

M=rand(4)
row=(1:3);
col=(2:4);
y=diag(M(row,col))
/Per