Re: Effectively selecting elements of a matrix

[Jonathan C. Webster]

Thomas Weber wrote:
> Hello, 
>
> my problem is as follows: 
> I have a matrix A (size: 50x500) and two vectors (vec1, vec2) of length
> 32000. The elements of the vectors are indices for the matrix, i.e. if
> vec1(4) = 5 and vec2(4) = 10, I want to select A(5,10) into a new
> vector results_vec at position 4:
> results_vec (4) = A(5,10)
>  
> My problem is how to create this new vector effectively:
Ravel the matrix A out into one long vector, A_vec, using fortran indexing if 
you want.
Then combine vec1 and vec2 into a single index vector into your A_vec.
 .
vec3 = vec1 + size(A)(1) * (vec2 - 1)

Finally, your new reordered vector A3_vec

A3_vec = A_vec(vec3)

I think thats right, but check it.

Jonathan






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

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