wolffd@0: function inds = som_sub2ind(msize,Subs)
wolffd@0: 
wolffd@0: %SOM_SUB2IND Linear index from map grid subscripts.
wolffd@0: %
wolffd@0: % ind = som_sub2ind(msize,Subs)
wolffd@0: %
wolffd@0: %  ind = som_sub2ind([10 15],[4 5]);
wolffd@0: %  ind = som_sub2ind(sMap,[4 5]);
wolffd@0: %  ind = som_sub2ind(sMap.msize,[4 5]);
wolffd@0: %  inds = som_sub2ind([10 15],[4 5; 3 2; 1 10]);
wolffd@0: %
wolffd@0: %  Input and output arguments: 
wolffd@0: %   msize  (struct) map or topology struct
wolffd@0: %          (vector) size 1 x m, specifies the map grid size
wolffd@0: %   Subs   (matrix) size n x m, the subscripts of n vectors
wolffd@0: %
wolffd@0: %   inds   (vector) size n x 1, corresponding linear indeces
wolffd@0: %
wolffd@0: % See also SOM_IND2SUB.
wolffd@0: 
wolffd@0: % Contributed to SOM Toolbox vs2, February 2nd, 2000 by Juha Vesanto
wolffd@0: % http://www.cis.hut.fi/projects/somtoolbox/
wolffd@0: 
wolffd@0: % juuso 300798
wolffd@0: 
wolffd@0: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
wolffd@0: 
wolffd@0: if isstruct(msize), 
wolffd@0:   if strcmp(msize.type,'som_map'), msize = msize.topol.msize; 
wolffd@0:   elseif strcmp(msize.type,'som_topol'), msize = msize.msize;
wolffd@0:   else error('Invalid first argument.'); end
wolffd@0: end
wolffd@0: 
wolffd@0: % check off-limits
wolffd@0: [n d] = size(Subs);
wolffd@0: offl = find(Subs < 1 | Subs > msize(ones(n,1),1:d)); 
wolffd@0: Subs(offl) = NaN;
wolffd@0: 
wolffd@0: % indexes
wolffd@0: k = [1 cumprod(msize(1:end-1))]';
wolffd@0: inds = 1 + (Subs-1)*k;
wolffd@0: 
wolffd@0: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%